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BSTRACT. A simple method of photographing Fabry and Perot interferometer fringes 
outlined, and some examples of the results obtained with mercury radiation are given. 


FTER some experiments on the photography of the fringe systems produced 
by a Fabry and Perot interferometer it was thought that a few examples 
such as those which accompany this note, while possessing no originality, 
ight not be devoid of interest as roughly indicating the results to be obtained by 
e use of a simple optical system. The only auxiliary apparatus employed, together 
ith the necessary camera and source of light, was a convergent lens. The arrange- 
ent adopted is shown diagrammatically in Fig. 1. 


Fig. 1. Arrangement of apparatus. 


- The condensing lens L caused light from the source SS to pass through the inter- 
rometer plates FP to form an image of the source approximately at the front of 
e camera lens C, the fringe system being observed on the receiving screen N. A 
nvenient arrangement was obtained by placing the camera lens and the source of 
ght at distances of about 40 cm. and 80 cm. respectively from the condenser L, 
hich was a single uncorrected lens of approximately 25 cm. mean focal length. 
his lens was of large aperture (14°5 cm.) and thus possessed the dual advantage of 
llecting a considerable amount of light from the source while subtending a large 
lid angle at the camera lens. The interferometer plates were placed as close to the 


mera lens as the interferometer mounting would allow (about 6 cm. away). By 
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this means the plates, which had an aperture of 2-6 cm., were prevented from cutting 
off an appreciable fraction of the cone of light converging on to the camera lens 
from the condenser. Hence the angle of the cone of illumination falling on to the 
receiving screen (shown approximately by the dotted lines in the diagram) was 
maintained fairly large. The visible portion of the mercury vapour lamp employed 
as the source of light being rectangular (3-5 x 2 cm.), the condensing lens L was 
turned until the patch of light on the nearer interferometer plate was approximately 
square and just about covered the plate so that a minimum quantity of light was 
lost. 

With the optical system arranged approximately as described, the quarter-plate 
screen of the camera could be almost entirely covered with fringes. The camera lens — 
possessed a focal length of 29 cm. and for the photography was used with the stop ; 
at j/I1. 

Hi use of a medium aperture condensing lens of short focal length, placed 
closer to the interferometer plates to obtain good covering power on the receiving 
screen, was prevented by the presence of the interferometer bed which carried the 
movable plate. This bed extended for a distance of nearly 20 cm. beyond the plates — 
towards the condensing lens. All the apparatus rested on shelves fixed to the 
laboratory wall, the interferometer being on a separate shelf. Mechanical insulation | 
of the interferometer was further improved by the insertion of a thick layer of — 
cotton wool underneath the baseboard of the instrument. 

Ilford Special Rapid Panchromatic plates were employed in taking the photo- 
graphs shown in Fig. 2. The details were as follows: 


Table. Details relating to Fig. 2 


Approximate plate ' Exposure 


Photograph Fringe system separation (mm.) | (sec.) 


(a) Full mercury as o-7 « os 
(b) Full mercury 2°0 bare) 
(c) Mercury green o-9 75 
5461 
(d) Mercury green 2°8 75 
A 5461 


Angstrom unit from the main line*. Calculation shows that with the plates approxi- 
mately 3 mm. apart such a satellite should give fringes about equally spaced between 
those of the principal system. The edges of the main fringes on the sides of higher 
wave-length are diffuse. These results are in very fair agreement with the measure- 
ments given by H. G. Gale and H. B. Lemon in a discussion of the fine structure 
of the salient mercury linest. In their paper a grating photograph of A 5461, 
obtained from a Cooper-Hewitt lamp, is given. 

A filter to isolate this radiation was conveniently made by mixture of dilute 
solutions of malachite green and potassium bichromate. The glass cell containing 

* This satellite is on the side of lower wave-length, 

t Astrophysical Journal, 31, 78 (1910). See also J.C. McLennan, Proc. R. S. A. 87, 276 (1912). 


| 
The satellite shown in Fig. 2 (d) appears to be situated at about a quarter of an | 
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ey Fig. 2(a). (Original width 6-2 cm.) 
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Fig. 2(b). (Original width 5:8 cm.) 


| 
| 
| 
| PHYSICAL SOCIETY, VOL. 42, PT 3 (GREEN) PLATE III 
| 
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Fig. 2(c). (Original width 5-8 cm.) 


PLATE IV 
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the liquid was placed close to the interferometer plates on the side towards the 
incident light. The manner in which such a filter may be adjusted to isolate the 
desired green line is approximately as follows. A fresh dilute solution of malachite 


_green crystals absorbs the orange and red portions of the spectrum roughly between 


the wave-lengths 6400 and 6000 A.u. By increase of the concentration this absorp- 
tion band may be made to extend into the red at 6700 and the green at 5100. The 
edge of the band in the red is sharp and varies in position with the solution-strength 
much less than does the more diffuse edge of lower wave-length. When the con- 


_ centration of the malachite green is adjusted so that the range of elimination stretches 


from about 6600 to 5500 the strong pair of yellow lines \A 5791, 5770 are cut out 
while the required green radiation is transmitted. A solution of this strength also 
exhibits absorption in the violet, and in the blue as far as about A 4500. In this 
region, however, the absorption is not strong enough to remove entirely the strong 
violet line A 4358, so that the addition of a little potassium bichromate is necessary 
both to complete the absorption in the violet and to carry the elimination sufficiently 


_ far into the blue to remove the green-blue radiation A 4960. Incidentally this filter 
- cuts out a group of orange-red lines in the neighbourhood of 6200. In the case of 


the mercury lamp used in the present experiments these radiations were relatively 
strong only when the lamp had just been switched on when cold. Otherwise they 
were practically invisible. 

Figures 2(a) and 2(5) exhibit some lack of uniformity in illumination. In 
these cases, the necessary exposures being short, this uniformity could have been 
improved by placing a circular stop 2 cm. in diameter immediately in front of the 
lamp and keeping the condenser L perpendicular to the optic axis of the system. 
The resulting diminution in the amount of light falling on the receiving screen 
would not have mattered greatly, being corrected for by longer exposure. ‘The 
much more extended exposures necessary owing to the presence of the filter when 
the green fringe system alone was being photographed made the arrangement 
actually adopted preferable for figures (c) and (d), since thereby little light was lost. 

I should like to take the opportunity of expressing my gratitude to Dr L. J. 
Freeman, in collaboration with whom some earlier experiments on the subject were 


_ performed, for his very extensive assistance and interest and for being so kind as 


to lend me the camera with which these photographs were taken. 


DISCUSSION 


Dr J. S. ANDERSON. I think the author is to be congratulated on obtaining 
such beautiful interferograms with simple apparatus. It would enhance the value 
of the paper if he would add information regarding the type of mercury vapour 
lamp used and its current consumption, and also the degree of enlargement of the 
photographs from the originals. Has the author tried any other type of plate? 
In work with the Hilger universal lens interferometer I have obtained satisfactory 
photographs with exposures of about 15 sec., using Ilford rapid chromatic plates 


and a Wratten and Wainwright mercury green filter, and overrunning the mercury 
i Nf 
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vapour lamp (a quartz one). In photographing Fabry and Perot fringes the con- 
ditions are somewhat different and it would not, I presume, be advisable to overrun 
the lamp. For my work it would be of great advantage to obtain a more sensitive 
plate, but the grain must be of small size as the original photographs are small and 
require to be enlarged. 


Dr L. C. Martin. It is usually the experience that a quartz mercury vapour 
lamp is unsuitable for interferometric work on account of the internal pressure 
developed when the lamp grows hot. I take it that Mr Green used a Cooper- 
Hewitt lamp for such beautiful photographs as those which he has shown us? 


AUTHOR’S REPLY. The mercury lamp used was a Cooper-Hewitt lamp, having a 
glass tube 2 cm. in diameter and carrying a current of 3-3 amp. It was inadvisable 
to obtain a brighter source of light by overrunning the lamp, in view of the broad- 
ening of the spectral lines consequent upon the increase of temperature and pres- 
sure. The preliminary set of photographs taken in earlier experiments and shown 
at the meeting have enlargement factors from the original negatives varying between 
18 and 2:5, while the positives in Fig. 2 have corresponding values 2-2 in the case 
of (a), (b) and (c) and 3:5 in the case of (d). No systematic test was made with other 
types of photographic plates. 

With regard to the effect of the conditions in the source on the interference 
pattern produced, perhaps I may mention that on attempting to analyse the — 
mercury radiation from a vertical capillary arc operating at a pressure of about — 
go cm. of mercury and carrying 0-3 amp. I found that the fringe system in a Fabry — 
and Perot étalon had vanished almost completely, whereas the Cooper-Hewitt lamp 
gave satisfactory fringes. This was presumably due to the excessive widening of the 
spectral lines at the high temperature and pressure of the capillary arc. It is pos- 
sible that the fringes given by the Cooper-Hewitt lamp might have been further — 
sharpened by reduction of the current strength below 3-3 amp., but any advantage — 
gained would have been largely offset by the necessity for longer exposures. | 
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CHARACTERISTICS OF DISCHARGE TUBES UNDER 
“FLASHING” CONDITIONS AS DETERMINED BY 
MEANS OF THE CATHODE RAY OSCILLOGRAPH 


BY Wa Ay GEYSHON, PH.D. 


ABSTRACT. An account is given of a method of determining the current-voltage 
characteristics of discharge tubes under flashing conditions, a cathode ray oscillograph 
being used. Results are given for neon lamps of the commercial Osglim type. Current/ 
time and voltage/time curves are derived from the oscillograph records. It is suggested 
that the method might prove useful in investigations of intermittent discharges in various 
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gases under different conditions of pressure and disposition of the electrodes. 


§r. INTRODUCTION 


which a flashing neon tube circuit has been used for the excitation of nerve 
: preparations in order to produce limb movements of natural form. It was 
thought desirable that some information should be obtained about the nature and 
duration of the condenser discharges which were utilised for this purpose, and 
accordingly a number of observations were made with a cathode ray oscillograph. 


A very brief account of the experimental method and of the results obtained has 
(1) 


Posen have been carried out during the past two or three years in 


been given in a recent paper 

In the course of the experiments results were obtained relating to the character- 
istics of flashing neon tubes which were not relevant to the original purpose of the 
investigation, but seemed to be of general interest. It is these results, not previously 
tecorded, which are set out here. 

Clarkson) in a series of papers has given accounts of experiments on the 
characteristics of discharge tubes in flashing circuits ; the method of experimentation 
and the circuits which he used differed from those described in the present paper. 
For the purpose of the physiological experiments it was necessary to include a 
resistance or an inductance in series with the condenser or the neon lamp in the 
now well-known flashing circuit first described by M. D. Hart and also by Pearson 
and Anson. The characteristic curves given here are therefore those for a neon lamp 
in series with an ohmic resistance, or an inductance, or both, through which a 
condenser discharges at regular intervals. 

Observations have been made on about ten different Osglim lamps, chiefly of 
the “beehive” pattern; the latter were used in the nerve excitation apparatus, since, 
with other circuit conditions the same, larger quantities of electricity passed during 
discharge through this type of lamp than during discharge through other types. 
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§2. DESCRIPTION OF APPARATUS USED 


Fig. 1 at (a) and (6) shows diagrammatically the arrangement of the flashing 
circuit, The neon lamp N, from which the series resistance 11 the cap has been re- 
moved, is in parallel with the condenser C which has the non-inductive resistance r 
in series with it, Fig. 1 (a), or alternatively the resistance 7 Is In series with the neon 
lamp, Fig. 1 (6). The charging circuit consists of a battery of dry cells or high-tension 
accumulators and a resistance R which may vary from 10,000 to 100,000 ohms. 


C C 


R 


(6) 


N neon lamp; R, r non-inductive resistances ; C condenser; B battery. 


(c) 
B, B copper blocks; S copper sulphate solution; M copper wire. 


Fig. 1. The flashing circuit. 
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(a) Beehive lamp N: r=7600 ohms; (b) Lamp B5: r= 1000 ohms; 
C=o-ors xh. C="05 pF. 


Fig. 2. ‘Typical oscillograph figures recorded by tracing. 


In the physiological experiments a non-inductive potential divider of the type 
shown in Fig. 1 (c) has been used in the discharge circuit. This consists of two copper 
blocks, B, B, in a solution of copper sulphate; the copper wire M dipping into the 
solution acts as the movable point. For the purpose of obtaining the oscillograph 


records, this arrangement has occasionally been replaced by a copper sulphate 
liquid resistance of the pattern described by Jordan®), { 


' 
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The potential differences developed across R and r are applied to the two pairs 
of deflecting plates of a cathode ray oscillograph. The figures seen have been re- 
corded either by tracing or photographically using bromide paper. 

It will be observed that in the arrangement of Fig. 1 (a) both the charging and 
the discharge current of the condenser C pass through the resistance 7, whereas 
in the arrangement of Fig. 1 (b) the discharge current only passes through r. 


Fig. 3. Apparatus for obtaining discharge-current/time curves. 
C,, C, condensers; B,, B, batteries; N,, N, neon lamps; 
Ca cathode ray oscillograph; R, r resistances. 


Some typical results are shown in Fig. 2. The parts marked XY are rather faint 
in the actual oscillograph figures owing to the rapidity of travel of the oscillograph 
spot, and are therefore more difficult to trace accurately than the rest of the figure. 
Discharge-current/time curves have also been obtained with a neon-tube diode 
flashing circuit, shown diagrammatically in Fig. 3. The lamp JN, (usually an I-type 
lamp) is in series with a diode D (a triode with anode and grid connected) and 


Fig. 3 (a). Discharge current (circuit resistance 2000 ohms). 


battery B, the condenser C, being in parallel with N,. Since the charging current 
to the condenser is by this means kept constant (the filament temperature of the 


diode being low enough to ensure that the saturation current is taken throughout) 


the potential difference across the condenser (and therefore that across the diode) 
changes linearly with time during the charging period. This latter potential differ- 
ence is applied to one pair of plates of the oscillograph, and affords a time base. 'To 


Milliamperes 
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the other pair of plates are connected the two ends of the resistance 7 in the discharge 
circuit of the second neon lamp N,. The lamp JN, is made to flash at a frequency much 
higher than that of N,. In order to keep the oscillograph figure steady, a small 
coupling is introduced between the two circuits. 

The discharge-current/time curves, recorded photographically on bromide 
paper, are found to be of the same shape as those deduced (see § 3) from the records 
like those of Fig. 2. A typical record is shown in Fig. 3 (a). 

It may be noted that if the voltage Ri, were applied to the oscillograph plates, 
a record of condenser voltage and time could be obtained, since condenser voltage 
— B — Ri,. From this it would be possible to construct the characteristic of the 
lamp when there is no resistance in the discharge circuit. 


§3. RESULTS 


The curves of Fig. 2 show the relation between r (¢—7,) and Ri, Fig. 2(a), or 77 
and Ri,, Fig. 2 (b), where 7 is the discharge current from and 2, the charging current 
to the condenser. 
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(A) Lamp beehive N: C=o-015 nF’; R=70,000 w; ° (B) Lamp J, (I-type): C=0-0022 uF; R=100,0¢ 
B=174:4 volts; (a) r=7600; (b) r=2800w; B= 162 volts; (a) r=3000; (6) r= 1850 
(c) r=880 a, (c) r=1200 w. 


Fig. 4. Characteristics of neon lamp under flashing conditions. 


The voltage across the lamp is (B — Ri,), Fig. 1a, or (B— Ri, — ri), Fig. 1 (d), 
where B is the battery voltage and i the current through the lamp, so that if B, R,r 
and the sensitivity of the oscillograph are known the characteristic of the lamp anda 
flashing conditions can readily be derived from the figures. The latter, in fact, for 
the arrangement of Fig. 1 (a), give directly the general form of the characteristic, 
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Milliamperes 


110 120 130 140 150 160 
Volts 


(A) Lamp B5: B=156; R=100,000w; r=1000w; (a) C=0'003 pF’, n=614; 
(6b) C=0'05 pF, n=72; (c) Static characteristic. 
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(B) Lamp J,: B=162 volts; R= 100,000 w; r=3000 w; (a) C=0'0022 pF, n=750; 
(b) C=0:0222 nF, n=103; (c) Static characteristic, 


Fig. 5. Characteristics of neon lamps: (A) beehive; (B) I-type. 
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The voltage B was measured under working conditions by means of an electro- 
static voltmeter, which was also used in determining the sensitivity of the oscillo- 
graph. From the figures and the circuit constants the relation between condenser 
voltage and resultant current (7 — i,) during the discharge can be found ; the general 
form of the curve is shown in Fig. 2(b). Since (i—%) =C.d V /dt the time relations 
in the circuit can be deduced if C is known. . 

Fig. 4 shows some typical characteristic curves of a beehive lamp, for which the 
circuit conditions are stated, for different values of the discharge circuit resistance. 
Fig. 5 illustrates characteristics of a beehive and an I-type lamp at (A) and (B) 
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Fig. 6. 


(a) I. Discharge-current/time curve for conditions corresponding to curve a, Fig. 5 (A). 
II. Lamp-voltage/time curve. 

(b) I. Discharge-current/time curve for conditions corresponding to curve 6, Fig. 5(B). 
II. Lamp-voltage/time curve ae 


” » 


respectively. For comparison, the static characteristics of the lamps are also shown. 
The arrangement used was that of Fig. 1 (d). Fig. 6 shows the time relations cal- 
culated from the curves of Fig. 5 and the circuit constants. 

The method of calculation is as follows: The curve showing the relation between 
condenser voltage V and current (¢ — 7) is first drawn. Since (i — 4) = C.dV /dt 
we may write At = C.AV’/[i — i,], where [z — ¢,] is the mean value of (i — 4) in the 
interval AV. Starting with the highest value of condenser voltage, small changes of 
voltage AV are marked out on the (é — i,, V) graph, and the corresponding times 
At calculated. The sum of these time-intervals gives the time which has elapsed | 
since the initiation of the discharge ; thus the relatioriship between condenser current 
and time can be determined. From the circuit constants and the figures orginal 
traced corresponding values of condenser current, lamp voltage and discharge — 
current are found, so that the time relations for all the quantities concerned may be 
traced out. The results may be checked by calculation of the quantity of electricity — 
lost by the condenser during discharge (i) from the area of the {(i — i,), t} curve, 
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(ui) from C(V,—V,), where V, and V, are the highest and lowest condenser 
voltages recorded, respectively. 


Current (milliamperes) 


O 5 10 15 20 DS; 
Time x 10-> seconds 


Fig. 7. Typical discharge-current/time graphs for two beehive lamps. 
(a) Lamp B1: B=174 volts; R=60,000w; r=1000w; C=0'05 uF 
(6) Lamp B6: B=192 volts; R=70,000 w; r=1000w; C=0'05 uF 
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Volts 
Fig. 8. Characteristic curve for lamp Bs. 
i R=100,000, 7= 72; C=0:05. 
Lamp Bs J x R= 80,000, n=128; B=156. 
+ R= 50,000, n=192; 41 =1000. 


Fig. 7 illustrates the difference in the current/time relationships for two beehive 
lamps, one, B6, having a cathode consisting of a more open spiral than the other, Br. 
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. isti caused — 
Fig. 8 shows that the difference, if any, in the characteristic of : ae 
by alteration of the charging circuit resistance K is within the a fo) Fol —— 
error in the tracing of the curves. This result is of importance from p 


Milliamperes 
oo 


120 130 140 150 160 170 180 
, Volts 
Fig. 9. Effect of inductance in discharge circuit on characteristic curve. 
Lamp beehive N: C=o:01 pF; B=178 volts; r= 5240; (a) No inductance in discharge 
circuit; (b) Two No. 400 Igranic coils in discharge circuit; R = 70,000 w; (c) Two No. 400 with 
iron core in discharge circuit; R = 60,000 w. 


Discharge current (milliamperes) 


0 2 4 6 8 10° 12° 14-16 18 90 99 94 96°98 


Time x 10-5 seconds 
Fig. 10. Effect of inductance in discharge circuit on (z, t) curve. 


view of the physiological application of the circuit, since it shows that the rate of 
flashing can be altered without materially affecting the quantity of electricity 
passing during each single condenser discharge. 
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: It is clear from the results given, and from a large number of similar observa- 
tions which have been made, that (1) The range of alteration of voltage across the 
lamp is decreased, and the shape of the characteristic is altered, when the condenser 
capacity is raised (see Fig. 5(A) and 5(B)). This range is not affected to any great ex- 

tent by changes in battery voltage or in series charging resistance—at least within the 
limits of these experiments, and provided that the applied voltage be greater than 
the normal striking voltage. (2) With large values of condenser capacity the voltage 
across the lamp remains almost constant during the greater part of the discharge, 
vafter the initial rapid drop. (3) The presence of inductance in the charging circuit 
is found experimentally to lower considerably the critical resistance for flashing. 
This has been previously noted. 

The effect on the lamp characteristic of inductance in the discharge circuit is 
shown in Figs. 9 and 10. With a large inductance the current may rise very sharply 
while a small drop in voltage occurs: i.e. the negative slope resistance may be very 
great. 

§4. CONCLUSION 


The method given for obtaining characteristic curves of neon lamps involves 
the use of no more apparatus than the flashing circuit itself and the cathode ray 
oscillograph. So far it has been applied only to the commercial type of lamp, but 
it is clear that it could readily be applied to the study of intermittent discharges in 
pure neon and other gases, under different conditions of pressure and disposition 

of the electrodes. The use of the cathode ray oscillograph avoids the errors intro- 
duced by other types of recorders, such as that used by Anthouard ‘S) in his work on 
intermittent discharges through air. 
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DISCUSSION 


Dr D.W. Dye. I wish to congratulate the author on the execution of a very 
useful piece of work which has given definite results, while her methods have shown 
promise of usefulness in allied investigations. I have not had much experience of 
discharges of this nature in neon, but some years ago I tried such tubes giving 
intermittent discharges to see if they could be used as sources of discontinuous 
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oscillations from which high-order harmonics could be selected to serve as radio 
frequency standards. I was compelled to discard the type of tubes then available 
because it was found that within one cycle the discharge presented discontinuities 
which were irregular, and so clouded over the precise selection of a harmonic. 
It appears, however, from the author’s diagrams that the cycle is smooth and 
regular so that it is difficult to see how the results I obtained could agree with this. 

Again recently I tried the use of intermittent discharges from neon tubes to 
serve as a stroboscopic source of illumination of interference fringes produced 
between two surfaces, one of which was vibrating rapidly. In this case also no 
sharpness of definition was possible, but whether owing to luminous persistence 
of the discharge or to irregular discontinuities I could not say. 


Dr A. B. Woop. In her introductory remarks Dr Leyshon mentioned that the 
object of the experiments was to determine the wave-form of the oscillations in 
the neon lamp circuit. The voltage/current characteristic, whilst serving this purpose 
indirectly, does not appear to me to afford the most satisfactory method. Of course 
if a non-photographic oscillograph, such as the Western Electric type, is used, this 
method is the only one possible. Some years ago (1924), however, I obtained a 
series of cathode-ray oscillograph records of the oscillations of an Osglim (neon- 
filled) lamp, using the circuit shown in Fig. 11. This is the Anson-Pearson flashing 
circuit in which two condensers C,, C, are used in series to form a potential divider 
for oscillograph purposes. With this arrangement voltage/time curves of the com-— 
mencement of oscillations in the circuit were recorded in a single traverse of the 
cathode spot across a photographic plate in the vacuum chamber of the oscillo- 
graph. Such records are shown in Fig. 12. The frequency of oscillation of the circuit 
in these cases varied from 1000 to 3000~. The following points are of interest: 
(1) The voltage across the lamp at first rises (during the condenser charge) to a 
high striking potential, at which point it suddenly falls, recovers, and continues 
to oscillate. The mean voltage of the sustained oscillation is, as would be expected, 
considerably lower than the striking voltage. (2) The higher the frequency the | 
smaller the voltage-amplitude of the oscillations. Ultimately it was found that the | 
lamp would not oscillate at frequencies of the order of 100,000 ~ in the circuit 
shown in Fig. 1. (3) As Dr Dye has already suggested, the initial rise of voltage, 
before the lamp begins to glow, is not always regular. Some of the records show 
distinct irregularities in the rising voltage. (4) The wave form is typical of B. yan 
der Pol’s “relaxation” oscillations—a succession of slow rises and sudden falls 
of voltage. I cannot understand why the previous speaker failed to observe the 
presence of harmonics in such a wave-form as that indicated in the records. 


The PrestDENT. The use of the cathode ray oscillograph in examining inter- 
mittent discharges might help in the elucidation of some of the problems of — 
conduction of electricity through gases left unsolved by the more orthodox modes _ 
of experiment. I understand that Dr Leyshon is proceeding with these investiga-_ 


tions but employing simpler electrode systems, more amenable to calculation than 
those in the commercial forms of neon lamp. 
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Fig. rr. Anson-Pearson flashing circuit—oscillograph recording voltage/time characteristic. 
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Mr L. G. CarPenter. In view of the growing importance of the cathode-ray 
oscillograph as a means of investigating transient phenomena, it would enhance 
the value of the author’s very interesting paper if she would add some details as 
to the characteristics of the cathode-ray oscillograph she used. Any information 
as to the best material for the fluorescent screen in visual observations of the 
spot made by a rapidly moving cathode ray stream would be of especial interest. 
With regard to the bromide paper, where precisely was it placed, and was it 
specially treated in any way? 


Mr N. L. Harris. Dr Leyshon’s interesting paper sets out a useful method of 
obtaining dynamic discharge tube characteristics. The striking and extinguishing 
potentials of the lamps as given by Figs. 4 and 5 are some 15 or 20 volts lower than 
it is usual to find amongst a normal batch of Osglim lamps. In particular, the value 
of 115 volts for the extinguishing potential of lamp Bs in Fig. 5 (A) is very unusual: 
possibly the readings of the electrostatic voltmeter are subject to a correction of 
this order. The constant voltage across the lamp during discharge under certain 
conditions is probably due to the fact that the cathode fall of potential is “normal,” 
with the cathode incompletely covered with glow. It is well known that in these 
circumstances the tube voltage is almost independent of the discharge current. In 
connection with Dr Leyshon’s recent demonstration of periodic movements of the 
cathode glow, and Dr Dye’s remarks upon the apparently irregular growth of the 
flash, it might be advisable to emphasise the important part played by surface films 
on the electrodes. The striking of the discharge is intimately associated with the 
electronic work function of the cathode surface, and this in turn with any adsorbed — 
films of monomolecular nature at this surface. A small patch of such a film of an 
electropositive element is sufficient to initiate a local discharge at a lower potential 
than that required for the clean electrode surface, and small traces of an electro- 
negative element have the reverse action. It is wellnigh impossible to maintain — 
a discharge tube free from the traces of impurity sufficient to cause such effects. 
The presence of such films, which change constantly under the influence of the 
discharge, plays an important part in bringing about many of the phenomena 
described. The existence of these films is demonstrated very simply by the increase 
of striking potential noticed after a discharge tube has stood unused for some hours. 
An adsorbed film of gaseous impurities will then have formed on the cathode, and 
the striking voltage will be raised. The passage of the discharge cleans the surface, 
in part at least, and the striking potential falls. 


AUTHOR'S REPLY. I should like to thank the various speakers for their en- 
couraging Comments on and criticisms of the paper. 

The discontinuities in individual discharges mentioned by Dr Dye were not 
observed in my experiments—no doubt because each recorded figure represented ; 
an average of some thousands of discharges. The distinctness of the figures seen — 
showed, however, that the variation of the greater number of individual discharges 
from the average could not have been very marked—at least after the initial rapid — 
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rise of current. It is possible that the resistance in the discharge circuit may have 
had a stabilising effect. 


The original object of the investigation was to obtain some idea of the nature 
and duration of the condenser discharges through the neon lamp, these repeated 


discharges having proved to be a suitable form of stimulus for excitation of nerve 


preparations. 
In Dr Wood’s interesting records the voltage changes across the condenser 
during discharge are shown by the faint lines, which illustrate the difficulty in 


_ recording discharges of such brief duration by a method involving a single 


traverse of the oscillograph spot across a photographic plate. 

The explanation for the change in striking voltage which is shown so clearly 
in Dr Wood’s records has been given by Mr Harris in his contribution to the 
discussion. 

The electrostatic voltmeter I used gives readings which correspond closely 
with those of a Weston instrument and of a Cambridge Instrument Company’s 
voltmeter. The lamp B5 for which records are given in Fig. 5(A) certainly has a 


_ lower minimum voltage than any other lamp which was examined. I have lately 


redetermined the minimum voltage for this particular lamp, and have confirmed 
previous observations. The shapes of the characteristics of this lamp are, however, 
like those of other beehive lamps. 

In reply to Mr Carpenter, the cathode ray oscillograph used was of the standard 
Western Electric type (No. 224 A). Its sensitivity was about 1 mm. movement of 
the oscillograph spot for 1 volt potential-difference between the deflecting plates. 
The bromide (gas light) paper used was not specially treated. It was held down on 
the top of the oscillograph tube for about half a minute. This method of obtaining 
photographic records with bromide paper was described by Dr Dye in a paper 
read before the Physical Society*. 


* Proc. Phys. Soc. 37, 158 (1925). 
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ABSTRACT. Gas “ degeneracy” (the effect of non-Maxwellian distribution of molecular 
velocities) is distinguished from gas “‘imperfection”’ (the effect of intermolecular forces), 
and it is shown that specific heat measurements cannot be used to verify theories of the 
former, unless account is taken of the very rapid increase of the latter at the lowest tempera- 
tures and highest pressures. Accordingly a thermodynamic method is used to calculate 
the sum of degeneracy and imperfection at 4° and 5° Abs. for helium. Theories of de- 
generacy can then be roughly tested by the extent to which they allow a residual value of 
the imperfection to be calculated, and determine a consistent extension of the course taken 
by isothermals down to 15°. The method is applied to Fermi’s expressions for degeneracy, 
and also to Berthelot’s equation of state. One of the equations alone satisfies the required 
condition; it is shown that on this basis degeneracy would comprise 15 per cent. of the 
total departure of helium from the ideal gas laws at 4° and 5°, the remainder being due to 
true imperfection or intermolecular forces. 


§x1. STATEMENT OF THE PROBLEM 
ie the ideal gas be defined, for the present purposes, by its obedience to the 
equation 
aio p-kkT Eee (1), 


where p is its pressure in dynes/cm.*, T its absolute temperature, n the number of 
molecules per c.c., and « is 1:372 X 10718, 

Any actual gas departs from this law owing to the finite size of its molecules and 
the forces of attraction and repulsion between them, and we follow normal usage 
in referring to this departure as the ‘‘imperfection” of the gas. But for many years 
in thermodynamics (Nernst’s heat theorem) and more recently in kinetic theory 
(the inclusion of translatory motion under quantum restrictions) it has seemed 
probable that even a gas which had no imperfection in the above sense would never- 
theless still show departure from law (1) at high pressure and low temperature. In 
Statistical mechanics this accompanies a departure from Maxwellian distribution 


of molecular velocities, and has been called the “degradation” or “degeneracy” 
of a gas. 


=< 
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To include both the imperfection and the possible degeneracy of a gas, its true 
equation of state would have to take the form 


ome ec Tee Le FLD) enh peed i (2), 


| where F(T, p) is a function of temperature and pressure representing the effect of 
molecular dimensions and forces, and f (7, p) represents the effect of departure 
from Maxwellian statistics. 

For brevity we shall refer to these two functions as A and B respectively. 

The very large number of empirical equations of state which have from time 
to time been proposed, beginning with that of Van der Waals, deal exclusively with 
A, and there is no doubt that B is only likely to become important under the most 
extreme conditions of pressure and temperature realisable in practice. Attempts 
to verify the predictions of degeneracy under those conditions made by thermo- 
dynamic and kinetic theory have mainly been based on the following argument. 
Even the most complex equation of state must satisfy the general relations 


E = p/n, 
C, = (dE/aT),, 


where E is the energy of the m molecules and C,, the specific heat at constant volume. 
Hence degeneracy must affect C,,, and anomalies discovered in calorimetry may give 
information about B if the substance can be made to remain in the gaseous phase 
at low enough temperature. Hydrogen and helium obviously present the most 
favourable cases, because of their comparatively slightimperfection and consequently 
low boiling points. At temperatures so low that the contribution of rotational and 
vibrational degrees of freedom to the specific heat have vanished, both these gases 
should, in the ideal state, obey the law 


ae, 63 
C= gk. 


The measurements of Eucken, and more recently of Meissner, reveal a distinct 
fall below this value for both gases at very low temperatures, and hence these 
experiments have often been quoted as experimental verification of degeneracy. 
For hydrogen this interpretation is complicated by the recently found difference 
" between the rotational specific heats of the ortho- and para-molecular varieties, at 
temperatures just above+hat at which those quantities vanish (3). However, it is true 
that any survival or anomalous excitation of rotational degrees of freedom could only 
increase C,,, and so cannot be responsible for the facts; further, all theories of de- 
generacy do predict a tendency in the observed direction by their requirement that 


z Geo as) > 0. 


Nevertheless, for any such experiments to afford a convincing test of theoretical 
_ forms of B,.or even a decision as to whether B is large or small under the most 
extreme practical conditions, the course taken by the function A would have also 
to be known in this region. 
Such a requirement is not easy to fulfil; we shall indicate below some incon- 


sistences produced by extrapolating, to temperatures at which degeneracy is likely 
; 12-2 
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to be important, empirical equations of state known to represent A satisfactorily . 
under less extreme conditions. We are left with the suspicion that effects attributed 
to B may be due to an unknown increase in A, and any independent clue as to the 
relative importance of imperfection and degeneracy becomes worthy of enquiry, 
Hence in the present paper the problem of the possible place of degeneracy in 
experiment is investigated from a different standpoint, as follows. 


wr 


§2. METHOD OF THE PRESENT PAPER 


A method is devised whereby the vapour pressure and latent heat data of the 
Leiden experimenters are made to yield values of (A + B) through the latent heat 
equation of thermodynamics. Now direct measurement of helium isothermals has 
been carried down to 15° Abs. by the Berlin Reichsanstalt. From their data we 
plot a curve of the departure of the gas from the ideal law between go” and 15°. This 
represents the course of (A + B), but it is very probable that, even at 15° and the 
given pressure, B is small compared with A, and that the curve substantially re- 
presents the course of A. We show that the further points secured by our method 
at 4° and 5° lie on a smooth continuation of this graph of (4 + B). But whereas B 
was probably negligible at 15°, it is not necessarily so at the lower temperatures. 
Accordingly it becomes of interest to substract from our values of (4 + B) any given — 
theoretical value of B; it is then possible to see at once whether the residual value 
of A is or is not still consistent with the course entered upon by (4 + B) at the ; 
higher temperatures for which B was negligible. We can in this way make a rough — 
test of any theory of degeneracy, by requiring that the value of 4 which it extracts 
from our (A + B) must not lie in such a part of the diagram as to bring obvious 
discontinuity into the course of that function. 

This indirect test is capable of eliminating impossible theories rather than of 
distinguishing between those which satisfy the conditions, since one extrapolation 
can be shown to be wrong when another cannot be proved to be right. It does not, 
however, involve the assumption that an equation of state has equal validity on both 
sides of the critical temperature, and thus it may have some advantage over the only 
alternative method by. which theories of degeneracy can be compared with experi- 
ment. 

The method is applied to Fermi’s two expressions for degeneracy. A similar 
test is made of the validity of Berthelot’s equation of state in this region. 


§3. TOTAL DEPARTURE FROM GAS LAWS AT 4:29° ABS. AND 
520° ABS.,. DERIVED BY INDIRECT METHOD 


If equation (2) be written in the form 
p=nT(1+A+ B), 


where A and B are the imperfection and degeneracy functions of J’ and Pp, the total 


(A + B) can be calculated if pressure and density and atomic mass are known, 
since we have 


n= p/m, 
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Consider the latent heat equation 


L = (1/py— 1/p1) T (dpa), 
where L is the sum of external and internal latent heats, py, and p, are densities in 
‘the vapour and liquid phases respectively, and p is the saturated vapour pressure 
at temperature 7. Kammerlingh Onnes“ determined the density of liquid 
helium at 4:29° Abs., and the latent heat of evaporation at the same temperature. 
From the vapour pressure experiments of Onnes and Weber“ we construct a graph 
giving us the pressure and pressure gradient at that temperature. From all these 
sources we obtain the following: 
d==sA720-5 
L = 22 cal./gm. mol. 
2-20.50 10° erg; 5m, 
py = 07122 gm./c.c. 
dp|dT = 9:37 x 10° dyne/cm.? per degree, 
/ whence we calculate from the latent heat equation 
py = 070152 gm/jcc., 


a result which is not inconsistent with the statement from the other cryogenic 
laboratory” that at a slightly higher temperature the density of the vapour is about 
one-tenth of that of the liquid helium. 

For the mass of the helium atom we have 


m = 6°56 x 10-” gm. 
For Boltzmann’s constant 
ie 12372 10-**. 


Hence at 4:29° 
Hie] = 1301 X 10%. 


“From the curve plotted with the data of Onnes and Weber, at this temperature, 
p = 810 mm. Hg. 
= 1082 x 108 dynes/cm.?. 
2 A+ B=— 0204. 

The total departure from the ideal gas law can again be calculated at 5-20° Abs. 
This is the critical temperature at which L vanishes. But, using the critical co- 
efficient for helium, we have 

KT ,[PeVe = 3°13) 
Pe = 1718 mm. of mercury 
= 2:29 x 10° dynes/cm.”. 


- Ve= 99r4-x 10> Ce, 
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. . ce eee >) 
This is the critical volume for a single atom, since we have used the above atomic 
value for x. We accordingly combine it with the above value of the atomic mass 
and obtain 
p = 0-066 gm./c.c. 
Hence at 5:20° 
nicl =.7-20 % Is 


’, A+ B=— 0682. 


§4. TOTAL DEPARTURE FROM GAS LAWS AT HIGHER 
TEMPERATURES 


These indirect determinations of the total departure from the ideal state at 
4:29° and 5-20° must be compared with direct determinations at the lowest tempera- 
ture available. In 1926 Holborn and Otto® at the Reichsanstalt investigated 
isothermals for helium down to 15° Abs.; we calculate from their data the course 
taken by (A + B) at the two pressures we have been considering, obtaining two 
curves over the range go° to 15°. Beyond this range no results are obtainable except 
by indirect methods. The positions of points representing our values at 4° and 5° 
are therefore of some interest; we enquire later whether their fit on to the curves is 
improved or spoilt by assuming any given theory of degeneracy. 

Below 60° Abs. the change in the helium isothermals is so rapid that considerable 
care has to be taken in choosing any way, graphic or algebraic, of enabling their 
course to be shown without ambiguity. The usually adopted method of plotting 
a first-order correction term to Boyle’s law is not adequate to this case. We take 
(Table 1) the data of Holborn and Otto from their empirical expression, 


po= A+ Bp+ Cp?+ Dp*+ Ep*. 


Table 1. Coefficients in Holborn and Otto’s equation. 


Absolute 
temperature A Bx 10° Cx 16% D x 108 E x 10° 
gor 0°32992 062286 0°735 _ — 
65°1 0'23847 055080 2°384 —O°oI4I _ 
20°3 0'07460 —0°16424 18°529 —O°LIO5 — 
151 005558 —0°79740 54°370 —0°7513 3°796 


After unsuccessful attempts to make an unambiguous extrapolation of each 
coefficient separately, logarithmically and otherwise, it was found more definitive 
to calculate the integrated value of pv and subtract it from the normal value of the 
ideal gas. We plot the differences (pv — RT) as a function of temperature for the 
two conditions with which we have been concerned, namely, 810 mm. and 1718 mm. 
The result, transferred into the units of Holborn and Otto, is shown in Fig. 1, the 
unit of pressure being one metre of mercury, while the unit of volume is the volume 
of one gram of gas at that pressure and o° C. 

From the numerical equation of state which we have derived at the lower tem- 
peratures, we determine, transferring our values to the units of Holborn and Otto, 
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comparable expressions for (pv — RT) at 4:29° and 5-20° and the pressures 810 mm, 
and 1718 mm. 

It is seen that these points (+ in Fig. 1) fit a not unreasonable extension of the 
curves beyond the lowest temperatures at which direct data are available. Since the 
curvature of both graphs is increasing rapidly at 15°, more than one set of values 
at 4° and 5° might be regarded as agreeing with their extrapolation. Here again 
logarithmic plotting was not found to remove possible ambiguity, and the rigorous 
identification of any extrapolation would require an empirical equation involving 
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Fig. 1. Gas imperfection, degeneracy and temperature. 


higher order terms to which no clue is available at less extreme temperatures. Hence 
arguments from these curves are only valid in particular cases where decision has 
to be made between very divergent alternatives. 
So long as this restriction is regarded, we can say that, if a single curve from 4° 
~ to go° represents (A + B) and over the upper half of it B is negligible, a broad limit 
can at once be set to the allowable degeneracy. Either B at the lowest temperature 
must be small enough not to alter our extension of (A + B) appreciably, or else if 
B is large A must be supposed to have increased, with opposite sign, enough to 
‘prevent their algebraic sum from departing noticeably from the curve. In the 
particular case discussed below, the numerical values are such as to permit decision 


between the alternatives. 
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§s5. THEORETICAL VALUES OF DEGENERACY 


We turn to theoretical expressions for the function B, of which the basis has ‘| 
been twofold. ) 

(a) Degeneracy arises from the application of Nernst’s heat theorem to non- 
condensed systems, and is there inseparable from Nernst’s view of the existence of 
zero-point energy. Treated thermodynamically, the most complete exploration of 
its consequences has been that of Bennewitz, who, however, did not succeed in 
deriving in quantitative form a function which we could treat as B for any finite 
temperature and pressure. 

(b) Independently of Nernst’s heat theorem, Einstein, Bose, Dirac, Fermi and 
others have investigated degeneracy as a necessary consequence of statistical 
theorems on the non-Maxwellian distribution of molecular velocities. Confidence 
in these theorems is mainly due to their success, in the hands of Sommerfeld and 
others, in explaining otherwise intractable phenomena in other branches of physics. 
We select as a case for treatment by the present method two approximations derived 
by Fermi® for ‘“‘strong” and “weak”? degeneracy, and expected to describe 
phenomena under the most extreme and less extreme conditions respectively. The 
border-line between the regions to which they could apply, and indeed their 
validity anywhere, has not been otherwise tested. 

Strong degeneracy: Fermi predicted the following relation between p and , x, T, 


p = 0°077 h?n?!3/m + 7-6 m nh8x?T?/h?, 


in which I have reduced to numerical constants certain elements in his expression 
which are independent of physical data. The only symbol which has not previously 
been defined is h, Planck’s constant, equal to 6-55 x 107*?, 

Weak degeneracy : 


p = nT {1 + h8n/16 (7mxT)>?}. 
We insert in these equations numerical values appropriate to the cases we have 
considered, namely, 
T=4:29° p=0o-0152, 
T' = 5:20° p = 0°066, 
and compute the degeneracy function for an otherwise perfect gas in the form 


p= nT (1 + B). 
If weak degeneracy occurs, 
B=+ 0031 at 4:29° 
=+ 010 at 520°. 
If strong degeneracy occurs, 
B=-+ 3:48 at 4:29° 


=+ 1°30 at 5°20°. 


| 
| 


eri 
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§6. INFERRED VALUES OF “IMPERFECTION” 


f Subtracting these theoretical values of B from the values of (A + B) deter- 
mined by our indirect method, we find that: 
weak degeneracy being assumed, 
A =— 0235 at 4:29° 
=— O75 at 5°20°5 
strong degeneracy being assumed, 
A=— 3°68 at 4:29° 
“=== 1798 “at 5°20". 
Applying the method of test outlined previously, we see that strong degeneracy, 
if existent in this case, would cause the imperfection to be 18 times and 3 times the 
total departure from the ideal law at 4° and 5° respectively. Now the rates of in- 
crease of the curves drawn from the data of Holborn and Otto, regarded as sub- 
stantially graphs of A, are such that our values might still appear as extrapolations 
thereof, even after adjustments by, say, 50 per cent. or even 100 per cent. But 
| adjustment by factors of 3 and 18 would introduce into them a character having 
no connection at all with their courses at 15°. Further, it would involve the arbi- 
trary assumption that our values of (A + B) are due to a remarkable coincidence, 
the degree of neutralisation of + 3-48 by — 3-68 having the same proportional 
effect on the one curve as the degree of neutralisation of + 1-30 by — 1-98 had on 
the other. This would be a very forced explanation, compared with the simple one 
that B is small compared with A. Again, if strong degeneracy held at 4° and 5°, the 
above figures show that it would make the value of A at 810 mm. much greater than 
at 1718 mm. Reference to the curves, which are on suitably different scales because 
of the greater order of magnitude of A at 1718 mm. than at 810 mm., shows the 
extreme unlikelihood of this being so. The higher-pressure curve increases much 
faster than the lower, and we have no right to demand that they should converge 
again, since the mechanism of imperfection must be the same in both. 

On the other hand, the expression for weak degeneracy introduces a correction 
to (A + B) which is more convincing, in that it enables the Holborn and Otto curves 
to be extended in the same relationship to one another which they had already 
shown. The rate of increase of B relative to that of A is then similar on each curve. 
In fact at each pressure the degeneracy amounts to about 15 per cent. of the total 
departure from the behaviour of an ideal gas. The degeneracy correction to the 
pressure is about 3 per cent. and ro per cent. at 4:29° and 5-20", the pressures being 
rather over one and two atmospheres, while for the imperfection correction the 
corresponding figures are of the order of 20 per cent. and 8o per cent. and of 
opposite sign. In Fig. 1 the points x represent A, assuming weak degeneracy. 

Owing to this opposite sign, if conditions are reached at which the strong 
approximation becomes valid, we may expect a cancellation of imperfection by 
degeneracy, and a second “Boyle” point, at which again the behaviour approxi- 
mates to that of the ideal gas through the condition 

= —B. 
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§7. EXTRAPOLATION FROM BERTHELOT’S EQUATION 

It has been stated above that in this type of problem it is unsafe to assume that 
any equation of state is valid when approaching the absolute zero, however well it 
may fit the facts in another range. It is of interest to see what results would be 
obtained if this precaution were neglected and A calculated by extrapolating some 
equation known to be a good approximation at temperatures high compared with 
the critical temperature. Berthelot’s equation is the most successful of these; it 
may be stated in the form, 


eri 4 2 Bef, SM 
pomRT + Sea mie 


We calculate the values of A corresponding to those already deduced, inserting the 
appropriate values of 7, and p, and taking T= 4-29 and 5-20 and p = 810 mm. and 
1718 mm., and obtain 


A=— 0-314 at 429° 
=— 0351 at 5-20°. 
As might be expected for an equation designed for less extreme circumstances, this 


predicts a less rapid increase of departure from the ideal law than that found by . 
the indirect method, whether degeneracy is assumed or neglected. 


§8. CONCLUSION 


Since an unknown increase of gas imperfection at very low temperatures might 

be responsible for calorimetric effects attributed to gas degeneracy, the need arises 
for indirect methods of checking theories of degeneracy and of estimating the 
numerical importance which any of them would bear, relative to imperfection, in 
the most general equation of state. The method outlined in the preceding sections 
allows the total departure from ideal gas laws to be computed from vapour pressure 
and latent heat data for helium at 4° and 5° Abs. 
; In comparison, Berthelot’s equation of state fails to give a sufficiently rapid 
increase of imperfection. Fermi’s “‘strong’’ approximation to the degeneracy 
function also fails to justify the assumption that it can occur at these temperatures 
and pressures, since it would require too rapid an increase in the imperfection 
function, and would require the more rapid increase to take place at the lower 
pressure. On the other hand, Fermi’s “ weak” approximation provides a degeneracy 
consistent with curves constructed from direct data above 15° and indirect data 
below 15°. 

The proof that Fermi’s weak degeneracy is not inconsistent with the known facts 
does not, of course, establish its existence. If, however, it is regarded as worth 
consideration on this basis, we show that the share taken by degeneracy and im- 
perfection would be as follows: 


at 4:29° and 8 ‘ 
4°29 10mm p = nxT {1 — 0:235 + 0-031}, 


p = mT {1 — 0-78 + 0-10}, 


where the two last terms represent respectively the effect of intermolecular forces 
and of non-Maxwellian distribution of molecular velocities. 


at 5:20° and 1718 mm. 
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DISCUSSION 


Prof. W. Witson. I should like to ask the author how these low temperatures 
have been estimated. I had the impression that some equation of state such as 
that of van der Waals was assumed at the outset in order to determine them. It is 
interesting to note that Sommerfeld has recently applied the new statistical methods 
_to the electron theory of the conduction of heat and electricity in metals. The 

classical theory gave far too big a specific heat to the metal. This is now corrected, 
I believe, by the degeneracy of the electron gas in the metal. 


Mr J. H. Awsery. Professor Wilson has mentioned the use that Sommerfeld 
made of the theory of degeneracy. If I remember correctly, Sommerfeld supported 
his argument that electron gas should be degenerate at all temperatures known 
_terrestrially, by comparing its degeneration quantitatively with that of hydrogen 

and helium as found experimentally (using the well-known molecular weight and 
temperature relation). His only exact evidence of degeneracy at that time was, of 
course, from specific heat data which, as mentioned in Mr Johnson’s paper, give 
results that can only be attributed to degeneracy and not to imperfection. The 
fact that such uses have been made of the theory only increases the importance 
of an estimate such as Mr Johnson has made, and I should like to congratulate 
him on getting so much out of such scanty data. 


Prof. J. E. Lennarp-Jones. The new statistical theories developed by Bose 
and Einstein on the one hand, and Dirac and Fermi on the other, diverge most 
from classical statistics at low temperatures or high pressures, and calculation 
shows that it is only in the case of helium near its critical point that its behaviour 
is likely to be perceptibly different according as it obeys the classical, or the 
Bose-Einstein, or the Fermi-Dirac statistics. It is, therefore, a matter of great 
theoretical interest to determine the deviations of helium from the perfect-gas law 
near its critical point and to analyse the effect of the statistics and of the inter- 
atomic fields. Mr Johnson has made an interesting attempt on this problem, but 
- there are one or two points in his treatment which seem to me to require investiga- 
‘tion in any further work. 

In the first place, the assumption is made that the effects of interatomic fields 
and of the statistics are simply additive, equation (2). As, however, the calculation 
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of the effect of interatomic fields must itself depend on the statistics assumed, this 
simple separation of the two effects requires justification. In the second place, the 
author considers only the Fermi-Dirac statistics, whereas theory indicates that 
helium atoms should obey the Bose-Einstein statistics. It would add to the value 


of his work if the author could consider the effect of the latter statistics on helium — 


near the critical point. 


AuTHOR’s REPLY. Gas thermometry, at temperatures where imperfection and 
degeneracy are in question, is still accurate to a first approximation if the thermo- 
metric gas is at low pressure: this condition can be maintained while the gas whose 
temperature is measured may be at high pressure. Nevertheless, the question of 
Prof. Wilson and the first point raised by Prof. Lennard-Jones force us to recognise 
two or more stages of approximation in these problems. First assume that the gas 
laws hold and from the temperature readings deduce the amount of the consequences 
in imperfection and degeneracy. An equation of state can thence be constructed 
which may be used to correct the original temperature reading. With this corrected 
temperature a slight adjustment of the imperfection and degeneracy may then be 
made, according to the method, familiar in other branches of physics, of assuming 
that a certain error is zero in order to obtain a first approximation to its finite value, 
and using this to remove the assumption in a second approximation. Similarly, 
to the first approximation the imperfection and the degeneracy must be regarded 
as additive in order to determine their order of magnitude. Now that Leiden and 
Berlin are both obtaining data at the same low temperatures, the second approxima- 
tion suggested by Prof. Wilson’s remarks may soon be practicable, though I think 
the data are not yet accurate enough for the other second approximation which 
Prof. Lennard-Jones points out as ultimately necessary. 
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ABS TRACT. A magnetostriction constant K for the Joule effect is defined by the equa- 
tion p = KH where p is the alternating mechanical stress produced by a small alternating 
magnetic field H superposed on a steady magnetic field H,. It is shown that the alter- 
nating intensity of magnetisation J produced by an alternating strain 6/// is given by the 
equation J = K8//1; K being the same constant in both equations. These equations are used 
to calculate the motional impedance of a laminated ring toroidally wound and vibrating 
in its fundamental radial mode. K has been measured for various values of H, with a 
ring of the alloy Fe 60 per cent., Ni 40 per cent. It is shown that the magnetostriction 
effect for steady magnetising fields can be inferred from the measurements with alternating 

_ magnetic fields, if the reversible permeability and the differential permeability of the alloy 
are known. 


§1. INTRODUCTION 


ue need for accurate and constant substandards of frequency has led to the 
use of oscillating valve circuits controlled by very resonant mechanical 
vibrators. The valve-maintained tuning fork* is available for sonic frequencies 
and the piezo-electric crystal t for radio frequencies. These vibrators are not con- 
venient for frequencies between 5000 and 25,000 ~ and Pierce has developed the 
magnetostriction oscillator to. cover this range. The theory of the valve-maintained 
tuning fork has been given by Butterworth§ and others and that of the piezo- 
electric crystal by Dye|| on the assumption that the vibrators can be considered 
as mechanical systems having one degree of freedom. This assumption is equally 
valid for the magnetostriction oscillator, and the theory appropriate to it must be 
‘closely analogous to that for the tuning fork and the piezo-electric crystal. It is 
- essential in the three cases to know the nature and magnitude of the coupling 
“between the mechanical vibrator and the electric circuit. This coupling is defined 
by the piezo-electric equations and the piezo-electric constant for the crystal and 
there must be similar equations and a similar constant for the magnetostriction 
oscillator. It is the object of this paper to define this magnetostriction constant and 
to show how it can be measured. 
* W.H. Eccles, Proc. Phys. Soc. 31, 269 (1919). 
+ W. G. Cady, “ Pjezo-electric standards of high frequency,” Journ. Opt. Soc. of America, 10, 


475 (1925). 
| + G. W. Pierce, “ Magnetostriction Oscillators,” Amier. Acad. Proc. 63, 1-47 (1928). 


§ S. Butterworth, “The maintenance of a Vibrating System by means of a Triode Valve,” Proc. 


Phys. Soc. 32, 345 (1920). Enea 
|| D. W. Dye, “The Piezo-electric Resonator and its Equivalent Electrical Circuit,” Proc. Phys. 


‘Soc. 38, 399 (1926). 
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The magnetostriction oscillator used by Pierce* and Vincentt is generally a thin 
rod and the coupling to the electrical circuit is obtained by placing itina solenoi 
carrying a small alternating current superposed on a steady magnetising current. 
The stresses arising from the magnetostriction effect set the rod into resonan 
vibration in its fundamental longitudinal mode, with a node at the centre and anti 
nodes at the ends. The mechanical stress is greatest at the centre and falls to zero 
at the ends, with the result that the intensity of magnetisation must vary along the 
rod. There is a further complication due to the demagnetising effect of the magneti 
poles appearing at the ends of the rod. The theory is much simplified by considera: 
tion of a magnetostriction oscillator in the form of a closed circular ring which is 
magnetised by a uniform toroidal winding. A ring has a mode of vibration in which 
all parts of it move radially with the same amplitude and phase and, if this mode is 
excited by the passage of an alternating current of suitable frequency through the 
toroidal winding, the magnetisation and the mechanical stress cannot vary along 
the circumference of the ring. There is, however, a variation in these quantities over 
the cross-section of the ring, since eddy currents are induced in it. If the ring can 
be laminated this variation can be much reduced and the stress and magnetisation 
become practically uniform throughout the material. Before investigating the 
motion of the ring and its reactions on the electrical circuit, it is necessary to define 
the magnetostriction constant in a form suitable for the study of alternating effects. 


§2. THE MAGNETOSTRICTION EQUATIONS 


The Joule effect is usually defined as the steady mechanical strain produced by 
the application of a steady magnetic field. This definition is inadequate for a study 
of the alternating strains produced by a small alternating magnetic field superposed 
on a steady magnetic field, for these strains depend on the frequency of alternation 
and the possible resonance of the vibrating body as well as on its magnetostriction 
properties. It is preferable to define the Joule effect as a relation between the 
alternating magnetic field and the alternating mechanical stress produced by it. 
If the laminated ring described above is imagined to be rigidly constrained so that 
no motion is possible, the alternating magnetic field from the toroidal winding will 
produce the same alternating mechanical stress at all frequencies and we may write 

P= KO 9° «@ ke 4) eee (1) 
where p is the stress in dynes/cm.? H (or H, sin wf) is the alternating part of the 
magnetic field in gauss, w (or 27 x frequency) is the pulsatance and K is a magneto- 
striction constant. There is a steady stress, produced at the same time by the steady 
magnetising field, with which we are not concerned. On removal of the constraints, 
the ring vibrates and we may regard the stress p as the force maintaining the vibra- 
tion. ‘The equation of motion can now be derived in the usual form for a vibrating 
system with one degree of freedom. Let &, and = represent the instantaneous 


* Loe. cit. 


i . H. Vincent, “Experiments on Magnetostrictive oscillators at Radio Frequencies,” Proc. 
Vs. Soc, 41, 476 (1929); Nature, 120, 952 (1927); Electrician, 101, 729 (1928); 104, 11 (1929). 
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outward radial displacement, velocity and acceleration at time ¢. Let 7 be the 
radius, o the mean density and S the cross-section of the ring and let C be the 
velocity of sound in it and E the modulus of elasticity defined by the equation 
C? = E/c. The circumferential stress p is equivalent to a total force, tending to 
expand the ring, of 27Sp or 27SKH. 'The equation of motion is therefore 


anSKH = anrSot + RE + 2nSEE/r 
where the terms on the right-hand side represent in the order given the forces 


necessary to overcome the inertia, the damping and the stiffness. The symbolic 
solution of this equation for continuous and steady motion is 


g = anSKH 
k + jw (amrSo — Chr) meee 8 (3) 


where j = V — 1. The displacement is a maximum, that is the ring resonates, at 
a pulsatance w, given by 


2nrSoa — 27SE/ra,? = 0, 


from which wy = VE/o x 1/r = cfr 
and ariptc ee hae (4). 


Resonance occurs at the frequency for which the wave-length of sound in the 
material is equal to the perimeter of the ring. 
The velocity at resonance is 


£=anSKH/k=20SK (4nniy[R wees (5) 


‘where 7 is the number of turns per cm. in the magnetising winding and 7 is the 
alternating part of the magnetising current. 
“The motion of the ring having been found, it remains to calculate the electro- 
motive force induced in the toroidal winding as a consequence of the alternating 
Villari effect. An alternating mechanical strain ¢/r is accompanied by a change J in 


the intensity of magnetisation, where 
(ees or oe eee ees (6), 
K’ being the magnetostriction constant for the Villari effect. The electromotive 


force v induced in the winding is equal to the rate of diminution of the magnetic 
linkage through the winding, that is 


vy = — 2anrnS x 4nk from (6) 
— — anrnS x 4nK’€/r from (6) 
= — (27)? (47)? WS°KK 7k from (5) 

saree (4). 


This expression for v shows that it is in opposite phase to the current i, and the 
electrical power required to maintain the vibration is therefore 
(CDC gS Sal ee na (8), 
which must be equal to the power consumed by the dissipative forces in the ring, 
that is, to RE. Equation (5) shows that ké? is equal to 
CLS Cd soa) Say (9). 
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A comparison of (8) and (9) shows that K = K’ and therefore one magneto- 
striction constant is sufficient for the alternating Joule effect and the conve 
phenomenon, the alternating Villari effect. The magnetostriction equations beco 

p= KH for the Joule effect ...... (10) 
and T= Ké/r for the Villari effect ...... (x1). & 


This reciprocal relation between the two effects can be derived from a gene 
theorem given by J. J. Thomson*; it has been proved here from first principles i 
order to present a physical interpretation of the phenomena associated with t 
vibrating ring. 

An additional equation is required for the interpretation of the measurements: 
described below. Equation (3) shows that the velocity of vibration is half that at 
resonance at two frequencies given by ; 


4 V3 k= w(20rSo—20SEfrw®) wns (12). 
If the ring is very resonant, we may write w = w, + dw where wo/27 is the resonance 
frequency given by (4). On substituting this value above and neglecting (w/a )*, 
and higher powers, we find that 
Sw=4tV3k/4amrSo eee (13). 
Apart from the use of the reciprocal relation, this investigation differs from that 


of Piercet in its definition of the modulus of elasticity. He expresses the magneto- 
striction force{ in the form 


p= aS (Bi+ B, 
where Bz is the induction due to the periodic current 7, B’ is the induction due to 
the reverse magnetostriction effect, that is the Villari effect, and a is a magneto- 
striction constant for the Joule effect. The term aSB’ is a force opposing the motion 
which is included in Young’s modulus for the material and should not appear in 
his equation of motion§ as a term to be subtracted from Young’s modulus. In fact 
a new modulus of elasticity which he introduces later) is actually Young’s modulus ; 
he ae is of little practical importance and does not affect his subsequent 
analysis. 


§3. THE MAGNETOSTRICTION CONSTANT 


A quantity of the nickel-iron alloy containing 40 per cent. nickel was available 
in thin sheets 0-002 in. thick and, Pierce having reported that it vibrated well, an 
attempt to measure its magnetostriction constant was decided upon. A numba of 
circular rings of the alloy of 11-05 cm. and 8-55 cm. in external and internal diameter 
respectively, and a number of similar paper rings 0-0005” thick were prepared. ‘The 
metal rings were cleaned, covered on both sides with a thin layer of paraffin wax 
and piled alternately with the paper rings. The composite pile was clamped between 
two metal plates which were sufficiently hot to melt the wax and, after slow cooling, 


Ait sales LB ok : ‘cath ~ . : : 
rele} J omson, Applications of Dynamics to Physics and Chemistry, ch. 5. (Macmillan and Coy 


t Loc. cit. p. 35 and p. 40 + Loch ° 
J * Cc. Ctr. P. % X 
dec. Bit. py 56) ditectae Pp. 35 (2) § Loc. cit. p. 36 (5). 
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the finished laminated ring was taken from the press and superfluous wax cut away. 
The ring was placed on three rubber supports inside a circular wooden channel 
over which two toroidal coils of 227 and 215 turns were then wound (Fig. r). The 
227-turn. winding carried a direct current supplied by accumulators through a large 
‘inductance and a variable resistance, and the other winding completed one arm of 
the alternating current bridge with which the measurements of motional impedance 
were made. The ratio arms of the bridge were non-inductive resistances of 100 ohms, 
and the remaining arm contained a calibrated inductance in series with a calibrated 
resistance, both continuously variable. Alternating current was supplied to the 
_ bridge from a coil inductively coupled to a 20-watt valve oscillator, the frequency 
of which was continuously variable. An earthed wire grid interposed between this 
coil and the valve oscillator screened the bridge from the electric field of the 


SS 


— | 
sathsitttdla, 


Polarizing Winding 
Alternating Current Winding 


Fig. 1. Circuit for measurement of motional impedance. 


oscillator. The detector was a thermionic voltmeter reading from o to 1-5 volts and 
connected to the bridge through a tuned transformer. ‘These simple arrangements 
were quite adequate for the measurements of impedances up to 100 ohms at 18,000~. 
The stray capacities which existed between the elements of the bridge were certainly 
less than 50 pf. and, since this corresponds to a reactance of 178,000 ohms at the 
frequency used, the resulting error in the balance of the bridge was negligible. 
Some care was required to bring the laminated ring to a reproducible magnetic 
state. The ring retained some of its magnetism when the magnetising current was 
broken and the intensity of magnetisation therefore depended upon its previous 
_magnetic history. It was essential to demagnetise it before applying the required 
magnetising field. Residual magnetism was easily detected on the alternating 
bridge by the presence of motional impedance at frequencies in the neighbourhood 
of resonance and it was possible though not easy to remove it by momentary appli- 
cations of a magnetising field. The sound emitted by the ring provided another 
sensitive test for the presence of residual magnetism, for this sound could be easily 
heard when ‘any magnetism remained. ‘The completely demagnetised ring vibrated 
‘at twice the frequency of the alternating current, and the sound emitted was, of 
course, far above the audible limit. 
Measurements of the effective impedance of the alternating-current winding 
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were made over a range of frequencies in the neighbourhood of resonance.. Near 
the resonance frequency, a small continuous drift in the frequency of the valve | 
oscillator, as the filament battery discharged, was accompanied by a comparatively 

rapid change in the impedance and, as two adjustments were necessary to obtain 

a balance, the latter was somewhat elusive. It was simpler to set the variable in- 
ductance to a suitable value and to obtain a balance by adjusting resistance and 

frequency. The measured impedances are plotted in Fig. 2 as points on a rectangular 

co-ordinate diagram with reactances as ordinates and resistances as abscissae. They 
lie closely on a circle and the resonance point lies near the end C of the horizontal 

diameter BC. If BA and BD are drawn so that 

Z ABC = Z DBC= 60°, AB= DB=3B8C, 


we know, from the elementary theory of the circle diagram, that the amplitude of 

vibration at the frequencies corresponding to A and D is half that at the resonance 

point C. As will appear later, it is necessary to know the frequency increments 

corresponding to the arcs AC, DC in order to calculate the mechanical frictional” 
forces opposing the motion of the ring. The frequencies marked on the diagram, 

which were measured with a resonance wave-meter, were not sufficiently accurate 

for this purpose. Small changes of frequency were more accurately found by 

measurement of the beat note between the valve oscillator supplying the bridge and 

another valve oscillator constant in frequency. The bridge was set successively to 
the impedances corresponding to the points A, C, and D, and after it had been 
balanced by an adjustment of frequency, the increments of frequency corresponding 

to the arcs AC, DC were found by measurement of the change in the beat note. 

It was found that the arcs AC, DC corresponded to a frequency change of 

2I-5+1~. 

The circle diagram of Fig. 2 provides a sensitive test for the presence of eddy 
currents in the laminated ring. Pierce* and Blackt have shown that the effect of 
eddy currents is to depress the diameter BC, drawn through the resonance point, | 
below the horizontal. The angle of depression was about 80° for a solid rod and_ 
about 40° for a nickel tube split longitudinally to reduce eddy currents. The motional : 
impedance of the laminated ring at resonance is seen to be practically a pure re-_ 
sistance and it may be concluded that the eddy currents were too small to be 
appreciable. : 

‘The resonant ring thus complies with the assumptions in the theory given above, - 
which can therefore be used to interpret the experimental results. The motional 
impedance at resonance corresponds to the diameter BC of the circle diagram and 
we may write from (7) 

(277)? (47)? S°*KZ/k= 15-0 xX 1ro8eMU. eee (14). 
S, the cross-section of the ring, is 0-412 cm.*, and , the number of turns per cm. in 
HOE ices mie an k is found to be 13,700 gm./sec. by sub-— 
thence from (14) it a AG» Se a a ea = a 6-28 x sa 6 Beg (13); and 
cian PP tt A= 6:25 x 104. There is a correction for the fact 
aminated ring contains 85 per cent. alloy and 15 per cent. paper and — 


* Loc. cit. p. 37. + Loc. cit. pp. §5, 611 
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araffin wax by volume. The value of K is an average over the cross-section and, 
ince the paper contributes nothing to the magnetostriction effect, the value of K 
ust be increased in the ratio of 100:85. The final result is therefore 

IG ey Cis Me 

for the nickel iron alloy containing 40 per cent. of nickel and permanently mag- 
netised by a steady field of 14:8 gauss. Measurements with other magnetising 
fields are given in the table. 
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Resistance in ohms. ; 
Diameter of circle 15-0 ohms. The amplitude is half that at resonance at 21:5 ~ above or 
below resonance. Small figures in circle are frequencies in cycles per second. 
Fig. 2. Circle diagram for a vibrating ring. 
. Table. 
|} Magnetising field 3°63 727 LOL I4°5 17'8 21°8 250 31°9 
K for alternating fields | Very 
(x 104) 4 small | 177 6°7 74 76 76 74 5°3 
for steady fields Ve 
Ak 08 ard pee Le 53 42 30 2325 19'5 I2'2 
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§ 4. MAGNETOSTRICTION IN STEADY FIELDS Fj 
be shown that these measurements are consistent with the changes in 
length of the same alloy in steady magnetising fields given by Masumoto and Nara*. 
The alternating magnetic field superposed on a steady field produces changes of 
magnetic induction represented by an elongated hysteresis loop near the B/H curve 
and inclined to it, the average slope of which is known as the “ reversible”? permea- 
bilityt. The change of induction produced by a steady increment 5H of a steady 
magnetic field H takes place along the B/H curve and may be expressed in terms 
of the “differential” permeability which is the slope of the B/H curve at the point 
corresponding to H. “ Differential” permeability is always larger than “‘ reversible” 
permeability, and therefore a steady increment of magnetic field 5H produces a 
larger change of induction than an alternating field 6H sin wt. Now magnetostric- 


Tt can 


Masumoto and Nara 


Author 


oO 


E8l/l x 10-8 


ie) 10 20 30 40 
H in Gauss 
Fig. 3. Magnetostriction in steady magnetising fields, 


tion is a phenomenon resulting from the change of orientation of molecules by 
' 


magnetisation and depends upon intensity of magnetisation. It follows that the 
change of length produced by a steady increment of magnetic field 3H must be 
greater than that produced by an alternating field 5H in the ratio of the “differ-_ 
ential” susceptibility to the “reversible” susceptibility. This ratio was found by 
measurement of the slope of the B/H curve for the differential susceptibility and by 
measurement of the effective inductance of the ring on the alternating current 
bridge for the reversible susceptibility. The magnetostriction constants for steady 
increments of magnetic field are obtained by multiplication of the constants for 
alternating magnetic fields by the corresponding ratio of susceptibilities, with the 
result shown in the second row of the table above. These figures are the slopes at 
a number of points of a curve connecting mechanical stress with the steady mag- 
netising field producing it, and this curve can be found by graphical integratail 


; ae Nara, Sci. Rep. Tohoku Univ. 16, 338 (1927) 
pooner, Properti d Testi . : = 
ee perties and Testing of Magnetic Materials, p. 11. (The McGraw Hill Book Co. 
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from the slopes. The result is shown in the lower of the curves of Fig. 3. The upper 
curve is obtained when the changes of length of the same alloy as given by Masu- 
moto and Nara* are multiplied by Young’s modulus. The shapes of the curves are 
similar and they are roughly asymptotic at the stronger magnetic fields. The differ- 
ences at the weaker magnetic fields can be accounted for by a corresponding 
difference in the B/H curves. The alloy used by Masumoto and Nara was soft and 
that used by the author was somewhat hardened by cold working. Soft material 
magnetises more readily than hard material in weak magnetic fields and it should 
show a correspondingly greater magnetostriction effect in such fields. In strong 
fields, however, both materials tend to approach similar values of saturation in- 
duction and therefore should have similar magnetostrictive change of length. ‘The 
relation between the magnetostriction curves for the two samples is in fact qualita- 
tively the same as that between the B/H curves, and, moreover, the quantitative 
agreement at the higher magnetising fields is good. 
It may be concluded that the method of measurement which has been described 
is practicable and sufficiently accurate for the comparison of the magnetostrictive 
effects in alloys subjected to alternating magnetic fields. 
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DISCUSSION 


Dr J. H. Vincent. I think that the Physical Society is to be congratulated on 
the presentation of this paper, which marks a distinct advance in the subject of 
the application of magnetostriction to the production of mechanical oscillations. 
I suppose that the formula given for the period of oscillation of an anchor ring in 
its fundamental radial mode is only approximate. ‘The period of oscillation for an 

‘infinitely thin ring is the circumference divided by the velocity of sound in the 
- material. If the anchor ring is supposed to have its orifice filled up with the same 
material we approach the case of a wire vibrating radially. Whenever a current 
flows in a wire of magnetostrictive material the radius of the cross-section must 
change in consequence of the magnetic field which the current produces. ‘This 
should entail a number of interesting consequences, some of which might be capable 
of experimental demonstration. I should like the author to tell us a little more 
about the behaviour of the anchor ring when vibrating under the effect of alter- 
nating current in the absence of polarising current. 


Dr D. W. Dye. I have read the present paper with much interest and wish to 
mh congratulate the author upon a difficult piece of work in a field not hitherto explored. 
I admire the author’s courage in obtaining the results he has worked out, in spite 
of the admixture of paper, paraffin wax and nickel-iron alloy investigated. ‘These 


* Masumoto and Nara, loc. cit. 
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must form a specimen of rather poor vibratory properties. There are three points 
to which I would like to draw attention. (1) In discussing the Young’s modulus 
of the alloy the author states that Pierce should not have used the term aSB’ 
as a term to be subtracted from Young’s modulus. Pierce is quite consistent in 
what he has done; the difference between him and the author is simply one of 
definition of Young’s modulus, which is different in the two cases. One is free 
to choose either. (2) The dimensions of the built-up ring, as shown in the diagram, 
indicate a ratio of 1:3 between outer and inner diameters of the ring. This will 
result in a corresponding ratio of magnetising fields between inner and outer edges. 
The results obtained are, therefore, those referring to an average magnetisation 
corresponding to this range of H. (3) The values obtained from the circle diagram i 
and by the theory given pre-suppose a pure mode of radial vibration in which the 

stresses are purely compressional and tensile in the direction of magnetisation. — 
From my experience of quartz specimens of all shapes and sizes I should not be 
at all surprised to find that a ring of the shape shown had overtone flexural modes ~ 
such that the vibratory energy was divided up between two or more modes of ! 
vibration. This would profoundly affect the magnetostriction constant as here | 
determined. In the analogous case of the quartz plate, the effective piezo-electric 
constant nearly always comes out at about one-third only of its theoretic value when 
the motion is complex, simply because energy is divided up between various modes, 
only one of which reacts piezo-electrically upon the maintaining circuit. I would 
suggest that a ring of greater thickness should be tried to see if the same result 
is obtained. ; 


Mr W. A. Benton. I consider the paper of quite outstanding interest and 
importance. There was one remark, however, made by the author in presenting his 
paper to which I felt exception might be taken. Mr Smith spoke of the “alternating 
current producing stresses unaccompanied by changes of dimensions in the magnetic 
material”. It is difficult, if not impossible, to conceive of any stress not accompanied 
by strain, and, as a matter of fact, the mathematical reasoning would appear to hold 
good if small dimensional changes accompany the alterations in stress. It may 
interest Fellows to know that over 50 years ago Dr Du Moncel, and, I believe, 

' others, described telephones based on magnetostriction effects. A rod of magnetic 
material attached to the centre of a diaphragm reproduced music and speech when 
acoustically modified currents circulated round the rod. It would not appear, 
however, that anyone has used the reverse, or Villari effect, as a telephonic trans- 
mitter. ‘The very considerable changes of stress which could be produced by 
acoustic vibrations in a rod or wire made of a magnetic material especially susceptible 


to the Villari effect could be used to modify the strength of a current and thus 
form a transmitter. 


AutHor’s REPLY. The simple formula given in the paper for the fundamental 
mode of radial vibration of a ring is, as Dr Vincent points out, only approximate. 
An exact calculation which takes into account the finite radial width of the ring 
shows that the error involved in this approximation is negligibly small for the 
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comparatively thin ring used in the experimental work. I am unable to give 
Dr Vincent any further information on the behaviour of a ring subjected to an 
alternating magnetisation only. 


I agree with Dr Dye that the theory given by Pierce for the vibrating rod is 


“accurate provided that Young’s modulus for the material is suitably defined. Pierce 


does not, however, define this modulus in his paper, and the usual simple definition 
is therefore implied, whereas another rather complicated definition is required if 
his theory is to be consistent. There is no doubt in my mind that the predominant 


mode of vibration of the laminated ring is radial. The sound in air-emitted by the 


vibrating ring is audible and it is possible to explore its surface with a small 
stethoscope. The sound intensity appears to be uniform over the whole of the 


outside and inside cylindrical edges. The flexural forces in the vibrating quartz 


disc arise from the fact that the diameter of the disc is usually large in comparison 
with the wave-length of sound in quartz at the frequency of the vibration. The 
expansions and contractions of the quartz in the plane of the disc, depending upon 
Poisson’s ratio, are partially prevented by a lateral constraint from the inertia of 


the disc, and radial stresses are produced. These stresses are not quite symmetrical 


with respect to the quartz disc because the piezo-electric properties are not, in 
general, quite uniform throughout the disc, and stresses are therefore produced 
which can excite flexural modes of vibration. In contrast to the quartz disc, which 
is usually large in two directions in comparison with the wave-length, the laminated 
ring is small in two directions in comparison with the wave-length. Expansions 
and contractions in the directions of the thickness and the radius of the ring, 
depending upon Poisson’s ratio for the material, are practically unconstrained. 
Moreover, the ring is built up from a large number of thin stampings assembled 
at random, and any appreciable variation of the magnetostriction constant along 
the circumference is improbable. It is difficult to imagine how any stresses can 
arise which might excite a flexural mode of vibration. 

The conception of a material in a state of stress without strain, which Mr W. A. 
Benton finds difficult to accept, can perhaps be most clearly realised in the following 


case. Let a long thin rod of a material, having the property of magnetostriction, 


be magnetised so that its length and cross-section change. The rod can be restored 
to exactly the original dimensions by forces applied to its ends and surface, so that 
it is completely freed from strain. These forces are in equilibrium with the forces 
arising from stresses in the material, so that a state of stress without any strain 
whatever is produced. 


192 


FIELD INTENSITY MEASUREMENTS AROUND 
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ABSTRACT. The field intensities of wireless transmitting stations have been measured 
in order to determine some of the factors controlling the propagation of electromagnetic 
waves, and the effect of the surface conditions on their attenuation. A simple loop, 
condenser and valve voltmeter circuit has been employed, and careful tests show that 
this method is available for the measurements of intensities as low as 1 mv./m. The grid 
rectification voltmeter described is found to measure e.m.f.’s down to o-or yolt, to be 
free from frequency error, reasonably stable in calibration, and not to absorb power 
from the circuit to which it is connected—the last being contrary to current theory. 
The advantages of this intensity set, which agrees over a wide range of field intensities 
with a Western Electric set, in respect to simplicity, convenience and portability are 
pointed out. 

The field strength contours of three broadcast stations have been determined, from 
which the following conclusions have been drawn: (i) Very rapid attenuation of the 
signal is caused by Australian forest areas; this enormously curtails the areas for which 
a satisfactory service is provided. (11) The effective conductivity of the various types of 
ground surface met with varies from 4 x 10-18 to 0-07 x 107 e.m.u., according to the 
nun er of trees in the areas covered. (iii) The use of a longer wave-length gives a marked 
increase of intensity at distant points beyond forest areas. (iv) For daylight transmission 
over sea water up to a distance of 85 miles, after the application of curvature-corrections 
to the intensity, Sommerfeld’s formula is correct, to within the limits of experimental error. 
(v) The efficiency of radiation of the three aerials examined ranges from 48 to 60 per cent. 

Finally, the limits of signal intensity required for satisfactory reception under Australian 
conditions are discussed, and it is suggested that atmospherics and other disturbances 
are less prevalent in Victoria than in Europe or America. 


§i. GENERAL DESCRIPTION OF THE METHOD AND APPARATUS 


F we have a loop aerial tuned by a condenser, then the ratio of the signal intensity 
] to the voltage which it induces across the loop is constant for all.values of the 
intensity. ‘The value of this ratio depends on the geometry of the loop, its 
resistance and reactance, and can be determined for a given loop once and for all. 


Hence if this ratio is known, and if the induced voltage is measured at any point, 
the field intensity at that point can be deduced. 


In the set used by the author the centre point of the loop was earthed, as shown — 


in Fig. 1, in order to preserve the directional effect of the loop aerial by elimination 


of the “antenna-effect,” and the induced voltage across one-half of the loop was — 


na 
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measured with a valve voltmeter. The absolute value, E millivolts per metre, of 
the field intensity, was determined at this laboratory by an independent method 
and the voltage V across half of the field loop was measured at the same Bike. 
Hence the ratio E/V for this loop was obtained; it was found to remain constant 
by similar comparisons carried out at regular intervals. Its value was 


E/V = 86. 
The valve voltmeter was calibrated after each reading of V in the field by means 


of a high-frequency oscillator, a vacuum thermocouple and a non-inductive 
resistance box. ; 


Fig. 1. Circuit diagram. 


The measuring and calibrating equipment was mounted on a board in the back 
of a small motor car, while the fixed loop was carried in a case fastened on to the 
running board, and was erected on the side of the car for the taking of observations. 

“Measurements made with the apparatus in the car, and on the ground with the car 
removed, showed that the car had no appreciable disturbing effect. 


§2. DETAILS OF THE APPARATUS 


The valve voltmeter. The circuit is shown in Fig. 2. The points of interest of 
the voltmeter are: (a) the absence of the artificial grid leak resistance, so that the 
grid of the valve is apparently “‘free.”” This greatly increased the sensitivity of the 
instrument and diminished the power absorbed. There was no sign of instability 
of the valve, probably because the conductivity of the base of the valve, the valve 
socket and the valve pinch prevents the grid from blocking. (b) The size of the 
grid condenser, viz. 0°002 uf. (c) The balancing-out device for eliminating the 
steady plate current from the microammeter. This affords a very sensitive means 
of adjusting the zero of the instrument to any point of the scale. (d) The loop 
tuning condenser; as neither set of plates was at earth potential, in order to avoid 
body-capacity effects it was necessary to enclose the condenser and voltmeter in an 
earthed shield connected to the middle point of the loop. Although the batteries 
and the tuning condenser handle were not shielded, there was no disturbing hand- 
"capacity effect during tuning. The condenser was a double-spaced. transmitting 
‘condenser of maximum capacity 85 pp. 

The small link between the earth and grid terminals of the voltmeter was 
found to be of great practical convenience when frequent interchanging of the 
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connections was necessary. The grid terminal was very susceptible to large a 
tions of potential while connection was being made, causing momentary paralys: 
: fe aero of a voltmeter. These are: (1) sensitivity and range, (ii) frequency 
error, (iii) wave form error, (iv) power used, and (v) stability of eRe a _ 

The relative importance of these characteristics will depend on the circur ) 
which the voltmeter is connected and the quantity that is to be measured. When — 
a small high-frequency current is to be measured by means of the voltage <a 
an inductance tuned to resonance with the applied e.m-f., it is essential that “7 
voltmeter should absorb no power from the circuit. In the measurement of field . 
intensities by the method outlined above, any power taken by the voltmeter will | 
appear as an increase in the resistance of the loop. 


Earthed Shield 


Fig. 2. Valve voltmeter circuit. 


Performance of the grid-type valve voltmeter. (i) Sensitivity. With a portable 
microammeter having a sensitiveness of 10 mm. per microamp., the lowest measur- 
able voltage was o-o1 volt (deflection about 2 mm.) corresponding to a field 
intensity of o-86 my./m. This sensitivity was obtained by the correct choice of 
valve (Mullard P.M. 1 H.F.), the operating voltages, and the absence of an added 
grid leak. With a grid-leak resistance of ro megohms the sensitivity was con- 
siderably diminished, especially for small voltages where the extra sensitiveness was 
most required, (ii) Frequency error. If the voltmeter is to be independent of 
frequency, the reactance of the grid condenser must be small compared with that 
of the grid-filament capacity of the valve and its socket. Taking the latter as not 
greater than 30uuf., the grid condenser was fixed at 0-002 uf. To check this point, 
the voltmeter was calibrated at 600 kc. and goo ke. and the agreement was within 
1 percent. In the field the calibration was carried out at 750 ke., which is about the 
middle of the broadcast band under investigation. If the grid condenser was 
larger than o-orf., the grid took an appreciable time to reach its final steady 
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potential under the action of the incoming signal, and there appeared to be no 
advantage gained by the use of a larger value than that adopted. (ili) Wave form 
error. In order that the voltmeter reading should be independent of the wave form 
of the applied e.m.f., it is necessary that the rectified current be proportional to 
the square of the applied voltage. For the calibration curve below o-1 volt this 


condition is fulfilled, as can be seen from Table 1, which is taken from a typical 
curve. 


Table 1 
Rectified : J 
current ae O67" 1°45 | £2°75 | 43 645 | 84 | 11-2 «| 17°71 pa. 
Applied 
parenaal (opens i o'02 0°03 0°04 0°05 0:06 0:07 0:08 or volt 
Current Tes. | 40-8 | 40: ; 
(potential)? 2 ASS 405 ALS aes 414 | 419 |41°4 


Above o-r volt the curve departs from its parabolic form and it becomes 
straight from o-15 up to 0-4 volt. Hence for applied voltages below o-1 volt, the 
calibration is independent of wave form, while above this point it will depend on 
wave form. This characteristic is important when the induced voltage due to a 
modulated signal is being measured. Above o-r volt, there were considerable 
fluctuations of the voltmeter reading due to the modulations of the signal, while 
below this figure there was practically no sign of such modulation, even though 
the deflections of the microammeter were as large as 80 mm. As only the carrier- 
wave intensity was required, it was necessary to wait until the modulation ceased. 
The short period of the instruments used made this cessation easy to detect, a silent 
period of 2 seconds being enough to allow the pointer to take up its position for 
the unmodulated signal. No special precautions were taken to get a pure wave-form 
for calibrating the voltmeter except that the coupling between the oscillator and 
the tuned circuit containing the thermocouple and resistance box was kept loose. 
(iv) Power used. As was mentioned above, the power absorbed by the voltmeter 
should be small and independent of the applied voltage. A valve which uses the 
curvature of the grid-current grid-potential curve for rectification must absorb 
some power, because the action of the rectifier depends on the presence of grid 
current. In broadcast receivers, the grid condenser is usually connected to the 
positive terminal of the filament, with a resistance of 1 to 3 megohms across the 
‘condenser, and such an arrangement increases the resistance of the input circuit 
enormously, necessitating the use of reaction to counteract it. This damping of the 
input circuit increases as the input voltage decreases, and largely counterbalances 
the gain in sensitivity due to the grid condenser and leak. In the voltmeter 
described above these defects are reduced to a minimum, while the sensitivity is 
increased, by the omission of a grid leak. This arrangement is unsuitable for the 
‘detector valve in a broadcast receiver. No appreciable effect of the voltmeter on 
the input circuit could be detected by the following test. ‘T'wo instruments of 
identical construction were connected in parallel across a loop tuned to resonance, 
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and a small voltage was induced in the loop by a local oscillator. Upon the switching 
out of the filament in one voltmeter, no change was detectable in the reading of the 
other meter. This test repeated for a series of voltages gave no evidence of power 
absorption. This only shows that the grid voltmeter has the same effect on the a 
circuit whether the valve is alight or not, but bridge measurements at 1000 ~ 
showed that in both cases the input of the voltmeter with tuning condenser behaved 
as a nearly perfect condenser of about 30 upf. capacity. The same test, repeated 
with an anode-bend valve voltmeter and the grid voltmeter, showed that the | 
former instrument absorbed power. Evidence to the same effect had previously ~ 
been obtained by other means. (v) Stability of calibration. The voltmeter was 
calibrated after each reading of the voltage, so that stability of calibration, although — 
it obtained fairly well, was never assumed. The variations during a day were | 
somewhat irregular, but the constancy was usually within ro per cent. The cali- 
bration curves for a single valve over a period of several months were all of the same 
form, showing an equal percentage change over the full range. Hence one or two 
points on the curve served to fix the whole of it with sufficient accuracy for the 
deduction of the signal intensity. 

The microammeters. For applied voltages above 0-08 volt, a Paul unipivot. 
galvanometer (1-8 a. per division) was used; this instrument when shunted also 
measured the steady plate current. For lower voltages, a portable reflecting 
Tinsley galvanometer (10 divisions per ja.) was used. The shortness of the period — 
of these instruments made them very suitable for the work. These two galvanometers _ 
were mounted together in a box which was suspended by six small coil springs and 
rested on rubber sponges inside another box. The efficacy of this suspension was 
shown by the fact that in over 5000 miles of travel, for a large part on rough and 
unformed country roads, no trouble was experienced with this reflecting instrument. 

No special mounting was necessary for the rest of the apparatus. 

The microammeters showed no instability due to the absence of an artificial 
grid leak, but there was usually present a gradual creep of the zero of the voltmeter 
if the batteries and valve had not reached their steady state. No difficulty was 
experienced in allowing for this creep of the zero, unless it exceeded the rate of 
I mm./sec. : 

Calibration of the voltmeter.'The voltmeter was calibrated by a potential measured, : 
as in Fig. 3, by means of a thermocouple and a non-inductive resistance box 
calibrated with high-frequency current. Care was taken to minimise stray capacity 
effects and good agreement was obtained between a high-frequency calibration and : 


a 50~ calibration of an anode-bend voltmeter designed to be independent of — 
frequency. 

The field loop. In order that the ratio of field intensity to induced voltage, i.e. 
E/V, may remain constant, it is essential that the loop aerial used in the field should : 
retain its shape and size, and also that its high frequency resistance should remain — 
constant. While it was desirable to reduce this as far as possible to get a more 
efficient aerial, it was even more important to secure permanence, even at the — 
expense of efficiency. It was therefore decided that solid copper wire should be — 
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used instead of Litz wire. The loop was a square solenoid (side 65-1 cm.) of 18 turns 
of no. 24 gauge enamelled copper wire, spaced 0-2 in. apart on thin bakelite strips. 
This construction gave an exact centre point, and should diminish any possibility 
of change in resistance during the field work which was carried out during fine 
weather. In the course of 3 months no evidence was obtained for any change as 
great as 2 per cent. An occasional tightening of the loop was the only attention 
required. 


Se 


Fig. 3. Calibration circuit, 


§ 3. DETERMINATION OF THE RATIO E/V FOR THE FIELD LOOP 


The relation between the field intensity £ mv./m. and the voltage V millivolts 
induced between the ends of a tuned loop is 


V = 20AnwLE|RX 


where A is the area of each turn in metres’, 
n the number of turns, 
\ the wave-length in metres, 
w/27 the frequency, 
L the inductance in henries, 


R the resistance in ohms. 


All of these quantities can be measured. From the above relation it would 
seem to be possible to get the value of E/V for the whole loop, and hence for half 
the loop, but the voltage across the half-loop was found to be more than half that 
across the whole loop. This result was due to the circumstance that the natural 
frequency of the loop was comparable with that of the stations under investigation, 
and difficulties were experienced from the same cause in the measurement of the 
high-frequency resistance of the loop. While approximate corrections for this 
self-capacity effect could have been calculated, it was decided to measure the 
‘absolute value of the intensity with another loop, at the same time measuring the 

voltage across half the field loop. Hence the ratio E/V was at once obtained. It 
was E/V = 86. This ratio was checked at intervals during the field work by the 


same method. 


198 R. O. Cherry 


§4. DETERMINATION OF THE ABSOLUTE VALUE OF THE 
FIELD INTENSITY 


A standard loop was constructed for this purpose and was kept at this laboratory, 
the intensity being calculated from the voltage induced across the whole of the loop 
as measured by the valve voltmeter. It was a square (side 132-3 cm.) of 4 turns 
of & in. woven copper braid, spaced at 4 cm. apart. The inductance and resistance 
were determinéd separately and the value of wL/R was deduced for a wave-length — 
of 371 metres. This ratio was checked by the capacity-variation method. A special 
set of resistances was constructed and the measurements were repeated with both 
the voltmeter and a vacuum thermocouple as the detector in all determinations of 
R, L, or wL/R. The following values were obtained: L = 51-9 microhenries, 
R = 0913 ohms when A = 371 metres, so that wL/R = 287. 

By capacity variation wL/R = 291. The mean value of the step-up is therefore 
289, the high value being due to the low resistance of the copper braid. With this — 
loop, the absolute value of the intensity is considered to be accurate to 1 per cent. 

To check the absolute value of the intensity as determined with the standard — 
loop, simultaneous readings at the University were taken with it and a Western 
Electric field intensity set, the property of the Research Laboratory of the Post- 
master-General’s Department. The mean value of the intensity determined by 
both sets was 55 mv./m. The standard loop was similarly calibrated for a wave-length 
of 484 metres, corresponding to that of the other local broadcast station, 3AR 
Melbourne. A comparison with the Western Electric set gave an intensity of 50, 
as against 51 mv./m. read by the standard loop. The ratio of E/V for the field loop 
for this wave-length was E/V = 107. 

To test if #/V for the field loop was independent of EF, a series of readings over 
a wide range of intensities was taken in conjunction with the Western Electric set. 
With the value E/V = 86 obtained above, Table 2 gives the results obtained : 


Table 2 


4 
| 


W.E. set 59°5 30°5 13°4 | 10°4 75 2°8 2°15 1°46 | 0-55 mv./m. 


Valve volt- 
meter set 59°5 30°5 I4'0 


10-4 76 370 | 2°4 


It therefore appears that this method of measuring intensities is as accurate as 
is required in this class of work, 


§5. MERITS OF THE VALVE VOLTMETER FIELD INTENSITY 
MEASURING SET 


The apparatus described will measure satisfactorily field intensities as low as 
1 mv./m. This is a field too small for satisfactory reception according to English 
and American authorities. 

Advantages. The following advantages are claimed for the set: (i) Simplicity 
of construction. This is apparent from the details already given. The voltmeter and 
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oscillator are made of inexpensive wireless equipment, and the high-frequency 
decade resistance box, vacuum thermocouple and galvanometers are all easily 
obtainable. No elaborate shielding, such as is essential on sets making use of 
loop-injection methods, is required, and considerable weight is thus saved. Other 
difficulties, e.g. the attenuation box or the carefully balanced input circuit, are also 
avoided. (ii) Simplicity of operation. The loop tuning condenser is the only 
tuning control, and with a vernier drive this is easy to operate. ‘The voltmeter 
deflection for several broadcasting stations can be quickly determined, and a single 
frequency of calibration serves for all of them. This is an important simplification 
for rapid and accurate work. (iii) Speed of work. The whole of the field work, 
‘including driving, observing, booking and reduction of the results, was done by 
the author. The observations used for plotting the contour map were obtained in 
25 days in the field. In this period 497 individual observations were taken on two 
stations during their normal daylight working hours, no special transmissions 
being arranged. The area surveyed was over 10,000 sq. miles, and the distance 
travelled was 3067 miles. The time taken in an observation on two stations was 
about 7 minutes, from stopping to starting again. During the 9 months that the 
apparatus has been in use, not'a single delay due to failure of any part of it has been 
experienced. (iv) Portability. The loop and voltmeter alone can easily be carried 
by one man without a car, and this enables the apparatus to be set up in places to 
which a car cannot get. This property enabled the detailed effects of hills on the 
field intensity to be studied. 

Disadvantages. The great disadvantage of this type of set is its inability to 
measure very weak intensities. The rectifier used should obey a square law, Le. 
the rectified current should be proportional to the square of the applied voltage, 
which holds for this instrument below o-1 volt. For an intensity of 1 my./m. the 
deflection is about 2:5 mm., and half of this intensity would only give 0-6 mm. 
By slight alteration of the construction of the loop, however, without increase of 
its dimensions, the ratio E/V was reduced from 86 to 63 and a slightly more 
sensitive detecting valve (P.M. 12) has been added. The deflection for 1 mv./m. 
is now 7mm. No simple way of increasing the range of the instrument further, 
while still retaining its portability and other advantages, has been found. As there 
is only a single tuned circuit, lack of selectivity would also be a difficulty under 
certain conditions. 


§6. ROUTINE IN THE TAKING OF AN OBSERVATION 


For taking an observation, open sites likely to give the undisturbed value of 
the intensity in that locality were always chosen. The chief disturbing factors 
were: (i) Telephone wires, especially if present in large numbers. In one case a 
bank of 50 wires diminished the field intensity for over 50 yards on each side of 
“the line. In another case, two wires doubled the intensity directly under them 
and produced a directional error of nearly 50°. The field was normal again 15 yards 
on either side of the wires. (ii) Crests of hills. On a steep slope facing toward the 
broadcast station, the intensity was always greater than at the foot of the hill or 
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ce behind the crest. (iii) In timbered country, an open place - er 
so yards from any large tree was always selected. The site having been se a + 
the car is turned side-on to the direction from which the signal is Semin . 
valve, which is carried separately, is inserted into the voltmeter, and the ona 
set up. The loop is tuned and rotated to give the maximum deflection on ; 
microammeter, which is noted. The voltmeter terminals are then short circuite: 
with the link, the loop terminals are disconnected, and short leads from the re- 
sistance box are connected. The link is opened, the local oscillator started and the 
current is adjusted to give a convenient reading on the thermocouple metre. The 
voltmeter deflection and thermocouple currents for known resistances are noted, 
and the calibration curve is plotted. In the field this curve is only taken twice a day. 
For most of the observations, a single point near the observed voltmeter reading or 
two points on either side of it are sufficient, on account of the similarity of all the 
curves, as has already been explained. From the calibration curve the voltage due 
to the incoming signal is found and the intensity is at once deduced. 
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§ 7. RESULTS OBTAINED IN VICTORIA ON 3LO AND 3AR MELBOURNE 


Field intensity map of 3LO Melbourne. The observed intensities were plotted 
on a map, and the 15, 10, 5, 2-5 and 1 my./m. contours were interpolated, as 
shown in Fig. 4. Only three observations out of 330 that are plotted are not in 
agreement with the contours as shown; two of these were suspected at the time 
of observation to be affected by ground contours but they were taken at the only 
available sites. 

Effects of trees and mountains. The shape of the contours indicates that very 
marked absorption of the signal takes place in certain directions. In every case, 
these areas are heavily timbered and are usually associated with mountains or 
steeply sloping hills. In Victoria mountains and forests usually occur together, 
and it has not been possible to discriminate completely between them as to the ; 
cause of the absorption. In a few cases where the timber is dense while the country 
is only undulating, the rapid absorption has still been found, but it is not possible 
to state the effects of bare mountainous country. Only in a few directions does the 
transmission take place over fairly clear ground and the area given a satisfactory 
service is far less than was expected from the extent of the 10 my./m. contour which : 
was located before the present apparatus had been developed. In heavy timber 
the intensity falls from 10 my./m. to 1 my./m. in a distance of only 10 miles. 
Such areas exist to the east behind Mount Dandenong, to the north-east behind 
Whittlesea, and to the north-west in the Daylesford forest. Barfield in England* 
found that the signal was attenuated more rapidly when it passed over forest — 
areas than when it passed over clear country, but in no areas were the effects so_ 
pronounced as in the present observations. This discrepancy is due to the difference _ 
in size and the number per unit area of the trees, in the two cases, for in the areas 
dealt with in the present paper the forests cover many hundreds of square miles, 


* Journ, Inst. E, E. 66, 204 (1928) ; 67, 253 (1929). 
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while the density of the trees usually exceeds 10,000 per km?. In many cases the 
number is probably greater than 30,000. 

The height and leafiness of the trees was found to play an important part in 
the attenuation. Thus, in the areas mentioned, the trees are mostly over 100 feet 
high and have a fairly luxuriant foliage. In contrast to this, there is an area about 
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_ Fig. 4 Field intensity contours of 3LO Melbourne, wave-length 371 metres, power 5000 watts.* 
Shaded areas represent mountains and forests. 

* This is the power in the plate circuit of the power amplifier coupled to the aerial. 
LO covered with small trees, between 20 and 40 feet high, 


little foliage. The effect of this belt of timber is plainly 
t is not so marked as in the 


40. to 50 miles west of 3 
and with comparatively 
seen in the form of the contours in this direction, but i 
former areas. This very rapid attenuation of the signal as it passes over forest 
areas explains the difficulty experienced in receiving Melbourne broadcast stations, 
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in the eastern parts of Victoria, where, during daylight, the signal is very weak but 4 
of uniform intensity, while, after dark, fading and distortion are so bad as to spoil — 


the programmes completely. These reception difficulties start at Warragul, about | 
60 miles east of Melbourne, where the daylight intensity is considerably less than 


o-s my./m. The conditions existing in Gippsland are exactly similar to those in 
certain areas around New York City and investigated by Bown, Martin and 
Potter*. In the latter case the rapid absorption due to steel-framed buildings in 
the city near the transmitter produced areas within 35 miles of the station in which 
the daylight intensity was as low as 50 pv./m., and within these areas fading and 
distortion were very pronounced at night. It now appears that large areas of 
mountainous and timber country can produce exactly similar conditions at com- 
paratively short distances from the station. 

The land area within the 1 mv./m. contour is 9400 miles*. If the ground had 
all been of average conductivity, so that the 1 my./m. contour was a regular curve 
with its limits at the existing contour to the south-east, north and south-west, the 
land area within it would be over 20,000 miles?. 

Determination of the ground conductivity for various surface conditions. This can” 
be deduced from the field intensity contours by applying Sommerfeld’s theory of - 
ground absorptiont. The convenient curves published by the British Broadcasting 
Corporation{, which give Sommerfeld’s formulae in graphical form, make it 
possible to write down the intensity at a point distant from a transmitter, given” 
the ground conductivity; or, as in the present case, to reverse the process. No- 
curvature corrections have been applied, and it is considered safe to use the curves 
up to 100 km, only. The three directions in which the ground conditions are 
approximately uniform are: (1) South-west. The country is as uniform as could 
be desired, as there are but a few low hills and very few trees. The five contours 
(15, 10, 7°5, 5 and 3:mv./m.) agree excellently with those predicted from Sommer-_ 
feld’s theory for a conductivity of 4 x 10-!* e.m.u. or a resistivity of 2500 ohm-cm. _ 
This value is taken as the true “ bare ground” conductivity for Victorian conditions. 
The corresponding value determined by Barfield and others for English conditions 
Is 2°5 x 1078 e.m.u. (ii) South-east. The country is flat but there is a good deal 
of low scrub and some light timber is present in parts. The ground for a large part. 
has a tendency to be swampy in winter. The five contours (15, 10, 775, 5 and= 
2'5 mv./m.) agree fairly well for a conductivity of 3-3 x 107! e.m.u. (iti) North. 
The country tends to become hilly and timbered at Kilmore and Broadford, 
Ss miles north of 3 LO. There is however a narrow pass in which the signal is not 
rapidly attenuated. The country then becomes flat again, but a considerable amount 
of scattered timber still remains, a rough estimate giving from 200 to 500 trees 
per km.° The country is all comparatively dry, and the conductivity deduced from 
the 10, 7°5, 5 and 2:5 my./m. contours is 10-8 emu. In heavily timbered and 
mountainous country, the effective conductivity deduced is as low as 0-07 x 10719 


Proe. Inst. Rad. Eng. 24, 57 (1926). t+ Annalen der Physik. 28, 665 (1909). 
t The Service Areas of Broadcasting Stations, by P. P. Eckersley. 
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é€.m.u., which is only 1/50 of its value under favourable conditions. The value 
given by Eckersley* for this type of country is 0-05 x 10o-P e.m.u. 

All these values are deduced from the contours which are interpolated between 
readings of the intensity taken at suitable points. More accurate determinations 
could be made from special observations taken at selected stations, with checks on 
the station output, but in only a few directions is the ground of sufficient uniformity 
to justify such an undertaking. In the most favourable direction, viz. south-west, 
the five points taken along a radial line agree excellently with the value stated. 

Recent determinations of ground resistivity made by Mr Bruckshaw, of the 
Imperial Geophysical Experimental Survey, by Gish and Rooney’s method give 
the value 1500 to 2000 ohm-cm., which agrees satisfactorily with the “bare ground” 
conductivity determined above. 
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Shielding effects of the masts at the transmitting station. In the field strength 
contours around 2LO London, as shown by Barfield and Munrof, large crevasses 
occurred owing to absorption in the masts supporting the transmitting aerial. There 
is no sign of similar effects in the case of 3 LO Melbourne. The masts at the latter 
station are 200 feet high and are earthed at their bases. The aerial is in the form 
‘of T' with the horizontal “roof” pointing to the north. 

Effect of change of wave-length on the intensity at distant points. Information 
on this matter was obtained from observations on 3AR Melbourne on a wave- 
length of 484 metres. This station was “on the air” only during the morning, and 
results for it are not so complete as for 3 LO Melbourne. ‘The general form of the 
contours is fairly reliable, and Table 3 summarises the results of 160 observations. 
In it is shown the decrease of the product intensity x distance with increase of 
distance for various types of country. For purposes of comparison, all the values 
are so reduced as to make the product exactly 100 at the shorter distance. 

The decrease of Ed in the above table will partly be due to the curvature of 
‘the earth, the effect of which at more than 50 miles from the transmitter is appre- 
ciable. It has not yet been possible to separate the losses into true-absorption and 


curvature losses. 


* Eckersley, loc. cit. p. 8. + Barfield and Munro, loc. cit. 
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The aerial system at the station 3AR Melbourne does not even rea 
to a vertical uplead with a horizontal roof, there being a certain ee ) enetem 
directed upwards at an angle of from 30° to 60° in certain directions. e observa- 
tions contain evidence that some of this energy 1s returned at distances 7 
than 4o miles, and so it is not Pei to apply Sommerfeld’s formulae to deduce 

nductivities from this station. 
i Table fice the decrease of intensity due to a small extent of heavily 


timbered and hilly country. 
Table 4 


: Decrease of intensity 
Width of timber 
in miles A = 371 metres A = 484 metres 
6 12°4-1°6 my./m. | 14°5-5°7 mv./m. 
is 86-20. 5, | 8-0-3°6 Sy, 
6 1a seh eeos 66-411 ,, 
13 10°5-2 a / 105-376, 


All the quantitative values for the conductivity of the various types of ground : 
surface dealt with above and the variation of the attenuation with wave-length 
are deduced from values of the intensity extending over a wide area of country, in 
order to obtain the general features of the transmission. More detailed investigatio: 
of these values will be undertaken. 

Transmission over sea water. In order to investigate the decrease of intensity 
due to the effects of the curvature of the earth alone, a series of readings was taken 
over sea water where the surface has a uniform high conductivity (10! e.m.u.). 
The readings were taken during the daylight part of the journey from Melbourne 
to Burnie, ‘Tasmania, which extended to go miles from Melbourne. The motion of 
the boat made the readings of the unipivot galvanometer uncertain of 1 my./m. in 
all the readings at distances of more than 40 miles. A large number of readings 
was taken and the time for each was recorded. The position of the ship at that 
time was computed from its measured speed over a distance of 30 miles before 
entering open water. The results are marked with centred circles in Fig. 5, in which 
log Ed is plotted against d. 

If the intensity is corrected for curvature from results given by Ratcliffe and 
Barnett* for similar wave-lengths, the points marked with crosses are obtained; 
they lie about a horizontal line. This result indicates that there is no attenuation 
due to the finite conductivity of sea water up to a distance of 85 miles, and is in 
agreement with Sommerfeld’s formula for transmission over a plane surface of 
conductivity 10! e.m.u., which gives the departure from the horizontal line at 
that distance as less than 2 per cent. To test this point further, readings were 
taken on 3 LO and 3AR on both sides of Port Phillip Bay at the edge of the water. 
Simultaneous checks on the station output were obtained by separate intensity 
measurements at a fixed point, and large-scale ordnance maps were used to make 


* Proc. Camb. Phil. Soc. 28, 288 (1926). 
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sure that the measuring sites were in line with the transmitter. Three separate 


lines from each of 3LO and 3 AR were taken, and the length of sea water varied 
from 31:6 to 35:1 miles. 


Log (intensity x distance) 


@°%500 ZR ae eee <—=¥ ae x 
fe} 
fe} ° 
° 
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Fig. 5. Transmission over sea. 


In every case the product Ed for A = 371 metres showed a decrease, the mean 
being 5 per cent., while for A = 484 metres, the product increased, the mean value 
of the increase being 4 per cent. This increase is attributed to the return of the 
energy initially radiated upwards from the aerial. It may be mentioned here that 
a slight increase in the product Ed with a small increase in d has been found 
frequently, even in the case of flat land transmission, for both stations. 


§8. RESULTS OF OBSERVATIONS ON 7ZL HOBART 


The wave-length of this station is 516 metres and its rated power is 800 watts. 
The field-intensity contours for 10, 7°5, 5, 3 and 1 my./m, are shown in Fig. 6. 
These are drawn in from about 140 observations. Nearly the whole of Tasmania 
consists of wooded hills and mountains, and the effect of these is seen plainly in 
the form of the contours. The effect of Mount Wellington (4000 ft.) to the west of 
Hobart is very marked. Only to the north of the station does the transmission 
‘take place over undulating country without very much timber. The conductivity 
of the ground in this direction is = 0°5 x 107% e.m.u. Tasmanian observations are 
in complete agreement with the conclusions stated for above Victoria. 
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Fig. 6. Field-intensity contours of 7 ZL Hobart. Wave-length 516 metres, power 800 watts. 


§9. EFFICIENCY OF THE AERIALS OF 3L0, 3AR AND 7ZL ; 


‘The intensities at points on a circle of small radius d km. from the transmitter — 
are averaged, giving the mean intensity at that distance, for which loss due to — 
absorption will be negligible. If E volt/metre is this mean intensity, the power 
radiated by the aerial is 


P = E*d*105/9 watts. 


Table 5 gives this power and the corresponding efficiencies of the three stations 
investigated, 


gee ng ead 
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Table 5 
Distance Mean intensity | Power radiated | Input 
d km. volt/metre watts mat ee Efficiency 
BIO Ag 0'066 2680 5000 6 9 
3 AR 45! O'0104 2430 5000 43-6 é 
7 Ll, 20°8* oop Ke) 483 800 60:4 % 


1 'This is the mean radius of the 10 my./m. contour for these two stations, and a small 
ae for attenuation, based on a conductivity of 4 x 10°" e.m.u., has been applied 
to 3AR. 


§10. DISCUSSION OF STANDARDS OF SIGNAL INTENSITY FOR 
AUSTRALIAN CONDITIONS 

English and American authorities mention the value of 1o mv./m. of intensity 
as constituting a continuously reliable service. Also that reception becomes 
unsatisfactory below 2:5 mv./m. From information supplied by local people in 
the areas investigated during the present work, and from personal observation 
wherever possible, it appears that the service limits quoted above are unnecessarily 
high for local conditions. Crystal sets were found in use down to 2 my./m., while 
satisfactory valve reception was possible below 1 mv./m. Below 2 mv. /m., however, 
considerable variations at night were noted, although not sufficiently large to be 
classified as ‘‘fading” by the local listeners. The signal never became so small as 
to disappear entirely. 

The results of a series of transmission tests held by 3LO and 300 amateur 

listeners in Victoria, in February, 1928, showed that (i) fading, although not 
serious, was observed in many districts within 100 miles of Melbourne. The 
intensity at most-of these places has since been found to lie between 1-5 and 
0:5 mv./m. (ii) Within this area, atmospherics did not cause serious interference. 
(iii) Reception was unsatisfactory, in general, at places over 100 miles from 
Melbourne, except in Gippsland, to the east, where reception was only satisfactory 
up to 60 or 70 miles. The intensity at these distances is seen from the map to be 
less than 0-5 mv./m. 
We conclude therefore that the area of satisfactory reception extends to the 
-2*5 mv./m. contour and that no serious interference is present down to 1 mv./m. 
No measurements of the “‘noise level” or of the prevalence of atmospherics have 
yet been made in the south-eastern part of Australia, but the facts quoted above 
indicate that both these forms of interference are probably less prevalent here 
than in England or America. 
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DISCUSSION 


Mr R. A. Watson Watt (who read the paper in the absence of the author). 
I have an uneasy feeling that the author has made unjustifiably heavy concessions — 
to portability in his apparatus. In aiming at a high step-up wL/R for the field loop 
he has been driven to working the loop dangerously near its natural frequency, 
and to using a tuning capacity less than three times that of the voltmeter. I would | 
suggest that this capacity ratio has an additional disadvantage in vitiating thea 
attempt to eliminate antenna effect. The addition of 30upf across one half of a 
loop tuned by a capacity of the order of 85 yf should have been compensated by 
an equal capacity across the unused half. The discrepancy between half-loop and — 
whole-loop e.m.f.’s (p. 197) is probably due in part to this omission. It would, too, 
have been desirable to include provision for measurement of loop resistance in the” 
field, a provision which has been found necessary in most measuring sets. The 
very wide departure of the E/V ratios from the square of the wave-length ratio is 
disturbing in this connexion, for 107/86 = 1-24 while (484/371)? = 1-70. 

As regards the voltmeter, the unorthodox use of the valve cap and socket as a_ 
grid leak leads to a fear that the need for calibration at each measurement may 
have been mainly due to the variations of this composite resistance of large surface. j 
Further information on the gain in sensitivity as compared with a well-insulating _ 
socket and cap and a leak of, say, 1o megohms, would be of interest, as would also _ 
information on the mode of action of such a rectifier. Mr W. G. Hayman, of the 
University of Western Australia, who is at present attached to the Radio Research 
Station, Slough, has been good enough to make for me some tests on a voltmeter 
set up in accordance with the author’s figures. With a P.M.1.H.F. valve taken at 
random from stock, an arbitrarily chosen signal of frequency 1 megacycle per 
second gave a decrease from 0-7 ma. to 0-3 ma. in anode current with the author’s 
free grid arrangement. The grid pin was then isolated by the cutting of the valve 
cap while the socket was left unaltered, and the anode current change for the same 
applied signal was now from 0-6 ma. to 0-06 ma. With a 5 megohm leak across the 
grid condenser the current changed from 0-56 ma. to 0-065 ma. on application of 
the signal. It would thus appear that it is more important that the grid leak proper 


should not be swamped by casual leakage than that it should be superseded by 
such leakage. 


The grid current in the case last described was certainly under 0-02 pa. The 
mechanism of the moderately effective rectification obtained under these conditions, 
with 7,/e, and z,/e, characteristics as in the accompanying diagram, is not clear. 
In any case it is difficult to see how a grid-bend rectifier can absorb less power than 
an anode-bend rectifier properly designed for mitigation of Miller effect; that this 
is a grid-bend rectifier is shown by the decrease of anode current on application of 
signal. 

The phrase “field too small for satisfactory reception” (p. 198) should, of course, 
be qualified by the insertion of “domestic” before “reception”. 
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It might with advantage be pointed out that the relative attenuations on two 
wave-lengths shown in Table 3 are roughly consistent with Sommerfeld’s analysis 
and with Barfield’s conception of a “‘pseudo-ground-conductivity” o, which 
includes absorption by trees. The ratio of attenuation on the two wave-lengths over 


the same terrain is, in each case, close to the ratio of the squares of the wave-lengths. 


Sommerfeld’s ‘‘numerical distance” is, to a crude approximation, and for con- 
siderable distances, inversely as the square of the wave-length. This relation does 
not, however, hold in Table 4, which might well be expanded from its present 
somewhat obscure form. § 9 might be amplified to show that a specialised definition 


of efficiency is implied in the assumption that the polar diagram of radiation is a 


torus of zero aperture. 

I do not think that the author’s conclusion as to relative noise levels in Australia 
and in other places can be supported by the combination of psychological and 
physical evidence which he adduces. It is probable that the Australian listener is 
more tolerant in respect of noise level than is the hypothetical listener portrayed 
in the very stringent specifications of service field strengths used in Great Britain 
and America. 

These criticisms are, however, of but slight importance as against the great 
value of having such measurements, in another hemisphere, put before us in a form 
which shows very satisfactory agreement with measurements in Europe. As a 
worker for the Radio Research Board in Great Britain, I appreciate the privilege 
which I enjoy to-day of deputising for a colleague working for the Radio Research 
Board of Australia. 


Mr R. L. Smrru-Rose. Mr Cherry’s paper is interesting, both from the practical 
point of view in providing data which are useful to engineers responsible for 
the design of broadcast transmitting stations, and from its more scientific aspect 
in providing additional knowledge of the electrical constants of the earth at radio- 
frequencies. The method adopted for making the field intensity measurements 


"appears to depend, among other things, upon an accurate knowledge of the radio- 


frequency resistance of the receiving loop. It is perhaps doubtful whether this 
resistance can be measured with sufficient accuracy to warrant the statement on 
p. 198 that “the absolute value of the intensity is considered to be accurate to I EGE 
This criticism does not, however, in any way influence the value of the results 
recorded in the paper, since a much lower order of accuracy would be adequate 
to illustrate the results and permit of the various deductions drawn therefrom. 

I note that the conductivity obtained from the measurements made for open 


ground conditions in Australia +s of the same order as that obtained by Mr Barfield 


and myself in some measurements carried out in the south of England a few years 
ago, by an entirely different method*. This value of the conductivity is decreased 
when the ground is covered by trees, as was shown by Barfield’s later work which 


3 is referred to by the author. The much lower conductivity values obtained by 


Mr Cherry are undoubtedly due to the greater density of trees prevailing in certain 


* Proc. R. S. 107, 587 (1925). 
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parts of Australia. Similar results have been obtained by Lemoine and Rolf in 
the course of their investigations of the field strength from broadcasting in Sweden. i 
With reference to Fig. 5 of the paper, it may be pointed out that there appears ‘ 
to be a definite undulation in the course of the results there plotted. This undulation — 
appears to be sufficiently regular to warrant the dismissal of the idea that it is due 
to unsystematic experimental errors. No reference to this fluctuation is made in 
Mr Cherry’s paper, although I notice that in one of his previous publicationst he 
has drawn attention to this characteristic of the graphical representation of the 
results of field intensity measurements obtained for transmission over sea water. 
In the paper referred to the author suggested that the fluctuation was due to the 


interference of two waves, one of which was transmitted along the earth’s surface, — 


while the other was due to reflection at the lower boundary of the Heaviside layer. 
Since all the measurements were made in daylight, and it is known that the intensity 
of such a downcoming wave is practically negligible for the wave-lengths in question, 
it is unlikely that this will provide a true explanation. I should like to suggest as an ~ 
alternative that the second wave is obtained by reflection from the coast line on one 
side or the other of Port Philip Bay. At a distance of about 30 miles from the trans- 
mitting station, the ship on which the measurements were being made would be 
somewhere in the neighbourhood of St Leonards, and it is just possible that at 
this point direct and reflected waves would be obtained with a path difference 
sufficient to provide the ripple having a ‘‘ wave-length” of about 5 miles, as in- 
dicated in Fig. 5 of the author’s paper. At a distance of 50 miles the ship would be 
well outside Port Philip Bay, and it is likely that the distance between the points 
at which the two waves come into phase would be greater. This is consistent with 
the results recorded in Fig. 5, since the distance between the peaks increases to 


11 and 16 miles as the distance of transmission is increased from 50 to 80 miles. — 


The possibility of the reflection of waves at a coastal boundary has been pointed — 


out previously in connexion with some experiments on wireless direction-finding 


carried out in this country*, while confirmatory evidence has also been provided ' 


from the results of some experiments in Palestine. 

It is a well-known occurrence in wireless direction-finding that when the path 
of transmission lies nearly parallel to the coast line, a deviation of the waves may 
take place of the order of 3° or 4° on wave-lengths in the neighbourhood of 450 
metres. Where such refraction effects are observed it is reasonable to suppose that 
reflection may also take place. The effective coefficient of reflection will be small 
for the usual values of the conductivity of land and sea water, but this is consistent 
with the small amplitude of the ripple shown in Fig. 5 of the paper. 

If such experimental evidence of coastal. reflection could be confirmed by the 
carrying out of further experiments, a very useful purpose will have been served, 
since at the present time there is some discrepancy between the experimental data 
on coastal deviation and the theoretical explanation thereof. Perhaps the author 
¥ The Daylight Transmission of Wireless Waves over Sea Water (Broadcasting Company of 

ustralia, 1920). 
t Nature, 116, 426 and 408 (1925). } Nature, 121, 35 (1928). 
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could supply some further information or data obtained from his experiments 
which would assist in the further elucidation of this point. 


Dr D. Owen said that the last speaker’s suggestion was an extremely interesting 
one, but the effect of coast-line reflection, if real, appeared to be of small magnitude. 


The fluctuations of intensity shown in Fig. 5 formed only a small percentage of the 
total intensity. 


The PresIDENT said that wireless engineers would be very glad that this work 
was being done in Australia. Designers would be helped by the knowledge that 
‘the effects on transmission of soils and land surfaces are of the same order in that 
continent as in England. 


The author’s reply will be published in a subsequent part of the Proceedings. 
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A METHOD OF EXAMINING STEREOSCOPIC 
PHOTOGRAPHS 


By J. M. NUTTALL, D.Sc. anp E. J. WILLIAMS, Pu.D. 


Communicated by Prof. W. L. Bragg, October 22, 1929. 


Read by Mr T. Smith and discussed fanuary 24, 1930. q 
ABSTRACT. The paper describes a method of analysing stereoscopic photographs which 
has been found suitable for the examination of Wilson cloud tracks. The method consists” 
essentially in replacing the photographic plates in the cameras, illuminating them and 
tracing out the common image which coincides in space with the object originally photo- 
graphed. A system of movable pin points is used to trace out the contour of the im : 
the adjustments of the cameras being kept exactly the same as when the photographs w 
taken. The chief advantage of the method is its simplicity—practically no computation 
is required and it is not necessary to know the stereoscopic angle, the magnification, etc. 
The method has been used by the authors to measure the ranges of secondary 8-rays 

in investigations on the initial directions of emission of photoelectrons, and it has been 
found convenient and accurate. , 


§1. INTRODUCTION | 


photoelectrons” and in a letter to Nature{, we have briefly referred to our 

method of examining stereoscopic photographs. The method was develo 
originally as the result of a suggestion by Professor C. T. R. Wilson, F.R.S., andi 
suitable for the examination of Wilson cloud tracks. 

Calculations from direct measurements of the photographs have been found 
give satisfactory results in the case of «-rays, where the tracks are approximat 
straight, as is shown by the work of Blackett. In the case of relatively slow f-ra’ 
as, for example, those produced by X-rays, it is well known that the trac 
have generally a zig-zag appearance, and calculations from such photographs would 
be both difficult and tedious. Moreover, in photographing f-ray tracks many 
observers use two cameras with their axes inclined at a small angle, instead of at 
right angles, and the difficulties of calculation in such cases are enhanced. Lough. 
ridge§ has described a special apparatus, the ‘‘stereocomparator,” designed fi 
examining such stereoscopic photographs. 

During the last few years we have taken a large number of photographs of 8-ray 
tracks. The method we have developed for examining them, which will be described _ 
in this paper, has proved both accurate and convenient. 


I: papers on “‘the ranges of secondary §-rays”* and “‘the spatial distribution 


* J. M. Nuttall and E. J. Williams, Phil. Mag. 2, 1110 (1926). 

+ E. J. Williams, J. M. Nuttall and H. S. Barlow, Proc. R. S. A, 121, 612 (1928). 
t J. M. Nuttall and E. J. Williams, Nature, 123, 799 (1929). 

§ D.H. Loughridge, Phys. Rev. 30, 488 (1927). 
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§2. THE PRINCIPLE OF THE METHOD 


Fig. 1 is a plan of the two cameras drawn to illustrate the principle of the method. 
iste that we have a pair of stereoscopic photographs P, and P, from which we 


esire to ascertain the relative positions of points on the object photographed or the 
ontour of the original object. The plates, after development, are replaced in the 
cameras which are adjusted so that they are in exactly the same condition as when 
he photographs were taken. The plates are illuminated by the lamps M,, M, and 
re viewed from the right through the camera lenses L,, L,. Suppose 1,, J, are the 
images on the plates of a common point T on the object photographed, and let 


Fig. I. 


O,, O, be the centres of the lenses L,, L,. Since the relative positions of the plates 
P,, P, and the lenses L,, L, are the same as when the photographs were taken, then 
the point 7”, which lies on both [,O, and J,O,, has the same position in relation 
to the camera as the original point T of which J,, J, are the photographs. Hence 
by ascertaining the point of intersection of [,O; and LO, we can find the position 
of the original point T in relation to the camera. Similarly the positions relative to 
the camera of any number of points 7), T;, ... on the original object may be obtained, 
and hence their positions relative to one another. This is the essential principle of 
the method. If the images on the plates are diffuse through insufficient depth of 
focus, the method.can still be used and will give accurate results provided the 


centre of the diffuse image be chosen as if 


§3. ACTUAL PROCEDURE IN PRACTICE 
An arrangement of pin points A, B, C, ..., mounted on arms with universal 
movements*, is placed in front of the camera (see Fig. 2) roughly in the position of 
the original object photographed. In order to locate and record the position 
relative to the camera of a point T on that object (the point 7 might, for instance, 


* Wireless dealers supply these for use with crystal detectors. 
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be the origin of a Wilson cloud track), one of the pin points is moved so that it 
coincides with the intersection of ,O, and /,0,, where 4, J, are the photographic 
images of T. The fixing of a point in this way may be done by alternate ac 
of the position of the pin on /,O, and /,O, until it lies on both. In practice, how 7 
the fixing of a pin can be done in much less time than this method of successive 


Pmis tA) By tigen. 


Y 


Fig. 2. Apparatus used in fixing the pins. 


adjustment would require. Actually, by looking into one lens with one eye EF, and 
into the other lens with the other eye E,, an observer can readily see the photographs 
in stereoscopic relief, and a pin point can be brought into coincidence with the real 
image at the intersection of [,O, and 1,0, as easily as it could be brought into contact 
with a point on a real object. The relative positions of a number of points on the 
original object may be recorded by the use of several pins, each pin recording the 


position of one point. In this way the contour of the object is mapped (for instance, 
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he contour of a Wilson cloud track), and when the pins have been “fixed” the 
ssential process of translating the stereoscopic photographs of the original object 
into a spatial picture of it is completed. 
Any information required about the contour of the original object may now be 
sbtained from measurements on the pin points. Such measurements can most 
easily be made by projection of the images of the pins on to a screen by means of 
an optical lantern. In this way any desired magnification can be produced, and 
measurements of distances and angles which would be difficult to make on the pin 
points themselves become relatively easy. For this purpose it is convenient to have 
the pins mounted on a rotating table (a graduated spectrometer table serves ad- 
mirably) so that the line joining two pin points may be brought into any desired 
plane relative to the direction of projection by rotation of the table*. ; 
Owing to slight variations in the sizes of the photographic plates the fit of a 
plate in the holder of the camera may be slightly loose, and even though the adjust- 
ment of the camera may be the same as when the photographs were taken the plates 
are not necessarily in exactly the same position as at that time. This defect can easily 
be corrected if a reference point R, fixed at a known distance from the camera, is 
‘photographed at the same time as the object investigated. Then when the plates are 
being examined a pin point can be placed in the known position of R and the plates 
moved until the position of the point fixed from the photographs coincides with the 
known position. The plates are finally clamped in this position. This preliminary 
adjustment is not only necessary to make the reproduction accurate, but is essential 
if the photographs are to be seen stereoscopically at all, since stereoscopic vision is 
impossible if there is any relative displacement of the plates perpendicular to the 
plane containing the axes of the lenses. This adjustment has always to be made 
before proceeding with the measurements of the photographs as described previously. 


§4. ADVANTAGES OF THE METHOD 


The method which we have described in this paper has been used for the 
measurement of Wilson cloud tracks. The advantages of the method may be briefly 
summarised as follows: (1) Very little additional apparatus is required apart from 

“the camera used in taking the photographs. (2) It is not necessary to know the 
stereoscopic angle, the magnification, etc., and the method involves practically no 
numerical or geometrical computation. Moreover, the effect of any glass apparatus 
between the camera and object photographed can be eliminated without calculation 
if such apparatus be simply replaced in the same position relative to the camera 
when the pins are being fixed. (3) Uncertainties due to slight variations in the 
stereoscopic angle, arising from the loose fit of the photographic plates in the holders 


* Alternatively the contour of the object could be mapped by the use of a single pin point 
’ mounted so that it could be moved by graduated screws in any one of three perpendicular directions 
a travelling microscope could be adapted to give the necessary movements. One could thus obtain 
the three spatial co-ordinates of any point by fixing this pin point on the point of the object to be 
located and reading the co-ordinates of the pin from the graduated screws. This method might be 
specially suitable for accurate measurements on a-particle tracks. 
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or looseness of the holders in the cameras, can be eliminated. This is most important 
when small magnification is used. (4) In view of the fact that the images can be see 
in stereoscopic relief, the process of mapping out the contour of an object takes very 
little time. The zig-zag appearance of the Wilson cloud track of a slow B-ray is well” 
known, and the measurement of its total length in space from stereoscopic photo-_ 
graphs would seem a very formidable task. Yet in our work on the ranges of 
secondary f-rays* the range of a B-ray was measured in about 10 minutes. (5) The 
accuracy of the method is limited only by the definition of the photographs. In~ 
order to test the accuracy of the method we photographed a piece of fine copper 
wire about 3 cm. long which was bent into a zig-zag shape somewhat re- 
sembling a Wilson cloud track of a slow f-ray. The cameras had a stereoscopic 
angle of about 15°, the focal length of the lenses was 15 cm., the aperture was f/5-8, 
and unit magnification was used. The bent wire was placed so that it pointed on the 
whole away from the camera, in order to test as far as possible the measurements 
depending on the stereoscopic effect. As the depth of focus was not great the 
photographs of certain parts of the wire were diffuse, so that the test was a stringent 
one in every respect. The developed plates were replaced in the cameras and the } 
contour of the wire was mapped out with pin points from the stereoscopic photo-— 
graphs, as described above, and the total length of the wire was determined by 
measurement on the pin points. The original piece of wire was then straightened 
out and its length measured directly with a travelling microscope. The length 
measured indirectly from the photographs agreed within about 1 per cent. with the ; 
length measured directly with the microscope. a 
§5. APPLICATIONS OF THE METHOD } 
Photography with two cameras at right angles is a special case of stereoscopic _ 
photography and the method is therefore still applicable in principle. But in the case 
of photographs taken at right angles the fixing of a pin point cannot be facilitated _ 
by stereoscopic vision, and this part of the process will therefore take a longer time. 
We have used the method outlined in this paper for the measurement of the 
ranges of secondary -rays, the initial directions of emission of photoelectrons, and) 
the nuclear deflections of B-rays. In all cases we have found the method convenient _ 
in practice and accurate. : s ofl 
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Finally we wish to thank Prof. C. 'T. R. Wilson, F.R.S., who originally suggested _ 
we should try the method, and Prof. W. L. Bragg, F.R.S., for his help and encourage- _ 
ment in the various researches involving its application. 7 
NOTE ADDED JANUARY 22, 1230 ’ 
Further checks on the accuracy of the method have been made by one of us | 
in collaboration with Mr F. R. Terroux at the Cavendish Laboratory, Cambridge, — 
and an account of the results was given in a recent paper by these authors+. 


MA YD ah ¥ 
t E. J. Williams and F, R. Terroux, “The Passage of fast 8-particles through Gases,” Proc. — 
R. S. 126, 294, 295 (1930). a 
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DISCUSSION 


MrT. Smiru said that it was important that the lenses employed in the method 
hould be free from imperfections, particularly coma, and that the centres of the 
upils of the eyes of the observer should be placed and kept in coincidence with 
he projected centres of the lens apertures. Parallax methods of focussing, though 
Imost universally recommended in text books, are carefully avoided in optical 
aboratories. A somewhat similar method has been used in air surveying*. 


Dr L. C. Martin. It is quite true, as remarked by Mr Smith, that problems 
f this kind have been very thoroughly and successfully solved in connexion with 
ir surveying, particularly in the “stereoplanigraph,” an account of which has been 
ublished in the Zeitschrift fiir Instrumentenkunde. In this apparatus two per- 
pective projections represented by photographs are mounted in “‘cameras” in 
ositions related to the camera positions in the original exposures. ‘The instrument 
hen allows of the presentation of a stereoscopic view to the user and with the aid 
f a ‘‘ wandering mark” a contour map of the country can then be semi-automatically 
rawn. Parallax observations are sometimes resorted to by persons who find a 
ifficulty in holding the stereoscopic view; they will open one eye and then the 
ther to see if there is any apparent displacement of the mark with respect to 
he image, but this is hardly likely to be a very accurate procedure. 


AurHors’ REPLY. Imperfections in the lenses are no more detrimental to our 
ethod of examining the photographs than they are to the production of accurate 
hotographs, and the error arising from this source is of the same order of magni- 
ude whichever method is employed for the examination of the photographs. 
n fixing a pin the observer uses the central rays through the lenses by adjusting 
he position of his eye so that the image is in the middle of the field of view. Errors 
introduced by the use of rays through the periphery of the lenses depend on the ratio 
f their apertures to their distance apart in the stereoscopic camera. As the diffuse- 
ess of the photographic images of point objects is due to the rays coming from 
ifferent parts of the lenses an error of the same magnitude will be incurred if, 
n determining the original position in space of the points photographed, measure- 
ents are not made on the centres of their diffused images on the plates. The 
xtent to which the conditions for accurate measurement by the present method 
an be adhered to is shown by the very satisfactory results quoted and referred to 
n the paper. The method used is in certain respects the same in principle as those 
sed in air survey and it may be regarded as an extension or perhaps a simplification 
f such methods. In it a full-scale reproduction of the shape of the original object 
is made, and no optical apparatus in addition to the camera itself is used. 


* J. C. D’Almeida, “ Nouvel appareil stéréoscopique,” Comptes Rendus de l’ Acad. Frang. 47, 
r (1858). T. Scheimpflug, “The use of the skioptikon for the production of maps and plans from 
hotographs.”’ Photographische Korrespondenz, 35, 114 (Vienna, 1898). O. von Gruber, “The Carl 
eiss stereoplanigraph,” Zert. fiir Instrumentenkunde, 43, 1 (1923). J. C. Kapteyn used the principle 
in determining the parallax of stars ; several papers on this work were published about 1890. 
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A DETERMINATION OF THE COEFFICIENT 
OF DIFFUSION OF WATER VAPOUR 


By W. E. SUMMERHAYS, B.Sc. . 


Communicated by Prof. S. W. }. Smith, December 25, 1929. Read and 
discussed February 14, 1930. 


ABSTRACT. Water vapour was allowed to diffuse down a vertical tube about 2 in. im 
diameter, a steady concentration-gradient being maintained from approximate saturation 
at the top to approximate dryness at the bottom. The gradient was measured with two 
katharometers at a definite distance apart and the mass of vapour passing in an observed 
time was weighed. The result given is 0-282 cm.2/sec. at mean temperature 16-1° C. 


§1. INTRODUCTION 


have been made for thirty years, and I have been unable to find some of the 
original papers on the subject. In those determinations of which full — 
particulars could be found the methods adopted seemed to be open to criticism, and 
it was thought that the accuracy attainable with a katharometer* used as a hygro- 
meter would enable a method free from serious error to be found. 
Two distinct methods were attempted, (1) a direct method, (2) an indirect 
method. Whilst the direct method gave consistent results with different experi-_ 
mental conditions, the indirect method had to be abandoned. | 


N determination of the coefficient of diffusion of water vapour appears to 


§2. THE DIRECT METHOD 


A steady water vapour pressure-gradient was established along a vertical tube — 
and the mass of water vapour that had passed in a known time was found. The 
gradient was measured by means of two katharometers inserted in the tube a known | 
distance apart. It was considered essential that the diffusion should take place 
downwards as water vapour is lighter than air, so the top of the tube was covered — 
with a damp surface and the bottom was kept dry with phosphorus pentoxide. 

The apparatus consisted of a brass plate A, Fig. 1, 20 cm. in diameter and 
1°5 cm. thick. This was placed on the base of a retort stand R, the shortened rod 
of which passed through a hole in the plate 4: the rod thus constituted a vertical 
axis about which the plate could be rotated by means of a ‘‘tommy-bar” inserted 
in holes drilled in its edge. Two wells 6 cm. in diameter and 1 em. deep were cut 
in the plate. A similar plate B, 1 cm. thick rested on A, and was prevented from 
rotating. In plate B one hole 4-85 cm. in diameter and one 6 cm. in diameter were 
cut. The larger hole was normally covered with a flanged brass cap C. Over the 


* Cf. Daynes and Shakespear, Proc. R. S. A, 97, 237 (1920). 
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ee hole was placed a copper disc FE, 10 cm. in diameter and 0-6 cm. thick. This 
. was perforated over an area 5 cm. in diameter with holes 2 mm. in diameter 
aving centres 4 mm, apart. On the copper disc was placed the long tube, 4°85 cm. 


7 SORE KON 


Ieee 


elevation. Right: plan, with tube, perforated disc and cap removed. 


in internal diameter, in which the gradient was established. This tube consisted of 
three pieces X, Y, Z, connected together by means of two brass collars F,, F, each 
4:6 cm. long into which the pieces were screwed. The top and bottom pieces of the 


tube were each 10 cm, long and were flanged at their free ends. The middle piece 
15-2 


Fig. 1. Left: sectional 


=" 
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was 125 cm. long for half the experiments and 17-5 cm. long for the other half. 
Through a hole in the side of each collar a katharometer K,, K, was screwed and 


220 


its position adjusted until its base was flush with the inside of the collar and tube. — 


A copper disc 15 cm. in diameter and 0-3 cm. thick (to act as a cooling fin) was 
threaded on each katharometer and clamped to the collar. On the top of the tube 
was placed another perforated copper disc E,, but this disc had wound round its 
edge eight turns of manganin wire (No. 29 s.w.g.) kept in Position by means of 
shellac. In the top and bottom pieces of the tube a short side tube was soldered 
2 cm. from the flange. These tubes were normally closed with brass plugs through 
each of which passed one junction of a copper-constantan couple ff. During an 
experiment an unglazed, porous soup plate P containing distilled water was placed 
on the upper perforated copper disc and was covered with a glass plate G. The tube 
and collars were coated internally with paraffin wax and all surfaces in contact 
(except the top of the porous plate) were ground and greased with a mixture of 
beeswax and vaseline. 

Preliminary experiments were made, without the use of katharometers, with 
a short tube of large diameter, and it was found that the cooling due to evaporation 
from the lower surface of the porous plate and the heating of the phosphorus 
pentoxide due to absorption caused the bottom of the cylinder to have a tempera- 
ture higher than that at the top. A small electric heater was put in the water at the 
top but it was found impossible to get the temperatures equal, the temperature at 
the bottom being always the higher. A consideration of these facts led to the adop- 
tion of the perforated copper discs. The katharometers were designed to take 
0:35 amp. and a continuous current of this value raised the temperature of the 


katharometer blocks in air about 2° C. above room temperature; so the katharo- — 


meters were placed on opposite sides of the tube and the copper fins were fitted. 
When the katharometers were in place in the tube and a water vapour pressure- 
gradient had been established, it was found that with a continuous current the 
thermo-couple showed that the temperature at the top of the tube was nearly 1° C. 
above that at the bottom. As this temperature difference was in the right direction 
it was hoped that it would not give rise to convection, and a series of experiments 
was made with’a middle piece of the tube 17-5 cm. long. Results from these 


experiments were moderately consistent, giving a value of 0-32 cm.*/sec. for the’ 


diffusion coefficient. When the middle piece of the tube was made 12-5 cm. long 
and the experiments were repeated, the results were no longer consistent. The 
absence of convection being an essential condition for an accurate determination, 
it was decided to place in the cylinder along its whole length two strips of metal 
bent and soldered together to give a cross-shaped cross-section so as to divide the 
tube into four compartments, ‘This metal was well greased and put in position and 
the value obtained for the coefficient was 0-28 cm.%/sec., showing that convection 
had been present before the introduction of the metal strips. It was therefore 
decided to pass the current through each katharometer for 30 seconds when a 
reading was required, and to take readings at intervals of not less than half an hour. 
When this method had been adopted it was found that the temperature at the top 
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of the tube was 0-3° C. lower than that at the bottom. By the passage of a current 
of 0-12 amp. through the manganin wire wound round the top copper disc EF the 
temperatures at the top and bottom of the tube were sensibly equal (less than 
o-r C. difference). 

To prevent the initial drop in e.m.f. that occurs when a current from accumu- 
lators is switched on, it was arranged that except during the 30 seconds when a 
reading was required the current should pass through a resistance approximately 
equal to that of either katharometer. The room chosen for the experiments was one 
in which temperature changes were slow, and it was found possible, by taking 
advantage of suitable meteorological conditions and by manipulating a very small 
circular window near the ceiling, to keep the room temperature constant within 
0-2° C. during the whole experiment. When it was found impossible to keep the 
temperature within this limit the results of the experiments were rejected. 

Calibration of katharometers. The katharometers were calibrated in situ by 
the pushing of air of known vapour pressure through the tube, the porous plate 
P being replaced with a ground glass plate. An aspirator of large capacity was 
filled with air from outside the building and was then made to push this air through 
a ‘“saturator”* to the lower side-tube of the apparatus. A short rubber tube con- 
nected the upper side-tube to a ‘‘ bubbler,” so that the rate of flow could be observed 
and regulated. The saturator consisted of a copper plate 16 cm. square, on the 
under side of which a long spiral copper strip 1-8 cm. wide was soldered edge-on. 
This plate rested on the top of a large copper beaker containing water which reached 
half way up the spiral. Air was led into the spiral on the outside, and passed out 
through a tube in the centre. ‘This central tube passed through a solid copper 
cylinder soldered to the top of the plate. The temperature of the air as it left the 
saturator was observed on a thermometer inserted in a hole drilled in the copper 
cylinder. 

The rest of the top of the plate was covered with a thick piece of wood. An 
ordinary pail was lagged with felt and the copper beaker was placed on three corks 
on the bottom, the pail being then filled with water or ice. It was found that the 
temperature as shown by the thermometer could be kept constant within o-1° C. 
by means of small pieces of ice dropped into the pail. A temperature about 2° C. 
below that of the tap water was the highest chosen. Readings of the katharometers 
were not considered satisfactory until they had remained constant for at least an 
hour and a half. 

A “dry zero” for each katharometer was obtained by the passing of air from the 
aspirator through phosphorus pentoxide tubes and a long sulphuric acid tube. The 
dry zero readings are not of great value as it is found that with prolonged drying the 
katharometers register a gradually increasing dryness. The dry zero points shown 
on the calibration curves are those obtained after the passage of dry air through the 
"apparatus for three hours, but all values of vapour pressure taken from the curves 
lie between points obtained with the saturator, or close to one of them. The effect 


* This “saturator” was designed by Dr Shakespear and had been used by him in calibrating 


katharometers for some experiments on humidity. 
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on the katharometer readings of variations of atmospheric pressure ids ici : 
but were found to be very small and were sensibly equal in i two katharome 
used; no correction was therefore necessary except 1n the c ses P 
To start the experiment the porous plate P containing water an — 
the glass plate G was placed on the tube and a Petri dish containing P,O; was 


Katharometer deflections 


Vapour pressure in mm. of mercury 


Fig. 2. Katharometer calibration curves for water vapour. 


in the well under the tube. Another dish of P,O, was placed in the other well anc 
covered with the cap C. When the katharometer readings showed that the gradient 
was steady (after about two hours) another Petri dish of P,O; was prepared an 

covered with a greased brass plate of the same diameter, and was carefully weighed 
The plate adhered to the dish but could be detached by pushing it with a fine wi 
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hrough a small hole drilled close to its edge. The dish of P,O, in the well under the 
ap C was removed and the weighed dish substituted and its cover detached. ‘The 
ap C was replaced and the lower plate A rotated until the weighed dish of P,O; 
ame under the tube. 
Readings of the katharometers were taken at intervals of about half an hour 
: nd the experiment was allowed to run for from about 3 to 5} hours. The dish was 
then removed and reweighed and the diffusion coefficient calculated. 


Table of Results* 


Distance between 
katharometers (cm.) Temas Som 2 mee Sean ZO. Cen 20°0" "| 20°0 | 2000 
Tem t i- 
Bice me on I5°4 15°5 161 16°3 577) 16:0 16:8 Tejas 
Bottom katharometer : : ‘ : 
Heflection aye Bye 3) 5°4 4°3 4°4 48 4°9 
Vapour pressure at 
bottom katharometer | 4°7 4°7 4°9 50 41 4:2 4°5 46 
| in mm. of mercury 
| Top katharometer de- : : ; : ; 4 : : 
ae oreo casters Ss \17 AN 1X73 eT 7-8. sh 18:7 #1. 19-0 
Vapour pressure at top 
katharometer in mm. 9°5 96 10°05 | Io'l TO'O5 10°35 4) 10:85-, 1270 
of mercury 
Difference in vapour ; ; : ; ! : : : 
eecuces 4°8 4°9 See ae! 505 Wee 5 0) 60-35 | 04 


Mass of water vapour | 27-3 EO nao 200 | 25-6 16:8 | 325 | 25°4 
in m.gm. 

Time of experiment 6 
a 221 260 193 194 230 170 330 252 


Coefficient of diffusion 
in cm.?/sec. 


0274 | 07281 | 0°282 0:278 | 0:284 | 0:289 | 0:280 0:285 


Mean value of coefficient for 15:2 cm. distance: 0°279. 
Mean value of coefficient for 20'0 cm. distance: 0'2845. 
Mean value of coefficient for all experiments: 0'282. 
Mean temperature for all experiments: nope (Ee 


The experiments were made at room temperature and the results are for 
temperatures between 15-4° C. and 17:3° C. Four results were obtained with a 
distance of 15-2 cm. between the katharometers and four with a distance of 20-0 cm. 
They depend on the difference in vapour pressures obtained from the calibration 
curves, and the error possible here is thought to be not greater than 2 per cent. 
An attempt was made to weigh the water evaporated as well as that absorbed. 
For this purpose a small unglazed porcelain disc was cemented near the bottom 
of a short brass cylinder 5:5 cm. in diameter. Water was placed in the cylinder 
which was then covered with a greased ground glass plate. The cylinder was placed 
on another greased plate whilst being weighed and was then placed on the tube after 
the porous plate had been removed. It was found that the loss in weight was some- 


* The mass of 1 litre of water vapour at n.t.p. was taken as o'806 gm. 
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times greater and sometimes less than the gain, though the difference was always 

very small and could be accounted for by variations in the amount of grease adhering : 

to the lower surface of the cylinder. The possible disturbances set up in the tube 

by the changing of the water vessels led to the abandonment of further work in this 
direction. 

For steady diffusion the coefficient K is given by the equation ; 

mass 
~ area x time x concentration gradient 


For the particular application considered here this equation becomes 
mass x length x 760 x 1000 x (273 + #) 
~ 7 x radius? x time x 0-806 x 273 x vapour pressure difference 


where ¢ is the centigrade temperature. 


§3. THE INDIRECT METHOD 


The close analogy that exists between thermal conductivity and diffusion — 
suggested the possibility of finding the diffusion coefficient by the use of a mathe- 
matical equation previously worked out for conductivity. 

It has been shown that if a semi-infinite solid be at a given temperature and its 
bounding plane be then kept at zero temperature, the time ¢ taken for the tempera- 
ture at a depth x to reach its half value is given by the equation 

x/2,/kt = 0:477* 
where k is the coefficient of thermal conductivity. An attempt was made to use this 
equation to find the diffusion coefficient. 

With slight alterations the same apparatus as that used in the preceding method 
was adopted. The porous plate and perforated copper discs were removed and a 
tube 52 cm. long flanged at both ends was placed on the top of the original tube and 
was closed with a glass plate. The bottom piece X of the tube was replaced by a 
similar piece 20 cm. long. A Petri dish containing water was placed under the tube 
and left until the bottom katharometer A, gave a reading rather less than that 
corresponding to saturation for the temperature of the room. The dish of water 
was then removed and the apparatus left until both katharometers gave readings 
that did not alter for one hour, thus indicating uniformity of vapour pressure 
throughout the tube. A dish of P,O,; was then placed under the tube and the time 
taken for the bottom katharometer to reach its half value was observed. The experi- 
ment was first done before the inside of the tube was coated with wax. The values 
obtained in this case varied between 0-18 and 0-22. After the coating of the tube 
with paraffin wax the values were much more consistent with a mean of about 0°23. 
It appeared that the method was invalidated by the adsorption of water vapour on 


the inside surface of the tube and also by the unavoidable lag in katharometer 
indications. It was therefore abandoned. 


P x Gieenee Introduction to the Mathematical Theory of the Conduction of Heat in Solids, 
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DISCUSSION 


Mr J. H. Awsery suggested that a third katharometer, intermediate between 
he two actually used, might have been introduced for the purpose of ascertaining 
hether the concentration gradient was constant throughout the length of the tube. 


Mr R. S. WurePLe congratulated the author on the thoroughness of his work 
and suggested that temperature changes, which are the chief source of difficulty 
in this kind of investigation, might have been overcome by immersion of the whole 
apparatus in a bath of constant temperature—a method adopted by the Bureau 
of Standards. 


Mr RoLto APPLEYARD. With paraffin wax directly exposed to a damp atmo- 
sphere—as was found in the tropics with early forms of electrical condensers—there 
is a falling off in dielectric resistance, indicating the presence of moisture. I suggest 
that a material better suited to the purposes of the author could be found as a 
lining to the tube. 


Dr D. Owen. I feel some surprise at the observations of the last speaker. 
Experiments I have made have shown the small screening effect of an electric 
field due to the interposition of thin dried and waxed sheets of paper, and confirm 
the usual view as to the extraordinarily high insulating properties of parafhin wax, 
even in a moist atmosphere. 


Avutuor’s REPLY. The nature of the problem and the design of the apparatus 
made the adoption of a bath of constant temperature a matter of very great difficulty. 
If it had not been found possible to keep the room temperature within 0-2° C. 
some form of jacketing would have been attempted. In the direct method with a 
steady water vapour pressure-gradient adsorption of water vapour on the walls of 
the tube should not affect the results, as a state of equilibrium would be established. 
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A NEW TYPE OF DEWAR FLASK, FOR USE 
AS A CALORIMETER 


By W. L. WATTON, B.Sc., A.R.C.S., D.LC. 
Imperial College of Science and Technology 


Communicated by Dr H. R. Lang, Fanuary 11, 1930. 
Read and discussed, February 28, 1930. 


ABSTRACT. The water equivalent of a new type of Dewar flask, of which the inside 
was of copper, has been measured at laboratory temperatures, and found to be more 
constant that the usual type. The precautions necessary in using such a calorimeter are” 
fully discussed. 


§1. INTRODUCTION 


to be expected from calorimetric experiments which were to be made with 

a Dewar flask and an electric heating coil. As previous observers had found 
considerable variation of the water equivalent of such an arrangement*, experiments 
were made upon a vacuum flask in which the inside glass vessel was replaced by a { 
similar vessel of copper. 


T HIS series of experiments was carried out to determine the order of accuracy 


§2. PRELIMINARY EXPERIMENTS 


‘ 

Some preliminary experiments were made by Dr H. R. Lang in June 1928 with” 
a vessel similar to that used in the present experiments. A stirrer which — 
had a heating coil wound on it was caused to oscillate in a vertical direction 
A charcoal pocket was introduced to keep the vacuum high, and it was found 
advisable to use the same mass of liquid in successive experiments; a correction | 
was made for the heat due to stirring. It appeared that the water equivalent was” 
related to the amount of electrical power supplied. The writer set out to investigate 
this relation more completely, with the ultimate intention indicated above. The 
temperature/time results obtained in these early experiments have been plotted 
(§ 7) and found unsatisfactory. In the same way, an attempt to formulate hel 
thermal effect of stirring was not very successful, the values of the relevant constant 4 
(§ 8) differing by 50 per cent. The apparent relationship between the water equi- 
valent and the power input is due to (a) insufficient stirring and inaccuracy in 


the evaluation of the stirring correction, and (5) the loss of heat due to evaporation 
on the wetted surface of the oscillating stirrer, ; 


j 
* Experiments made by Dr H. R. Lang and Mr W. Goldsmith on an ordinary glass Dewar _ 


flask of about 800 ¢.c. capacity showed a sporadic variation of water equivalent of more than 
20 per cent. 
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§3. THEORY AND OUTLINE OF METHOD 


A heating coil is immersed in a liquid contained in the flask, along with a 
thermometer and a stirrer. The heat communicated to the liquid and calorimeter 
is evenly distributed by the stirrer, and the thermometer registers the increase of 
temperature of the calorimeter and its contents. 

Elementary theory shows that if 


K be the mechanical work done by stirring per sec. in joules, 

M the mass of liquid in gm., © 

s its specific heat in cal. per gm. deg. C., 

w the water equivalent of the calorimeter in cal. per deg. Gs 

(dT/dt),, the rate of rise of temperature due to stirring only, 

(dT /dt), the rate of rise of temperature due to stirring and electric heating, 

E the potential difference across the heating coil in volts, 

C the current through the heating coil in amperes, 

@ the observed temperature in degrees C. of the water bath, 

T the observed temperature in degrees C. of the calorimeter, 

i the time in seconds, 

¥ the mechanical equivalent of heat in joules per calorie, 

p, the cooling in deg. C. per sec. per deg. C. of temperature gradient during a 
stirring experiment, 

p, the cooling during an experiment with current flowing, 


then 
K/F = (Ms + w) {(dT[dt)m + pr(T— states (1). 
This can be written 
K/f = (Ms + w)/D, where D = {(dT dt) m + pi (T — )}* .--- ay 
When electrical energy is supplied to the coil, this becomes 
(K + EC)/F = (Ms + w) {(dT/dt), + px(T —)} sever (3); 
which can also be expressed thus: 
(K + EC)/F = (Ms + w)/F, where F = {(dT/dt), + p2(T — 9} ---(4). 
Therefore 
EC/¥ =(Ms+w)(3/F—1/D) tt (Ss 


Now if, during an experiment, (T — @) passes from a negative value, through zero, 
to a positive value, the value of F (or D in the case of stirring) can be found when 
7 —6. Hence errors due to heat loss should vanish. Even were this not so, if 
” p, and p, are very small and almost, if not quite, equal, the error due to the neglect 
of the difference between the effects in the two experiments is negligible. Therefore 


from 
ee) ECi} = (Ms + w) {(dT/dt), — (CARCASS tents (6). 
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§4. MEASUREMENTS 


(a) Electrical measurements. ‘The heating current was generated by 24 Chloride _ 
batteries, each of 40 ampere-hours’ capacity, connected in two banks of twelve cells, 
and these two banks were then connected in parallel. In series with the battery was 
a variable manganin resistance, a standard resistance, an ammeter and the heating 
coil. A manganin dummy coil of the same resistance as the heating coil could be 
substituted for it by a two-way mercury contact. The current was passed through — 
this dummy coil immediately before each experiment to enable it to become 
steady. The standard coil was of 10 s.w.g. manganin, immersed in @ water-cooled 
oil bath. The resistance of the coil was o-9961 ohm at 20° C. and its temperature 
coefficient negligibly small. 

The potential differences across the standard resistance and the heating coil 
were measured with a calibrated 4-dial potentiometer in terms of standard Weston- 
cadmium cells. A slab of bakelite, in which were mercury cups, enabled the 
various connexions to be made rapidly, a connector being placed in three different 
positions for the three quantities required. 

(b) Thermometry. The temperatures of the water-bath surrounding the Dewar 
vessel, and of the vessel and its contents, were taken with a platinum thermometer 
in conjunction with a Callendar-Griffiths bridge, which was carefully calibrated in 
terms of the largest coil. The thermometer had a fundamental interval of 4-775 ohms. 
The temperature of the laboratory was taken with a mercury thermometer, 
graduated in half-degrees; the ice and steam points were found to be at + 0-7° C. 
and roo-1° C. respectively, when corrected; readings of the mercury thermometer 
have been corrected accordingly. 

(c) Weight and volume. The volume of water necessary to fill the calorimeter up 
to 2cm. from the top was found by trial, and a file mark corresponding to this 
volume was made upon the neck of a flask. The flask was weighed with calibrated 
brass weights when filled to the mark and also after the contents had been emptied 
out. The volume of water used was 559-2 cc. at 17° C. It was found that, at the — 
same temperature, the same mass of water could be measured off by volume to 
0-2 gm. In most cases, the water was measured off in this way. 

(d) Time. The rate of rise was measured as described below (§ 6). The standard 
clock marking seconds on the chronograph was rated against the broadcast 
Greenwich time signals. It was found to have a variable rate, from zero to a loss 
of 20 seconds a day, This error, though large for such a clock, would still not be 
sufficient to affect the results. Other times were taken with a ship’s chronometer. 


§5. THE CALORIMETER 


The Dewar flask was about 50cm. long and 5 cm. in diameter, and had a 
capacity of about 800 cc. The usual inside of glass was replaced by a copper 
vessel of similar shape and size, the joint being made air-tight with sealing wax. 
The glass was made by the National Glass Industry of London. The joint (R, Fig. 1) 
was about 1 cm. from the top of the flask and an equal distance above the level V__ 
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f the water inside, when the apparatus was completely assembled. ‘The exhausting 
ube at the bottom A was connected to a Hyvac pump and to a discharge tube and 
harcoal pocket. The space was first evacuated by the Hyvac pump, and more 


ompletely by the charcoal pocket in liquid air, In less than 15 minutes a discharge 
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Fig. 1. Dewar flask with copper inner vessel, arranged as calorimeter. 


coil capable of giving a I cm. spark in air was unable to produce a discharge across 
the discharge tube. This tube was about 16 cm. long and was of the usual type. 
The modified Dewar flask was almost immersed in a lagged water bath, about 
2 cm. being above the surface. ‘The water was made to circulate by means of a pump 
driven by an electric motor. The temperature of the jacket did not change by more 
than o-5° C. per hour on any occasion, and the variation was usually less than this. 


230 W. ‘ee Watton 


The narrow and comparatively long shape of the calorimeter placed considerabl 
limitations upon the form of the heating and stirring arrangements. Particul 
was it difficult to find a piece of apparatus of the required compass which wo 
ensure that the water was thoroughly mixed and the heat distributed unifo nly 
throughout the whole mass. Several months were spent before an efficient stirring: 
system was designed, and even with this form of stirrer it was found necessary 
to stir at the rate of 35 revolutions a second if the results were to be consistent. 
A copper cylinder B (Fig. 1) was closed at the top and bottom with brass plates" 
G and C. A bearing in each held the propeller shaft. The bottom disk C was 
pierced with several holes so that it was a mere grid of metal. Equally spaced 
around its circumference were four supports D, portions of a cylindrical shell, 
inside which the multibladed propeller E revolved. The top disk G was coarsely 
serrated around the circumference which was a little larger than that of the copper 
cylinder. In the top of this cylinder, and about 1 cm. beneath the top baffle 
plate G, were bored a large number of holes F. The diameter of both baffle plates 
was a little larger than that of the cylinder, and was such that they loosely fitted 
the inside of the calorimeter. The propeller E was revolved in such a direction as 
to suck water down the inside of the copper cylinder and pump it up the outside. 
In each circulation the water had to pass through two baffle plates and suffer 
churning up in the propeller chamber. Moreover, should the water inside the 
copper cylinder have a slightly different temperature from that surrounding it, the 
high conductivity of the copper of the cylinder would promote an equalisation of 
temperature. The thermometer 7 was placed inside the copper cylinder through 
holes made to fit it tightly. : : 

Around the outside of the copper cylinder and equally distributed were six 
long, thin ribs H of bakelite; the heating coil ¥ was wound non-inductively in 
notches in these. In order that the heating effect might be distributed over a 
volume, a distance of about 2 cm. separated adjacent turns. For ease in tess 
later, the resistance was adjusted to be nearly 1 ohm, and about 1} metres of 
20 S.w.g. manganin wire were employed. In order to minimise errors due to a 
wrong position of the potential leads, the ends of this coil were slipped into a hollow _ 
copper cylinder K formed by boring a hole to fit the wire in pieces of 12 s.w.g. 
wire about 1 cm. long. The external current leads L were of 16 s.w.g. copper wire, 
the ends of which were filed down until they fitted in the open ends of the con-_ 
necting pieces. ‘The wires were soft-soldered in this position, the solder filling the 
inside spaces and making a solid joint. The potential leads M were of 28 S.W.g. 
silk-covered copper wire. ‘The covering was removed from that portion of the wire 
which was inside the calorimeter, and the wire was soldered to the end of the 
connecting piece K nearest the heating coil. All the bare wire in the calorimeter 
was then coated with an alcoholic solution of shellac and a current of some 5 amperes 
was passed through it so as to bake the shellac. Three coats were put on in this way. 

The shaft of the stirrer N was in two parts connected by a bakelite sleeve P. 
This served to diminish the conduction loss along the shaft. 

A lid of box-wood Q was pierced with holes to take the four leads, the stirrer 
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haft and the thermometer. The holes for the current leads tightly fitted the wires, 
he rigidity of which kept the lid in place. The lid fitted the top of the calorimeter 
losely, and evaporation during an experiment must have been small. The stirrer 


as driven by a shunt wound motor, a controlling resistance being in series with 
€ armature. 


§6. METHOD OF MAKING AN EXPERIMENT 


The pump circulating the water in the jacket was set going about 20 minutes 
before the experiment proper. The heating current was adjusted to the required 
value while passing through the dummy heating coil, and the measuring flask was 
filled with distilled water and cooled. The charcoal pocket was next immersed in 
liquid air, and by this time the temperature of the water bath was steady and its 
value could be determined, the time shown by the chronometer being noted. The 
temperature of the laboratory was observed immediately afterwards. 

The measuring flask was removed from its cooling bath and the contents, 
which were about 10° below room-temperature, were emptied into the Dewar 
vessel. The stirring mechanism was inserted and the leads were connected to a row 
‘of four conveniently placed terminals. The platinum thermometer was inserted 
after drying, a pulley and collar were screwed on to the stirring shaft, and the 
revolution counter was fitted to the collar. The experiment was not proceeded with 
until the nature of the vacuum had been tested and found to be “black.” 

The stirrer was started and its speed adjusted by trial to the required value. 

At first some trouble was experienced in keeping the speed of stirring reasonably 
constant, until the original motor was exchanged for another of a modern type, 
made by the General Electric Company, of 1/10 h.p. This ran fairly constantly, 
the total variation at the slowest speeds being no greater than 3 per cent. and at 
higher speeds very much less. This variation was probably due in some measure 
to the fluctuating voltage of the mains from which its current was derived. The 
- stirrer shaft revolved about three times as fast as the motor. 
The current was switched over from the dummy heater, and the chronograph 
“was set going. At some definite time, as read from the chronometer, the experimenter 
made a mark on the chronograph drum by tapping the key and allowing it to rise 
in unison with the half-second beats of the chronometer. In this way any particular 
record on the chronograph could be referred to events happening at the same time 
in other determinations. 

The bridge was set to be balanced at a temperature corresponding to 15 cm. 
of bridge wire (about 3° C.) below the jacket temperature, and the chronograph 
key was depressed as the galvanometer spot passed through zero, which was set to 
be at the position of balance. The galvanometer had excellent zero-keeping pro- 
perties and these were carefully preserved by fairly frequent reversal of the current 
- when the deflection was at all large. Any change of zero was carefully noted at 
‘the end of each experiment; but this occurred to so small an extent as to render 
any error due to it very small, Such errors, unless they were progressive, which 
is doubted, were eliminated in the method of reducing the results. In this way, a 
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mark was obtained for each 5 cm. of bridge wire until the calorimeter had reached 
a temperature corresponding to 15 cm. above the jacket temperature (as observe 

inni f the experiment). 
; eee cee of a ithe the p.d. across the standard resistance and 
across the heater were determined alternately by a second observer. The reading 
of the counter was taken periodically. At the end of an experiment, the vacuum > 
was tested again, and the jacket and laboratory temperature were redetermined, the 
time being noted. 

The water in the Dewar flask was sucked out by a filter pump. The flask was 
then dried with a mop, some methylated spirits were poured in and mopped out, 
and the inside was finally dried by the blowing in of air until it appeared dry by 
inspection and there was no smell of methylated spirits. 


§7. REDUCTION OF RESULTS 


After removal of the record from the chronograph the whole minutes, indicated 
by a special mark made by the clock, were numbered from zero. By reference to 
the mark indicating a known chronometer time, the chronometer time corresponding 
to the zero minute mark was known. This was called “‘zero time.’’ Hence if the 
chronometer time corresponding to any other mark was required, it could be 
obtained instantly by adding zero time to the chronograph time of that record. 

In a similar way, the marks corresponding to each 5 cm. of bridge wire were 
given serial numbers commencing with 1, and these numbers were denoted — 
collectively by the symbol G. The chronograph time for each of the marks was 
read off to the nearest tenth of a second, a procedure which was simplified by — 
adjusting the chronograph so that 1 cm. of record corresponded to 1 second of 
time. The time interval between each successive pair of marks was obtained by 
subtraction, and the intervals were denoted by the symbol F. (In stirring experi-— 
ments, the symbol D was used because the range of temperature between two — 
marks was 1/10 that used in the main experiments.) 

The intervals F as ordinates were plotted against the corresponding mean 
values of G as abscissae. ‘The points were found to lie on a straight line (referred 
to as a G line) as required by the approximate expansions of equations (2) or (4); 
for G bears a linear relationship to JT. Sometimes the points were very far from 
being collinear; this was almost invariably found to be the case when the stirrer was _ 
rotating at a speed lower than 35 revolutions a second. In most of the experiments 
described, the stirring was faster than this; other experiments were taken into 
consideration only if the G line was satisfactory, This linear relationship seemed to 
be a crucial test of efficient stirring, and was so regarded in these experiments. 
‘Typical G' lines for (a) good, and (6) fair, stirring conditions are shown in Fig. 2. 

The mean jacket temperature, during the period when records were being 
made, was calculated by simple proportion, and the G value corresponding to this 
temperature was worked out, This value of G when corrected for lag was that at 
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hich the jacket and calorimeter temperatures were the same and was denoted 
y G,. The corresponding value of F (namely /)) was read from the G line. 


But from equation (4) F, = 1/(aT/dt) 
o= , bey 


hence (dT/dt), is known. Stirring experiments giving the value of (dT/dt),, were 
orked out in the same way. The reciprocals of Fy were calculated and the stirring 
orrection was subtracted. The difference ‘was converted to platinum degrees and 
hen to the gas scale, the mean correction over the range of temperatures used in 


160 


355-7 revs. per sec. 
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Fig. 2. (a) Above critical rate of stirring; (b) below critical rate of stirring. 


the experiment being adopted. The electrical power supplied was calculated from 
the potentiometer readings on the assumption that the e.m.f. of the cells was 
r-0r83 volts. The change of e.mf. with temperature over the relevant range was 
$0 small that a correction was not necessary. This fact and the value of the e.m.f. 
were verified experimentally by comparison with N.P.L. standards. The tem- 
perature variation of the heat-capacity of the volume of water used was allowed for, 
the weight being corrected for the buoyancy of the air. The mechanical equivalent 


of heat was taken as 4-180 joules per calorie. 


§8. STIRRING 


Experiments on the correction for stirring were carried out in the same way as 
the main experiments, except that the temperature interval over which the time to 
‘tise was taken was diminished to 1/10 or 1/20 of that in the main experiments. 
Much greater difficulty was experienced in obtaining consistent results for the 
rate of rise with stirring than was experienced when electrical energy was used. 


The rate of rise in cm. of bridge wire per second at the mean jacket temperature 
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was derived from G lines. It was found that this rate 1/D) was connected with the 


rate of stirring by the law A=Dr; 


for, provided that the stirring was fast enough to mix the liquid thoroughly, when 
log D, was plotted against log r the points were on a straight line, from the slope 
of which the value of m could be obtained. It was found that, in general, when the 
stirring was below 35 rev. per sec. (i) a G line could not be drawn accurately, if at 
all, and (ii) the stirring did not obey the law quoted or any other ascertainable law, 
being quite inconsistent. 

Above 35 rev. per sec.*, however, the law was obeyed with fair accuracy, the 
value of m being 2:20 while the mean value of A was 3-03 x 10° when 7 was ex- 
pressed in rev. per sec. and D, in cm. of bridge wire per sec. 

Table 1 shows the accuracy obtained, the correction calculated being compared 
with that obtained experimentally. Owing to the difficulty of obtaining reliable 
values experimentally, it is considered that the calculated correction is more — 
accurate than a single observed value and it has accordingly been applied to the 
experiments. 

Table 1. Stirring corrections 


Rateya.p-Somer 31'0 33°8 34°4 35°2 | 38°5 ) 43°0 
Observed 0:0064 00068 0°0079 00089 0°0099 ) o°or3I 
Calculated 0°0063 0°0076 070079 | 070083 | o-or02 | O°OI29 


The stirring correction amounted to 1/5 of the total heating effect when the 
electrical power was 20 watts ; but diminished in importance as the power increased. 
Accordingly, great care was taken to obtain a sufficiently accurate value of the 
correction. 


§9. SOURCES OF ERROR 


Heat loss. Errors due to this cause were reduced to a minimum by the careful 
maintenance of a high vacuum between the vessels of the Dewar flask, and by 
details in the design of and method of working the apparatus already described. 

Lag of the thermometer, The platinum thermometer had an appreciable lag, 
though encased in the thinnest glass tube available. This lead to two errors in the 
value of £) (or Do) as read from a G line: (i) when the thermometer registered the 
same temperature as the jacket, the temperature of the calorimeter was higher than 
this: hence the condition T’ = @ was not satisfied; (ii) since the lag was assumed 
to be proportional to the rate of rise of temperature, the amount of lag would change 
with F and the G line would be turned through an angle. 

The first effect was found to cause an error of about 1 per cent. on the value 
of the water equivalent, but the error due to the second cause was about 0-1 per cent. 
and was therefore negligible. 

The thermometer was cooled and then placed in the water bath which was at 
room temperature, and was continually stirred. A series of readings was obtained 


: . ; oe 3 
Occasionally satisfactory stirring conditions obtained at lower speeds than this. 
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with the galvanometer and chronograph. ‘These were plotted on a time/temperature 
curve to test their regularity. 


Elementary theory shows that 
aL dica( Oi 1) (A EE Ps <a2 G7) 
where ® is the temperature of the liquid, 
T the temperature of the thermometer at any time ¢, 
da constant (with dimensions of time). 
We have also , 
© = ©, + mt, where m = dO/dt. 
By integrating equation (7) we have 
Netlog {Loe O)/(T =O) halts (8). 


The values of ¢ and log, {(T) — ©)/(T — ©)} were plotted and X was found to be 
15"0 sec. The value of m was obtained from the appropriate G line and hence the 
product mA, This gave the lag of the thermometer in centimetres of bridge wire, 
and was subtracted from the values of G, before the corresponding Fy values were 
read from the G line. 


§10. SUMMARY OF RESULTS 


The values obtained for the water equivalent are given in the seventh column 
of Table 2, which gives also the experimental results from which the water 
equivalent was calculated. 


Table 2. Water equivalent of Dewar flask 


I 2 3 4 5 6 u 
t 76:6 001240 35°70 29'2 18°23 19'0 130 
774 0°01223 35°15 29°4 17°96 18°4 129 
24°8 0:03976 114°7 38°0 18°93 19°4 131 
BarT 0'02887 83°51 aga 18°13 18°4 133 
342 002840 82°12 BaP ly 18°25 18°6 me 
44:0 o°02121 61°05 40°7 18°44 19°6 130 
57°38 oo1614 46°41 38°8 19°24 20'4. 129 
70'6 001271 36°21 37°6 18°30 20'5 123 
109°5 0'00764. 21°99 36°2 19°60 20°4 130 
65°2 001431 PATS 34:0 19'10 18°7 129 
87°4 0°01023 29°43 34°5 19°29 19°5 129 
71°6 O'01231 35°88 39°2 20°74 2135 138 
519 001809 52°43 36°7 19°46 183 134 
40°0 0°02395 69°42 37°5 19'09 18:8 134 


Column 1 gives the value /% as obtained from the G curve, in sec. per cm, of 
bridge wire; 2, the rate of rise of temperature due to electrical energy in deg. C.; 
3, the electrical power used in watts; 4, the rate of stirring in rev. per sec.; 5, the 
mean water bath temperature in degrees C.; 6, the mean laboratory temperature 
in degrees C.; and 7, the water equivalent in gm. calories per degree Cc, 
16-2 
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COMPARISON WITH RESULTS OBTAINED BY OTHER OBSERVERS 


a ory 
Griffiths* in his experiments on aniline found that the stirring correction was 
about 1/25 of the total heating effect. When water was used, he experienced trouble 
due to electrolysis; but his German silver wire was uncovered and there was a 
considerable p.d. across his heating coil. The resistance of the heating coil in the 
experiments herein described was constant to o-oo1 ohm. 
Griffithst experienced considerable difficulty in his experiments upon 7. His 


observation that the thermometer readings were not steady and that the water — 


equivalent was inconsistent with low-speed stirring are confirmed by the present 
experiments. He used a water-motor to stir his calorimeter, and was troubled 
by its variable speed. In a somewhat smaller degree, the author experienced this 
difficulty, and this fact may account for the difficulty in obtaining really good stirring 
corrections. In the ¥ experiments the rate of stirring used was of the order of 
2000 rev. per minute, and the heat due to it varied from 1/6 to 1/10 of the heat due 
to electrical energy only. The heat due to stirring under the same conditions varied 
by about 5 per cent. Griffiths attributes this to the variation in the heat generated 
at the bearings. This may have occurred in the present experiments also, though 
there is no direct evidence of it. 

W.R. and W. E. Bousfield{ used a Dewar flask of much larger capacity, about 
3000 c.c. A much larger amount of electrical energy was supplied to heat this up, 
yet only a single blade propeller revolving at 340 rev. per minute was used. The 
author has found that such apparatus does little else than cause the mass of water 
to circulate en bloc. Its large diameter renders it very susceptible to loss of heat 
by evaporation and similar causes, and thus makes for some uncertainty in the 
results. Indeed, the experimental evidence obtained by both Griffiths and the 
author casts considerable doubt upon the efficiency of the stirring and hence upon 
the accuracy of the results obtained. The water equivalent in the Bousfields’ 
experiments was 1/6 of the thermal capacity of the calorimeter and its contents, 
and the greatest discrepancy in the water equivalent was about 1 per cent. 


§12. DISCUSSION OF RESULTS 


The experiments herein described show how important it is to ensure that the 
stirring is really mixing the contents of the calorimeter thoroughly. It appears 
that the water equivalent is independent of any of the quantities measured, and 
that the mean value of a number of determinations will accordingly be the most 
probable value. The experiments serve to demonstrate that, provided the necessary 
precautions be taken, the water equivalent can be measured with an accuracy of 
about r per cent., so that calorimetric experiments should be accurate to o- 5 per cent. 


ne aie the accuracy depending upon what fraction the water equivalent is of 
the whole. 


* Proc. Phys. Soc, 13, 234 (1894). t+ Phil. Trans, 184 A, 361 (18 
} Phil. Trans, 211 A, 199 (1911). Seas 
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DISCUSSION 


Mr R. S. Wutepte congratulated the author on a thorough and useful piece of 
work. The efficiency of stirring was an important requirement which the author 
had successfully met; failure in this respect was a common source of error in such 
experiments. He agreed with the author’s finding that Dewar flasks are of little 
use for calorimetry; further, their heat-retaining properties vary very widely from 
one flask to another. 


Dr H. R. Lana. I wish to emphasise that the inner vessel is not copper-lined, 
but is itself actually made of thin copper. It was intended that the sealing wax 
joint R, Fig. 1, should eventually be replaced by a fused copper to glass joint, such 
as is made for large wireless valves. 

In many experiments on the specific heat of oils, made of recent years in America, 
Dewar flasks have been used, and, despite the frequent remarks by other observers 
as to the vital importance of efficient stirring, insufficient attention has been paid to 
this matter. I would go further than Mr Whipple and not only urge greater atten- 
tion to the point, but suggest that some such tests as the plotting of ‘““G” lines 
should be carried out in every case so as to ensure that the stirring is satisfactory. 


Mr J. H. Brinkwortu. Thanks are due to the author for this work which adds 
very materially to our knowledge of the vagaries in behaviour exhibited by vacuum 
vessels. May I emphasise the importance of always having an arrangement whereby 
the quality of the vacuum can be continuously observed during experimental 
measurements. To illustrate this I show some figures relating to heat losses observed 
in the determination of the specific heat of steam near 100° C. by the continuous 
flow method. The diagram indicates that in Callendar’s experiment the glass 
vacuum jackets always cracked. I experienced the same kind of trouble when using 
glass calorimeters myself, and in consequence I designed the silica flow tube. It 
should be noted that the effect of silvering is quite unimportant when there is no 
vacuum. 

With more special reference to the paper under consideration, it would be very 
interesting to know if sporadic variations in the value of the water equivalent 
would be observed if experiments were made on all-glass vessels in which the per- 
fection of the vacuum was maintained: and would the water equivalent vary with 


the depth of water in the inner tube? 
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Diagram illustrating Mr Brinkworth’s remarks. 


AUTHOR’S REPLY. I agree with Dr Lang’s remarks concerning the importance 
of actually testing the efficiency of stirring. Wherever possible, a quantitative test — 
should be applied, for I have found that visual tests, such as the dropping of ink 
into water which is being stirred, are unreliable. It was because of Mr Brinkworth’s — 
work upon the heat loss from vacuum jackets that I maintained a high vacuum 
between the walls of the calorimeter. My experiments have shown the necessity 
for constantly testing the perfection of the vacuum. It is possible that the variation : 
of heat-retaining properties among flasks of the same type, observed by Mr Whipple, 
may be due to variations in the degree of perfection of the vacuum. 

I have no reliable evidence concerning the effect of the depth of water in the — 
vessel upon the water equivalent, since this depth was fixed by the stirring apparatus. 


So far as I am aware, no reliable experiments have been made with an all-glass — 
flask in which a high vacuum has been maintained. 


239 


THE MASSES OF THE PROTON AND ELECTRON 
By H. T. FLINT, University of London, King’s College. 
Received February 4, 1930. Read and discussed March 14, 1930. 


ABSTRACT. An application of a principle proposed by the author and described as 
the principle of minimum proper time is made to obtain an equation for the ratio of the 
‘masses of the proton and electron. The result obtained is in agreement with the experi- 
mental value: The possibility of explaining the asymmetry in regard to the masses by 
means of a metrical consideration is discussed. 


exist between certain physical quantities which have hitherto been regarded 
as independent. This has led to a search for these relations and we have as 
‘a notable example Eddington’s attempt to express the fundamental charge e in 
terms of A and c. We will consider here another relation. 

It has been suggested* that the world line of an electron is to be regarded as 
rade up of fundamental units of length of magnitude h/m c, where m, is the elec- 
tronic mass; similarly the world line of the proton is composed of units of length 
h/M,c, and that in the study of these fundamental masses no shorter length asso- 
ciated with their world lines will ever be revealed. 

Perhaps the most convenient expression of this theory is by means of a principle 
of minimum proper time which states that in association with the electron (or 
proton) no proper time less than h/mc? (or h/Myc?) will be observed. This is, of 
course, different from the Bohr-Heisenberg uncertainty principle which states that 
the position of an electron can be determined as accurately as we wish but that 
the momentum can then be determined only to an order of accuracy given by the 
equation AgAp ~ h, where Ag denotes the error in the determination of position, 
Ap that in momentum. It has, however, been shown that in atomic dimensions 
the application of the principle leads to no contradiction to the uncertainty principle. 

In a recent paper Fiirtht arrives at the principle by using the uncertainty 
relation and modifying it by assuming that the electron cannot be located as 
accurately as we please. By means of a number of rather arbitrary assumptions 
he arrives at a value of the ratio M,/m) and it is the determination of this ratio 
that we make the subject of the present paper. We shall approach it purely in the 
light of the principle of minimum proper time. 

Fiirth states that the value h/mpc is to be regarded as the radius of the electron, 
- but as this value is 2000 times the accepted value he inclines to the belief that the 
principle can apply only to neutral masses. There is no reason to regard h/myc 
as the electronic radius nor to limit it to neutral masses. The establishment of the 


T= new quantum theory in its later developments indicates that relations 


* Proc. R. S. A, 117, 630 (1928). + Phys. Zeit. 30, 895 (1929). 
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principle of minimum proper time shows that it is true for charged or uncharged 
masses and there is no suggestion that the fundamental length means more than 
a length along the world line. Fiirth’s results are, however, very interesting and 
we can remove his assumption. ; 

Let 7, denote the minimum proper time h/m,c*. If it were possible for the mass" 
my to be converted to radiation we could presumably calculate the frequency of [ 
that radiation by means of the equation m,c? = hy, and if the proper time is not less 
than 7, the radiation arising from this transformation of matter will have v as a_ 
maximum frequency. The proper time is the ordinary time in the system in which 
the electron is at rest, so that the frequency is calculated in that system. 

It is not possible to contemplate this transformation in the case of an electron 
alone, since this would require the disappearance of charge, but we may suppose 
it to occur in the case of an electron and a proton. In this case we consider the 
two bodies as a single one of mass (M, + m). The new quantum theory teaches 
that it is necessary to consider phenomena as possessing a double-sided character, 
one aspect being that in which we contemplate particles and in the other we con- 
template wavés; so that we can consider the problem as that of a particle of mass 
(M, + m,) and as that of radiation of a certain wave-length. The principle of 
minimum proper time applied to this union of electron and proton shows that the 
maximum frequency possible is v) where | 


Vy = (M, + m) c/h, 
or the minimum wave-length is A, where 
Ay = h/(Mo + mp) ¢. ; 


Now the wave-length in radiation is a unit corresponding to the unit along the 
world line of the particle. If we enquire what this unit could correspond to in the 
case considered, the suggestion at once arises that it must correspond to the distance 
between the centres of the electron and proton. Our method of approach to this 
point is rather different from that of Fiirth but this suggestion is substantially due 
to him. 

The older theories give as the value of the radius of the electron 7, = ke®/mgc® 
and of the proton Ry = ke®/M,c?. The value of k depends upon the distribution of 
the charge, and the values k = } and k= 2 have been considered in the classical 
theory. If, however, the charge is regarded in the light of the new quantum 
theory it becomes difficult to consider an electron or proton as a sharply defined 
structure, and we can only speak of a radius as an ‘‘equivalent” one. This Firth — 
does, and with a rather arbitrary choice of equivalent radius he finds k = 1 5/32. 7 


We can, however, still speak of a distance between the centres of the two bodies 
and we may write for it the value: 


uk (L tg) | 


since the dimensions of the charges appear to be proportional to e?/c? x mass. : 
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We have thus: h es = ob, ng ) 
(M,+m)c  c? \my | My 
or writing ' p= M,/m, 


e+ (1/p) + 2 = helke. 


The occurrence of hc/e? is noteworthy. It occurs in Eddington’s calculation 
of e and in the determination of the number of chemical elements by the application 
of the principle of minimum proper time; and we note that it is a pure number, 
having no dimensions. 

We cannot proceed further to examine whether the equation gives the experi- 
mental value for « unless we know k accurately. If we take the classical view and 
write k = $ or 3, p is of the right order but the value is too small. The result is, 
however, striking, for we have an equation in » thus reducing the number of 
fundamental constants by one and the failure to obtain the correct value is explained 
by our ignorance of k. If we accept Fiirth’s assumption that the equivalent radius 
must be calculated on the assumption that it is the radius of the sphere within 
which all the charge must be enclosed so that it produces the same moment as 

the charge in the actual distribution assumed in the new quantum theory, k = 15/32 
and the value obtained for yz (1838-2) is then in agreement with the experimental 
value (1838-3). 

One of the most interesting and puzzling things in atomic physics is the 
asymmetry with regard to mass in the case of proton and electron. The question 
is, why is the positive charge associated with a mass so widely different from that 
associated with the negative charge? In theoretical considerations the asymmetry 
is very striking. In one line of attack, that of the five-dimensional theory, we meet 
with a constant @ which may have the values + e/m)c. ‘Thus the theory seems to 
account for the occurrence of positive and negative charges, but there is no sug- 
gestion that m, has more than one value. It would be a great step if we could in 
some way relate mand M,, and this makes it worth while to consider attempts in 
this direction such as that given above. 

We will consider the problem from another point of view, though the asymmetry 
is replaced by another. Nevertheless a change in the point of view may be useful. 
It has been suggested that atomic phenomena are indicative of a special kind of 
-metric preferred in nature. In this we are following the way opened up by Weyl 
and Eddington in their inclusion of electromagnetism with gravitation into a space- 
metric system. We take the view that it is not in electromagnetic phenomena that 
the metric is to be found, but in quantum phenomena. Now, if we take the view 
that there is one standard of measurement for the proton and another for the 
electron, we can do without the introduction of a second mass, and retain one mass 
for both electron and proton. But we introduce two scales of measurement and 
so do not alter the number of constants. Some may prefer to introduce different 
masses directly rather than different scales, but there is more in the change than 
this, for the metrical view has certain other advantages. We can only decide which 


assumption to adopt when we discover that one idea is more fruitful than the UNS 
NI D _ 
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We know from the general theory of relativity that the space near a proton must 
be much more strongly curved than that near an electron and we can readily 
appreciate that the metric of space also might be markedly different in the two- 


cases. 
This point may be explained by reference again to the world lines of the proton 


and electron. These may be regarded as made of elements of lengths /, ‘M,c and h/mye. 
We have, as it were, two scales, one more finely divided than the other. Now the 
view under consideration is that this difference in magnitude is only apparent, 
the length 4/M,c in the proton space being regarded as equivalent physically to 
h|mc in the electron space. We have a kind of squeezing together of space in the 
proton relatively to that in the electron. The reason for this statement of equivalence 
of these two lengths is that in a parallel displacement in the proton region from one 
end of the element //M,c to the other we have the same change in length per unit 
length as in the electron space from one end of the element h/mc to the other. 

If we judge the phenomena concerning electrons and protons from the same 
point of view and apply the same metrical considerations to each we shall conclude 
that a proton moves more slowly than an electron under the same conditions and 
so estimate its acceleration at too low a value and attribute to it greater inertia. 
We ought to regard the unit in the proton space not as of length h/ Myc but as h/mpgc, 
since the unit of length means physically so much more in that space. Thus our 
constant @ is + e/m) c, but with the change of sign we must change our metric. 

There is another point. If this view is right the proton space is a kind of miniature 
of the electron space, and we should expect the radii of the two bodies to appear in ~ 
the ratio of the units of the scales, i.e. inversely as their masses. This is in fact 
what is assumed about these bodies for, as we have seen, . 


2 2 
e e , 
and = k——.. 
Myc Ks M,¢e 
This view is offered as no more than a suggestion of a new way of looking at 
one of the problems occupying us nowadays. 


Taek 


DISCUSSION 


Prof. W. Witson: Dr Flint’s paper contains a further application of a principle 
which he proposed some time ago; he has already made some important applica-_ 
tions of it, one of the most interesting being an estimation of the upper limit for 
the atomic numbers of the elements. I think it is no mere coincidence that this 
turns out to be about 97, while the atomic number of uranium is g2. It is hardly 
surprising that the distance h/mpc should have a special significance. A similar _ 
expression, namely h/mv, where mv is the momentum of the electron, is well — 
known to represent the experimentally determined value of the wave-length of — 
the waves associated with the electron. 


As Dr Flint remarks, his principle differs from the Bohr-Heisenberg uncer- 
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tainty principle which asserts that the errors, Ag and Ap, in the determination of 
the position and momentum of one electron, are governed by the relation 


Ag Ap sh. 


I feel some doubt as to whether the latter principle is correct. It is true, of course, 
that any device for accurately determining the position of an electron at some instant 
will cause it to suffer a violent change of momentum; but what reason prohibits 
us from employing a device which gives us accurately the velocity (only) of the 
electron and, while this is in operation, from employing some other device, which is 
adapted to determine its position accurately? 

The asymmetry in the masses of the electron and proton is one of the most 
provoking of physical facts. It is one of the things we should like to be able to 
trace back to some still more fundamental fact, and Dr Flint’s tentative suggestions 
are of great interest. 


Mr J. E. Carrurop questioned whether the author was justified in adding 
together the masses of the proton and electron in his calculation, since there would 
be a change of energy involved in their near approach. 


Mr I. O. Grirritus. I listened with interest to Dr Flint’s account of his de- 
duction of a fundamental length connected with the electron. In Fiirth’s paper, to 
which he refers, a value of a certain expression is obtained which is 2000 times the 
accepted value of the radius of the electron. As this expression involves the mass of 
the electron linearly, it was perhaps natural to introduce the atomic mass in ad- 
dition to the electronic mass, thus introducing a factor of approximately 1840 
which goes some way to counterbalance the factor of 2000. The agreement which 
is claimed would however have been more striking if it could have been obtained 
with a simpler fraction than 15/32 for R. 


DrW.N. Bonn. I should like to express my appreciation of this further attempt 
by Dr Flint to relate theoretically experimental constants that were formerly 
unrelated. The comparative smallness of the number 1838, where numbers of the 
order 10! and 10-!° may be concerned, indicates that the relationship might be 
expected to be calculable theoretically. The calculation given in the paper deduces 
the relationship in a very interesting way, which I am not in a position to criticise. 
As regards the numerical value, according to Birge* the ratio M,/m, has two values. 
In any case, numerical agreement to nearer than one part in 1000 cannot be signi- 
ficant, as it is limited by the experimental errors in the value of hc/e?. Apparently 
the experimental value 27 (137:29 + o-11) has been used by the author, not the 
value 27 x 137 suggested by Eddington. Personally I feel that there is now enough 
evidence to accept the latter value. However, all lack of exact numerical agreement 
could be attributed to the uncertainty as to the value of k; and the doubt about 


this value seems the only inelegant part of the paper. 


. 


* R. T. Birge, Phys. Rev. Supp. 1, 47; 48 (1929). 
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Capt. C. W. Hume. In attributing the difference in mass between the proto 
and electron to a change of metric, the author has invoked the increase in curvatur 
of space-time in the neighbourhood of the proton. But is not the difference 
mass already implicit in the difference of curvature? The proposed change seems” 
to be only a verbal one. : 


Mr T. Smits said that he felt the greatest admiration for those developments 
of modern physics which skirted the border-line between the physical and the 
psychical, but he felt that a dictionary of the terms employed was badly needed | 
The author has used a phrase “physically equivalent,” for instance, which the — 
speaker was quite unable to understand. 


AUTHOR’S REPLY. It is not necessary to go into detail in this case with regard to” ) 
the process of union of an electron and proton to form radiation. We must for o 
present purpose take account of the whole mass concerned, wherever it may 
distributed. 

The reference to the curvature in the neighbourhood of the proton was 
by way of illustration; it is not desired to draw any conclusion from it, but only to 
point out that, as there are such great changes in the geometry in this region | 
space, we need not be surprised to find a metrical change also. ) 

The meaning intended in the use of the term “‘ physically equivalent” was tha 
any change taking place in the electron region over the space of one unit length — 
appropriate to this region may be regarded as equal to that which takes place in the 
protonic region along the unit appropriate to it. We may think of an observer _ 
attached to the electron and compare his experience with that of an cbectvell 
attached to the proton. In the same electromagnetic field the experience of each — 
would be the same along the units appropriate to them. 
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A HARMONIC ANALYSER 
By J. HARVEY, A.R.C.S., B.Sc., City and Guilds Engineering College 


Yommunicated by Dr E. Mallett fanuary 28, 1930. Read and discussed March 14, 1930. 


BSTRACT.- The paper describes a device for finding a given number of the Fourier 
omponents of a function by means of successive weighings. 


SNR OD Cl LOIN 


OURIER’S theorem states that if a function f(x) is periodic in intervals of x 
equal to 2/, then f (x) can be expressed as the series 


“al, 277 
OT oy ue COS ys XT A 
+ by sing «+ by sin |X + «ss 
here a, a, d,--. b;, b ... are constant coefficients. 


If the form of the function f (x) is not given, but the lengths yo, y,, v2, --+ Ven—1 


f the 2 ordinates of the curve y = f (x) for the values 0, - » «=. (22 — 1) : of x 


re known, a finite series 
T 277 T 
Ay + G, COS GH + Gy COS H+... + Ay COSN FH 


l l 


+b, sin” « + by sin & +... + dy sin (n — 148 


l l 


f 2n terms can be found representing a function ¢ (x), which has the same values 
oy Vi» Yas «+» Yon—1 a8 f (x) for the corresponding values of x. 


Then a = = [yo +41 t+ Vo + «+ Yona), 
I V7 2V7 VT 
a, = 7 lye +94 cos — + Jz Cos = — +... + Von-4 CO8 (22 — 7], 
here 7 =1,2,...(%— 1). 
a = [yo + 1 COS 7 + Yo COS 277 + ... + Yana COS (an — 1) 7], 


I HT Ora: ‘ VT 
= == —-+... pest (2h — 1) 
b, = Ely sin 7 + yy sin ap + Von—1 ( ike , 

* 
here = Ty 2, ao: (i — 1)" 


* See Whittaker and Robinson, Calculus of Observations. 
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By means of the harmonic analyser here described the summation of the series 
for the coefficients when 2n = 12 can be effected. On the principle employed ar 


instrument may be constructed for evaluating the series for the coefficients corre- | 


sponding to any value of 27. 


Bp 4 5h Os 8h co 
SBS pas 


§2. THE INSTRUMENT : 


From a bar A, supported on horizontal knife-edges, a bar V is suspendec 
which can rotate in a vertical plane called the principal plane, perpendicular to 
the knife-edges. Along the bar V’ there is a group of twelve (= 2) circular dises 
D coaxial with the rod and supported by a plate P rigidly fixed to V’. Each dis¢ 
is capable of rotation about V’. Fixed to each disc and projecting radially in opposite 
directions are two similar bars H and H’. The construction is such that any o: 
disc may be turned about V’ into a position in which its bars H and H’ make with 
the principal plane any angle which is a multiple of 30° (= z/n)*. 


* There are several methods of achieving this. One such is as follows: On the upper face ¢ 
the plate P and of each disc, twelve radial grooves at angles of 30° with each other are cut. In 
small cylindrical hole parallel with V cut in each disc upwards from the lower face, a plunger plays 
against a fine spring at the upper end. The lower end of the plunger is wedge-shaped and can ente 
any one of the twelve grooves in the disc beneath and give a slight check to the rotation of one dise 
with respect to the other at every 30°. The check is just sufficient to keep the discs, with their bars 
when loaded, relatively fixed during free oscillations of the whole system about the knife-edges. 
The plungers are placed in such positions in the’ discs that’ all the bars H and H’ can be brought 
together into the principal plane. The rotation of a bar from this plane through an angle which is 
multiple of 30° can then be made. 
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On the bar H attached to each disc is a ‘unit mass” U symmetrically disposed 
bout the axis of the bar, and having sufficient friction with the bar to keep it in 
ny position in which it has been placed. When a unit mass is situated in a certain 
osition near the disc, as in the diagram, called the zero position, the bar V being 
fertical, the force of gravity on the unit mass produces a moment about the knife- 
dges. This moment is counter-balanced by that produced by a mass m fixed to 

the bar H’ opposite. Hence when the unit masses on all the bars are at zero 
ositions, the system with the bars in any relative horizontal positions will be in 
equilibrium about the -knife-edges and V will be vertical. The vertical position 
of V is indicated by the coincidence of the prolongation of V, as an index J, with 
he zero of a fixed scale. 

Fixed to the horizontal bar A carrying the knife-edges, at right angles to them 
d to the bar V, are two similar arms B of the same length, opposite each other, 
onstituting a beam as in a balance for weighing. From these arms are hung scale- 
ans, S, which balance with respect to the knife-edges. A counterpoise C movable 
n the bar V can be adjusted in position so as to give sensitivity to the system. 

On the base of the supporting frame there is a mechanism by which the bar V 
carrying the discs and bars is fixed or released: by a horizontal movement of a lever 
© jaws are opened or closed on the bar, while a vertical movement of the lever 
a certain position locks or unlocks the jaws. 


§3. DETERMINATION OF THE COEFFICIENTS 


The twelve unit masses are moved along the bars, taken in order, from the zero 
ositions through distances yy, y,, V2, --- Yu, corresponding to the lengths of twelve 
quidistant ordinates of the curve y = f (x) in the range x = 0 to « = 2/. The bars 
are then brought into the principal plane so that lengths representing positive 
ordinates are in one sense from V while those representing negative ordinates are 
in the opposite sense. The bar V will then, in general, be deflected from the vertical 
position, positively or negatively. The vertical position is restored by placing mass 
in one of the scale-pans. The moment of the weight of this mass about the knife- 
edges is equal to the sum of the moments of the weights of the unit masses, which is 


Yor M+ 2+ +) + Iu 


Hence the mass added to the scale pan is a measure of a). 

In general, to obtain a coefficient a,, the bars in the order referred to are turned 
Bo as to make angles 0°, 7 x 30°, 7 x 60°,...7 X 330° respectively with the 
principal plane. Mass is then added to one of the scale-pans to maintain the bar V 
in the vertical position. The weights of the unit masses now exert about the knife- 


edges the moment 
Yo +31 C08 (7 X 30°) + Yo Cos (ar x 30°) + ... + 41 COS (IIT X 30°). 


Hence the mass in the scale-pan is a measure of a,. 
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To obtain the coefficient 5, the bars are kept in the same relative positions as 
for a,, and the whole system is turned through a right angle about the bar V, 
If V is maintained vertical by the placing of mass in one of the scale-pans, the 
moment about the knife-edges of the weights of the unit masses is now 


y, sin (r X 30°) + yp sin (ar x 30°) + ... + Yu Sin (IIT X 30°), 
and the mass in the scale-pan is a measure of 5,. 

On a model actually constructed each unit mass is of 50 grams and has a range 
of movement of 16 cm. on the bar H. Ordinates may be set off by callipers te 
within o-s5 mm. Each arm B is 30 cm. long, so that a movement of a unit mass 
through o:5 mm. along the bar produces a moment about the knife-edges w 
can be equilibrated by a mass of 0-08 gram approximately in the scale-pan. The 
instrument can be made sensitive enough to weigh to o-2 gram. If the ordinates 
be measured along the bars in centimetres, the mass placed in a scale-pan to main 
tain the bar V in the vertical position should be divided by ro to give the coefficient 
sought, except in the cases of a and a, when the mass should be divided by 20. 


§4. ACCURACY 
To test the accuracy of the model, the twelve ordinates of the curve 
y = 3 — cos x — 3 cos 2x + 3 cos 3x 
+ 18 sin x + 9sin 2x + 3 sin 4x 


for the values 0, 7/6, 27/6, ... 117/6 of x were calculated. These ordinates were 
set off along the bars H of the instrument and the coefficients of the cosines and _ 


By calculation 


1 By instrument By instrument 
from series 


(1) (2) 


Coefficient 


ay 3°00 2°04 2°93 
{ ay —4°50 —4°44 —443 
by 3°75 S75 3°68 
¢ —3'00 —2°94 —2°97 
by 2°25 2°20 2°23 
| as I*50 1°47 1-48 
bs 0°00 0°03 O-or 
15 at Sera 0°00 
by O75 o'75 o'77 
‘e 0°00 o-or o°or 
bs 0°00 0°00 ool 
0°00 0°00 O°or 


The instrument should be sufficiently accurate for the analysis of many 
of geophysical and meteorological phenomena. 
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5. EXTENSION TO MORE THAN TWELVE FOURIER COMPONENTS 


An instrument consisting of twelve discs and bars may be employed to give the 
oefficients of the terms of a sine and cosine series for a number 2” which is a 
aultiple of twelve, if checks are provided so that the discs can be turned relatively 
hrough corresponding angles. Thus when 27 = 24, twenty-four equidistant 
rdinates yy, V1, V2, +-- Y23 Of the given curve over the range of periodicity are 
neasured, and relative movements of one disc over another through multiples of 
5°, or =7a/12, must be made. For the determination of a coefficient a,, the unit 
nasses on the twelve bars taken in order are moved distances yo, 91,25 --- Vu 
rom the zero positions. The bars are then turned through angles 


Opti ese) (oP tS. )yesai( LIP ES +) 


espectively from the principal plane. A mass M, is added to a scale-pan to main- 

in the bar V in the vertical position. Next, the unit masses are moved so that 
hey are at distances yj, Vis, --» 23 along the bars taken in the same order as before 
rom the zero positions, and the bars are turned so that they make angles 


[eore ree), (iar 15); ..< (237 X 15.) 

espectively with the principal plane. Mass M, is added to a scale-pan to main- 
ain V vertical. Then M,+ M, is a measure of a,. The coefficient 5, may be 
btained similarly from the masses required to maintain V vertical when each 
f the two systems of bars, as arranged for a,, is turned through a right angle. 
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DISCUSSION 


Mr J. H. Awsery congratulated the author on his ingenious invention and 
asked how much time it would require for an analysis, as compared with the time 
taken by arithmetical methods. The data given in the paper showed that the 
results were sufficiently accurate for practical purposes. 


Mr J. E. Catturop. Speaking from the point of view of a teacher of physics, 
I should like to congratulate the author on a piece of apparatus which should be of 
great use and interest in a laboratory for teaching purposes. I suppose the use of 
callipers could be eliminated by suitable graduation of the arms carrying the 
movable weights. 


The PRESIDENT said that the instrument, by its direct attack on Fourier analysis, 
would give students a very clear view of the principles of the subject. Many persons 
17 
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must have interested themselves in the problem, but probably the instrumen 
shown at the meeting was the first that had been actually constructed, at all even 
with such excellent workmanship and accuracy. 


AUTHOR’s REPLY. The bars carrying the movable weights could be graduated. 
but it would, I think, be more difficult to set off the ordinates along the bars by t 
graduations than by the aid of callipers. The pair of callipers, with its vernier, can 
be turned in the hand towards the light and the length of an ordinate set oe on i 
accurately and rapidly. 

In the estimation of the coefficients of the function given in the table in th 
paper, the time taken by using the instrument was under twenty-five minutes, 
single weighing being made for each coefficient, and I do not think that much ti 
would be saved by this method in comparison with certain arithmetical meth 
in the case of twelve ordinates. 


251 


A METHOD OF COMPARING SMALL MAGNETIC 
BU SCHR UIBIEITIES 


By R. A. FEREDAY, B.Sc., East London College. 


Communicated February 26, 1930. Read and discussed March 14, 1930. 


BST: RACT. A method of the non-uniform field type is described. It possesses the 
specific advantage that the force exerted on the specimen is constant over a wide region 


of the field, so that a careful setting to a position of maximum force is rendered un- 
necessary. : 


Ste LN © DW COR Yi 


HE present communication gives an account of a method which has been 
developed at the suggestion of Prof. C. H. Lees, and in consultation with 
him, for the determination of small magnetic susceptibilities, and departs 
considerably from general practice in measurements of this type. The apparatus 
was needed for an investigation in which a high sensitivity was required when a 
specimen of but a few milligrams was being used, and for this reason, and since 
relative measurements only were required, a method of the non-uniform field type 
was chosen. 


§2. THE NON-UNIFORM FIELD METHOD 


The force due to a field H, acting in any direction s, on a small body of volume v 
and susceptibility x, that of the surrounding medium being kg, is given by 


r oH 
J, = (K — Ko) vH ~ ieee: (a): 
Thus, if two bodies are brought in succession to the same part of the field, 
and the forces on them are compared, we shall have: 


ieee OR ls Wht i. vhs (2). 
Ke — Ky fo 
In the derivation of (1), it is assumed that the force on all parts of the body is the 
same, a condition which it is almost impossible to realise in practice with a specimen 
of reasonable size. Equation (2) will, however, hold accurately for larger specimens, 
provided they have the same size and shape, and are brought to the same position 
in the field*. The method is therefore much better suited for relative than for 
absolute measurements, but each specimen must be placed at the same point in 
the field, and:this constitutes the chief difficulty of the method. 
, Curie and later observers used the magnetic field between two inclined poles ; 
the direction of the field gradient, and hence that of the force exerted on the 


* Cf. Stoner, Magnetism and Atomic Structure, p. 39+ + Ann. Ch. Ph. 7, 289 (1895). 
17-2 
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specimen, being at right angles to that of the field itself. A curve is drawn showing 
the variation of force along the direction of the gradient, and it is found that at a 
certain point it passes through a fairly narrow maximum, on each side of which 1 
falls off sharply. In the comparing of susceptibilities, each specimen 1s placed at 
this point of maximum force, and to secure this, the electromagnet producing a 
field is usually mounted on some form of carriage, so that it can be displaced 
relatively to the specimen until the latter is in the desired position. It is, however, 
by no means easy to impart micrometric movements to a heavy electromagnet, | 
especially when the specimen is suspended in the interpole gap and must not be 
unduly disturbed. 


(a) (b) 
Fig. 1. Lines of force in converging fields. 


| 


Another form of the non-uniform field method was introduced by Wiedemann* _ 
but has not found favour with modern workers. In this method, the field, its — 
gradient, and the force to be measured, are all directed along the axis of symmetry 
of the field, which is produced by the pointed pole of an electromagnet. A suitable — 
field must be of one of the types of Fig. 1, which represents sections of the field 
by a plane passing through the axis of symmetry. 

In the case of a field of type 1 (a), both the field intensity and its gradient 
increase from left to right of the figure; hence the force exerted on a body placed — 
on the axis will be very much greater at the right than at the left of the figure, and 
if measurements were made in such a field it would be necessary to attend very 
carefully to the setting of the specimens. With a field of type 1 (6), on the other 
hand, the field decreases as its gradient increases, so that the variation of their 
product along the axis will be small. If these opposite variations could be made tay 
balance exactly, so that H x 0H/ox had a constant value at all points of the axis, 
then the force exerted on a given specimen would be independent of its position — 
so long as it remained on the axis, and the great experimental difficulty of setting 


the specimen would entirely disappear. It is on this conception that the method 
to be described is based. 7 
t 


* Pogg. Ann. 126, x (1865). 
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§3. THE MAGNETIC FIELD 


_ In an attempt to realise a field satisfying the condition H.0H/0x = constant 
along the axis of symmetry OX, the assumption was made that the pole pieces 
f the magnet could be used with such a field that the permeability of their material 
ould be high; the pole faces would then very nearly be magnetic equipotential 
surfaces. The problem thus reduces to the determination of equipotential surfaces 
in a field satisfying the above condition. Further, since the field is symmetrical 
about an axis, it is only necessary to determine the curves of section of the equi- 
potentials by a plane through this axis, and in the sequel, the term “ equipotential ” 
is used to describe such curves. Two systems of co-ordinates are used in the plane, 
one a polar system, in which the axis of symmetry is @ = 0, and the other a Cartesian 
system, the axis of symmetry being taken as OX and a perpendicular at the origin 
as OY. In each case, the centre of the interpole gap is taken as origin. The equi- 
potentials were first determined by a simple semi-graphical method, and the results 
were used to design a pair of pole pieces. With these pole pieces in the apparatus 
to be described, numerical data were obtained which served for guidance in the 
‘more accurate analysis, which is given in the next paragraph*. 

The important constants determined experimentally were: 


Field at origin (Hp) ... ee 1000 gauss. 
H.dH/ox a mA or 21>. 


Since the field is symmetrical about the axis, the magnetic potential at the 
point (7,4) may be expressed as 


V = V+ A,7P, (cos 0) + A,7?P, (cos 6) + Agr3P3 (cos 8) + ... «+++ (3). 


For any point on the axis of symmetry, the Legendre polynomial P,, (cos @) is equal 
to unity, so that the potential at any point on the axis is given by 


V = Vok Agta Agr? + Agra Ake cet (4). 
T he field along the axis at the point (r, 0) is therefore given by 
4 _ H =0V er = A, + 2Agr + 3Agr? + 44g? + ... 

sez Ast ae g Agr th AA ses ) setoes (5) 


(H, being the field at the origin), 
therefore 


H? = H,? + 4H) Azar + (6H As + 442") 7? + (8H, A, + 124, As) 7° + «.- 
2 a) = 4H, A, + 2 (6HyAs + 442") 7 + 3 (8H, Ay + 12A,A5) 7? + eee (6). 
‘Now we require to satisfy the condition that H OH /ér, ie. 4.0H?/or shall be 


* T am indebted to Prof. Lees for suggesting this method of attack, and for assistance in its 
development. 
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independent of 7, having a constant value of, say, c. This condition will be satisfied 
if the coefficients of powers of r in (6) all vanish; i.e. if 


a : c 
os (H*) = 4A, = 2c; te 4, = 3H,’ 
: 2A." Cc 
6H,A; + 44,? = 0; i.e. A; = a, eee 
: A a 
8H, A, + 124,A; = 0; 1c. A,g= — aa As ra © SHS” 


Substituting these values of the coefficients A,, A,, etc., in (3), we have, as the 
expression for the potential at the point (r, #) in a field which satisfies the required 
condition along the axis, 


re: e 
V = V, + HorP, (cos 8) + ag r?P, (cos #) — 6He PP, (cosG) Per noe (7), 
V = VY — Bt 2 =, c 3 es 
ne 7P; (cos 7) + spas P, (cos @) 6H,A” P, (cos @) be edn (7a). 


Substituting the numerical values characterising the experimental apparatus, i.e 
H, = 1000, ¢= 3 x 10°, we have 
Re: a 
C4 Rennes 4 OF 
The series will converge least rapidly when r and cos # have their maximum values 


The pole gap being, in the apparatus, 3-4 cm., 7 cannot exceed 1-7, so that when 
the series of (7 a) is converging least rapidly, its terms will have values 


= O-OI5. 


1-7 + 0°43 — 0:073 + .... 


It is thus clear that for an apparatus not too different from the one tested, the value — 
of (_V — V,)/H) may be expressed by the first two terms of (7 a), without serious 
error. 

We thus have 


V-—YV, 


ae BY lees alee 
i, = 7B, (cs &) ok Spratt a (cost) = meaaea 


He r? (3 cos?@—1) ...... (8). 


If different values are taken for V’ in this equation, we have a series of equations 
to the corresponding equipotentials. It will be convenient to use the same equation, 
referred to the Cartesian system of co-ordinates; its form in this case is 
V-YV, c 
— ae ry 
Hi, x tHe (2x2 — +?) veeeee(8 a), 


The equation may also be expressed in a third form. Suppose that the equipotential 
for some particular value of V cuts the axis of symmetry at (X, O). Then we must 
have | 
—V c 
——9 = X¥ + — = 
Ti, X + 4He (2X2 — 0), 
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o that the general equation may be written 


c C c 
X + -X2? =x + —_ x? — —_ y?, 
Zidnn elds 4H?” 
: — Vy. c : 
we write i="; , = 6, then the different forms of the equation to the 
Hy 4f1y 
quipotentials become respectively 


G7. COS 0 OF" (2.Cos* 0-1) ae (8), 
Fics OY i'l | a nee (8 a), 
CX ep X Ve c 2bx7)— bye ee wees (8 b), 


rom which successive equipotentials are determined by giving successive values 
o either a or X. | 

The equations just derived define equipotential surfaces of infinite extent. In 
fashioning a pole face to be an equipotential surface it is, however, necessary to 
use a limited area. The condition to be satisfied is that if the field is produced 
between two poles shaped as equipotential surfaces according to the above equa- 
tions, and separated by the appropriate distance, then their areas must bear such 
a ratio to each other that the same total flux passes through each. 

We shall satisfy this condition by terminating the sections of the equipotential 
surfaces at points which are on the same line of force. We have the equation to the 
equipotential (equation 8 b), 

a=x +b (2x7 — y?), 
-- o=dx+b (4x.dx — 2y.dy), 
(4bx + 1) dx = 2by.dy. 
A system of normal curves (lines of force) must satisfy 
(4bx + 1) dy = — 2by.dx, 
Nl ACER 
aby © 4bx +1 
(2b) logy + (46)~1 log (4x + 1) = const. 
2 log.y + log (4bx + 1) = const. 


O, 


y? (4bx + 1) = const. 
Thus, if one equipotential is terminated at the point (x,,y,), the other must be 
terminated at a point (x2, ys), related by the expression 
yr? (4bo + 1) = 92? (4b%2 + 1); 


ch Ls ec (9). 
; yo? 4bx, + 1 

If the desired interpole gap (2X), and the radius of the larger pole (91) are known, 
then equations (8b) and (9) suffice to determine a pair of pole pieces for any 


assumed value of b. 
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§4. EFFECT OF FINITE SIZE OF SPECIMEN 


Having determined a field distribution such that a point specimen placed on 
the axis experiences a force which is independent of its position on the axis, we 
have to find to what extent this constancy is affected when a body of finite size is 
considered. 

For a finite body symmetrical about the axis of symmetry of the field, the 
resultant force will still be directed along this axis, but that portion of it which 
is due to the effect of the field on the outer parts of the specimen may vary with 
position along the axis, and thus destroy the constancy of the total force. 

We have, from (8 a), the potential at any point (x, y); 


V =V,.+xH, + a, (2x2 — y®), 
0 


therefore the field, H,,, along a direction parallel to the axis OX is 


H, = ingles 

x Hy, : Fe ase 
att © 
ox Hy 


Thus the product H.dH/0x, which determines the force on the specimen, should 
be quite independent of y. This approximation is not close enough for the present 
purpose, and a more accurate calculation must be made in which a further term 
of (7) is taken into account. Thus we have 
Pa = 
V = V, + HorP, (cos 6) + 3H, r°P, (cos @) — HS r°P, (cos @) 
= th 52 9 — 2 : 
=V,+ oes 4 r? (3 cos? @ — 1) — r2He r® (5 cos* 6 — 3 cos @) 


— 3 2 2 c 2 
=V)+ HAyx + iH, (2x = y ) = 12H,° (2x8 —= 3xy?), 


Deak ey 3 se ore + a7 
Ha = He + rast + ast + saat + ee — rat + yt at 
+ a — ap 
0 (He) = WS apt t pat +20 Fae Beata Say ae 
ny Set =a am 


wy (Se aoe ( ee 
oy | Ox} Ae? \? He"): 
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Substituting numerical values for the present experimental apparatus, we have at 
the centre of the gap (x = 0) 


0 H 0H) oc 

As an example, consider a point at 5 mm. from the axis. The mean gradient over 
the interval is that at y = 0:25, which is 13500. The change in force is y.dy/dx = 6250. 
Since the value of c on the axis is 3:10°, the percentage change in its value at 5 mm. 
is thus only of the order 5 per cent. Thus, the force exerted on the outer part of 
_a specimen of 1 cm. diameter is 5 per cent. different from that exerted on the centre, 
so that the force exerted on the whole body will differ by only one or two per cent. 
from that on a point at its centre, and will thus be constant to within this degree 
of accuracy as the body is shifted to different parts of the axis. 

It must be remembered that in practice the force at a distance from the axis 
will be seriously affected by “fringing” effects, unless the magnet poles are of 
large diameter. It is, however, clear from the foregoing analysis that if the specimen 
is small compared with this diameter it is quite unnecessary to take elaborate 

precautions to bring it to a definite position in the field, since H.0H/cx will have 
a value which is appreciably constant over a very large region of the interpole gap. 


§5. EXPERIMENTAL REALISATION 


As has already been stated, a pair of poles has been designed by a less accurate 
method; these have been used with a simple torsion balance for susceptibility 
“measurements, as it was found that their shape does not differ seriously from that 


Fig. 2. Section of pole pieces (full size). 


indicated by the more accurate theory. The force obtained in the pole gap is not 
accurately constant, but the constant nature of the susceptibility measurements 
obtained shows that it does not vary seriously. As a precautionary measure, how- 
ever, the specimens are set roughly to the same point on the axis when suscepti- 
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bilities are being compared; for this purpose, a small pointed plug has been let 
into the centre of one pole face, and the specimens are set to a standard distance 
from this. 

Any attempt to realise the theoretical conditions more accurately has been 
postponed until experience of the method has been gained from a series of measure- 
ments with the present apparatus. 

The pole pieces in use are shown in Fig. 2, from which it will be seen that the 
smaller one has been elongated in order to reduce the disturbing effect of the outer 
part of the pole core; this is necessary since the poles of the electromagnet are of 
equal diameters. 


§6. APPARATUS 


The electromagnet. Details of this are shown in Fig. 3. The coils have about 
2500 turns of heavy-gauge copper wire, and can be supplied with current up to 
10 amp. from accumulators. As will be seen from the figure, copper tubes, through 


"WOODEN SUPPORT 
FOR 


T-TUBE 
BRASS BAR 


uN 
SUPPORT 
FOR 


MAGNET 


“ 


COOLING TUBES 


Fig. 3. Electromagnet (partly sectioned) showing method of supporting torsion balance 


which cold water can be circulated, are disposed so as to prevent heating of the 
poles or the air gap by the windings. A massive brass bar is bolted across the top 
of the magnet to keep the poles at a fixed distance apart, and this serves as a support 
for the torsion balance used in measuring the force on the specimen. 
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The torsion balance. The whole of this system is contained in a T-tube of 1-in. 
glass tubing, each arm of which is about 12 in. long (Fig. 4); it is held in wooden 
clamps attached to the brass bar just mentioned. The tube is provided with a 
brass torsion head H/ fitting into its upper end, through which passes a brass rod R 
which is secured by the screw G and can be adjusted vertically to suit suspensions 
of different lengths. The lower end of this rod is filed down to a diametral plane 
to which the suspension is attached. 'To the lower end of the suspension, at O, 
is fixed the torsion arm A, which is made of glass rod about 0-5 mm. in thickness, 
except for a short thicker portion at O. The lengths of the two parts AO, OA, of 


MAGNET 
POLES 


Mi 


Fig. 4. The torsion balance. 


the arm are about 5 cm. and 20 cm. respectively. At the end of the shorter part is 
a small glass cup C, which takes the specimen; at the other end there is a pointer P 
of very fine glass, which is used in conjunction with a reading microscope to measure 
deflections of the system in a manner previously described *. 

The difference in weights of the two arms is balanced by a spiral of lead wire L, 
this material being chosen for its low susceptibility and high density. Small changes 
in balance of the system due to difference in weights of the specimens are com- 
pensated by shifting of the copper wire rider W. . 

The ends of the horizontal part of the tube are closed by two closely fitting 
paper caps Kj, Ky, provided with microscope cover slips as windows, through 
which microscope observations can be made without danger of error. 


* Fourn. Sct. Inst. 6, 302 (1929). 
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The substance examined is enclosed in a small glass bulb from 2 to 3 mm. in 
diameter which fits into the cup C. The clamps hold the T-tube at such a distance 
from the electromagnet that, when the tube is adjusted at the correct height, the 
cup lies on the axis of symmetry of the field. ‘ ae i 

Considerable difficulty was experienced in finding a suitable non-magnetic 
adhesive for the attachment of the suspension and lead counterweight, as any force 
exerted by the field on the adhesive, especially that used in attaching the suspension 
to the torsion arm, caused a deflection due to a motion of translation and not to a 
rotation about the suspension. Shellac varnish, applied thinly and “baked” by 
means of an electrically heated wire, was ultimately found satisfactory. 

The apparatus proved so sensitive that for the observations made so far, on 
somewhat strongly magnetic salts, it has been found necessary to provide a bifilar 
suspension of relatively high restoring couple. A silk fibre is joined at each end 
to the torsion arm, and passes over an accurately turned and freely pivoted pulley 
at the lower end of the rod R. The difference in weights of the loaded and unloaded 
arm is so slight that the restoring couple can be considered as unaffected by the 
differing weights of the specimens. The provision of two light paper dampers D, D, 
sufficed to damp out pendular oscillations, which at first were very troublesome. 


§7. EXPERIMENTAL PROCEDURE 


Small glass bulbs with short narrow necks are used to hold the powdered 
specimens. An empty bulb is weighed and placed in the cup at the end of the 
torsion arm; the rider is adjusted so that the arm balances horizontally, and the 
two paper caps are placed on the ends of the T-tube. The tube is then slipped up 
or down in its wooden clamp until the bulb, observed through the microscope M, 
(Fig. 4), 1s seen to occupy a position on the axis of symmetry of the field; this 
adjustment is made with ample accuracy by eye. The torsion head is now turned 
until, with the magnet energised, the bulb occupies a position at a standard distance 
from the brass plug in one pole piece. These simple adjustments ensure that the 
bulb occupies, with sufficient accuracy, a standard position on the axis. 

The pointer is now observed through the high power microscope My, and from 
observations of successive turning points of small oscillations, the positions of rest 
for no field and with the field applied are deduced. The deflections of the pointer 
never exceed 8 mm., corresponding to a movement of the bulb of 2 mm. These 
observations give the effect of the field on the empty bulb; they are repeated after 
the bulb has been filled with a powdered salt, and the effect on the salt alone is 
deduced by difference. The magnetising current is measured with a high grade 
moving-coil ammeter, and after adjustment to a standard value of 2:1 amp., the 
circuit is made and broken several times in order to obtain the true field corre- 
sponding to this current. 

Generally, the weight of the salt is of the same order as that of the bulb con- 
taining it (about 10 mg.) and the deflection in the magnetic apparatus is about 
40 times that due to the empty bulb. The weighings are conducted on a Bunge 


| 
. 
| 
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short-beam balance specially adjusted for high sensitivity in the weighing of masses 
of the order in question. 

With such small quantities of material, it is essential to carry out the work with 


the utmost care, since a very slight trace of any impurity having a susceptibility 
differing appreciably from that of the material under examination would suffice 


to invalidate the reading. This was at first a source of much trouble, but as the 


technique of the work was mastered the difficulties disappeared. In order, however, 
to detect any trouble of this nature, susceptibility measurements were made on 
three or four specimens of each substance, and the mean was only accepted if the 


‘individual determinations were in satisfactory agreement. The values of the sus- 


ceptibility deduced from observations on different specimens of the same substance 


rarely differed by more than 3 per cent. 


§8. RESULTS AND CONCLUSIONS 


In all the observations made up to the present time, pure heptahydrated nickel 


sulphate has been used as a standard of reference for relative measurements. 
Observations have been grouped, a separate calibration of the apparatus in terms 


of this salt being made for each group of four or five substances. As illustrating 
the constancy of the calibration, the results for four successive groups, extending 
over a period of eleven months, are given. The numbers tabulated represent the 
susceptibility of the sulphate in arbitrary units characteristic of the apparatus. 


October 1927 en ee 750 
November 1927 ... Ses 757; 
March 1928 Sie oe Gee 
June 1928 ... oe ie 761 


The extreme values differ by just under 3 per cent., and the constancy over the 
period of any one group was certainly of a much higher order, since precautions 
were taken to prevent disturbance of the apparatus over the period of each group, 
whereas adjustments and minor alterations were sometimes made between calibra- 
tions. 

As work on the measurements is still in progress, the publication of the full 
results, comprising measurements on a wide range of nickel salts, is left to a later 
date; but since in two cases results of another worker are available, they are given 
here for comparison. Jackson* gives the following susceptibilities : 


NiSO,.7H,O = 16:0 X 1078 at 292°2° abs. 
NiSO,(NH,).SO,.6H,O = 10°7 x 10-6 at 286°5° abs. 
The present measurements give for the susceptibilities of these salts at room 
temperatures, in arbitrary units, 750 and 517 respectively. Taking Jackson’s value 


for the sulphate as standard, this gives 11-0 x 10~® as the susceptibility of the 


double salt, which is in satisfactory agreement with Jackson’s own value. 


* Phil. Trans. A, 224, 1 (1923). 
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The consistency of the results, and the ease with which they can be obtained, 
have proved to be entirely satisfactory, justifying the further development of the 
method. Work has accordingly been commenced on the design of an improved 
apparatus, in which it is hoped to realise the theoretically desirable field distribu- 
tion more accurately. Attention is also being paid to other experimental details, . 
with the object of developing an apparatus capable of giving precision results with 
a greater ease than has been possible hitherto. 
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DISCUSSION 


Dr D. OweEN said that the author’s method represented a definite advance on ' 
previous methods. It was an important achievement to have moved in the direction — 
of making H.0H/éx constant, for variations in this quantity admittedly gave rise 
to considerable errors in Curie’s work. The speaker expressed a hope that the author’s 
method would be used in a thorough repetition of investigations such as the effect 
of temperature on magnetisation. For this purpose it would be necessary to measure — 
considerably lower values of the susceptibility than those referred to in the paper, — 
and this appeared to demand that the value of the product H.0H/éx should be 
made substantially higher than at present. Would it be feasible to design a magnetic 
circuit with a rapidly diverging field for this purpose? ; 


Dr J. H. Suaxsy: I also should like to congratulate the author on a very beautiful — 
method of using a non-uniform field to give a uniform product. Is it possible to 
adapt the method to measure the susceptibility of fluids? 


Mr J. Guixp raised the question whether the treatment of end effects in § 4 of 
the paper covered all the demagnetising effects of the ends. 


AUTHOR'S REPLY. In reply to Dr Owen, I would point out that although no mea- : 
surements have yet been made, with the apparatus described, of susceptibilities — 
smaller than 107", this is not due to any lack of sensitivity on account of field 
design; actually the value of H.cH/éx is comparatively large, and by a simple 
modification of the torsion balance it will undoubtedly be possible to attain avery 
high sensitivity. 

In attempting to increase the value of the product H.dH/éx by making the © 
field diverge more rapidly, one is faced with the difficulty that one of the poles must 
be made of small diameter; “fringing” effects then become more serious, and the 
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heoretical conditions cannot be realised so accurately, Also, an increase in 0H/0x 
y this means, for a constant magnetising current and available gap, reduces H to 
ome extent. 

I have not used the apparatus to measure the susceptibilities of fluids, but there 
is no reason why it should not be used for that purpose. 

The treatment of § 4 is intended as an investigation of the field configuration 
n the absence of the specimen. The correction to be applied for demagnetising 
effects of the ends of the latter would be of the order of only o-or per cent. and is 
therefore neglected. 


. 


264 


THE EQUATIONS OF MOTION OF A VISCOUS 
FLUID IN TENSOR NOTATION 
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ABSTRACT. The tensor notation, which is probably best known in its application in 
Einstein’s four-dimensional theory of gravitation, can be used in its three-dimensio 
form in obtaining the general equations in many branches of mathematical physics. In 
the present paper the various steps in the proof of the equations of motion of a visco 
fluid are translated into tensor form. An attempt has been made to develop the various” 
formulae in outline from first principles, with references where necessary to a standard” 
work, and also to give whenever possible a geometrical interpretation of the various vector | 
quantities met with in the course of the argument. In particular the fundamental assump= 
tions of Stokes, on which the viscous equations are based, are stated in the form of 
relation between the stress vector across an element of area and the vector element of area 
The equations of motion in tensor notation are afterwards translated into ordinary” 
notation for the special case of orthogonal curvilinear coordinates. These can be at once” 
reduced to the commonly useful forms of spherical polars, ellipsoidal coordinates, cts 
In a final section it is shown how the tensor equations can be used to deduce the equations — 
of viscous flow in terms of Stokes stream function for motion symmetrical about an axis. 
The stream function is shown to be identical with the surviving component of the co- 
variant vector potential of which the other two components vanish. ; 


§1. INTRODUCTION 

N all physical problems of a general nature and in hydrodynamics in particular 
|: considerable simplification arises if the various vector quantities are r 

presented by a single symbol instead of the three components of the vector i 
Cartesian coordinates. The result of doing so is to ensure that none but vect 
expressions occur in the equations; on the other hand, the vector expressio: 
(such as the vector product, divergence, curl, etc.) have to be manipulated b 
apparently arbitrary rules, many of which must be proved by returning to th 
Cartesian notation. The tensor notation, on the other hand, has the advantage that 
it makes use of none but the ordinary rules of algebra and the calculus, and, while 
the actual components of vectors are dealt with, it is possible by use of suitabl 
notation and by a special convention (the summation convention*) to obtain the 
same economy in the use of symbols as with the ordinary vector notation. 

In the present paper an attempt is made to develop the theory of the moti 
of a fluid in outline as far as the general equations of viscous flow by the use of th 
tensor notation in three dimensions. The actual proof of the equations of viscou 
flow is made to depend on the assumption that the resultant stress across an arbi- 
trary element of area is a function of the first derivative tensor of the velocity 
vector, is linear in this derivative and vanishes with the element of area. The 
resulting function must be a vector function, and any properties of the fluid con- 

* Appendix A. 4 
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tained in it must have the form of scalars, since the fluid can have no directional 
properties. The customary argument by transformation to the principal axes of the 
rate-of-strain ellipsoid is replaced by the assumption that the form adopted for 
the vector function representing the stress is the most general possible subject to 
the above conditions. 

Tt would have been possible to give the whole of this argument in a form re- 
stricted to Cartesian coordinates independent of the particular system, and in 
fact the form of most of the actual expressions representing physical quantities 
would have been absolutely unaltered. It seemed however worth while to use the 
more general form applicable to general curvilinear coordinates, as by this means 
it was possible to use the results to deduce the equations in general orthogonal 
coordinates in ordinary notation. 

An attempt has been made to develop the tensor calculus from first principles 
so far as it is required, though such proofs of the various steps as may be found by 
reference to standard works have for the most part been omitted. Reference is 
frequently made to The Absolute Differential Calculus by Levi-Civita (English 
translation)*. Definitions of the various physical quantities occurring in the 

‘equations of motion are given in some detail; after the proof of the equations of 
motion, the various tensor expressions, including the expressions for the stresses, 
are converted to general orthogonal coordinates in ordinary notation. 

The method of eliminating the pressures, given in Lamb’s Hydrodynamicst, is 
developed in tensor notation and expressed in terms of the vector potential. ‘The 

result is used to develop the equation for the special case of motion symmetrical 

about an axis in orthogonal coordinates, in terms of Stokes’ stream function 7%, 
by showing that % is identical with the third component of the covariant vector 
potential, the other two components being zero. 


F §2. THEORY OF TENSORS 
- ,. The essential feature of the tensor notation is that it applies equally to any 
possible system of coordinates. Starting with any Cartesian system X;,(X1,X2, X53), 
we take as the typical system of coordinates any three arbitrary functions of the X;,, 
il wy = fi; (M1, Xa, Xs). 
‘These generalised coordinates will be in general both curvilinear and oblique; but 
the tensor formulae might be restricted if required to the special case of the trans- 
formation of Cartesian axes by displacement or rotation, which would lead to 
_ important simplifications in detail but would leave the general results unaltered. 

In the tensor calculus we are concerned with transformations from any co- 
ordinate system x; to any other system *,, where either system is considered as 
defining the position of an arbitrary point P; the x, are considered either as func- 
tions of the x; or of the original Cartesian system. 

The simplest case is that of a function of position ¢ which has a definite value 
at every point P. It may be considered as a function of the coordinates x, and is 
obviously transformed by the relation of ‘‘invariance ” 

b (X41, Xe, Xs) = (ay , Xp 5 X3)- 
* Blackie and Son, 1927. + 4th edition, § 328, equation (8). 

PHYS.SOC.XLII,3 18 


5 6e C. N. H. Lock 


Before proceeding to discuss the possible ways of expressing a vector function 
of position it is desirable to remark that the tensor calculus is only concerned with 
“infinitesimal” vectors. Every vector used in physical problems can in fact be 
represented by a directed line of infinitesimal length, with the exception of the 
straight line representing the position of the point P; the position of P is not 
considered as a vector in the tensor calculus. 

We start therefore by considering the infinitesimal vector A, which is a function 
of the position of the point P, as being represented by the line joining P (coor- 
dinates x,) to a neighbouring point O (coordinates x; + dx')*. Consider also the 
points R,, R,, R, whose coordinates are (x, + dx, X2, x,) etc. and T,, T;, T; whose 


coordinates are (x, %, + dx,, 3 + dx,) etc. These eight points form to the first — 


order a parallelepiped which is in general oblique angled. 


One method of representing the vector A is by means of the three components 


ES) A 


ied 


dx' (dx, , dx, , dx3) which are proportional to the lengths PR, , PR,, PR,. An equally — 
reasonable method would be to make use of quantities proportional to the projec- — 


tions PN, PN», PN, of PQ on the axial lines PR, , PR,, PRs. 
Write ds for the length of PQ and consider the formula connecting ds with the 


dx’, ds may be considered as an invariant function of A and of the position of P. 
In Cartesian coordinates 


dst =dXj°+dX24+dX% 4 2 2 2 2 cua (2-1), 
which may be written dX,dX™, 
with dX, = dX, 
By direct transformation we have 
op Gg OA" hoe 
ds? = joahe ee dx'dx'*, 


*® See Appendix A. 
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or writing : 2: for are Ay tea (222), 
t 7 
Cenk Cn ee an ae eee ORE ERIE. er (222), 


with the relation. £i7 = Z;- By definition g;; are functions of position which depend 
only on the particular system of coordinates, not on the particular vector A. 
From the geometry of the parallelepiped, writing 0,3 for the angle R, PRs, etc., 


CE ge 2ail Rabks COSUsne ees (2°4). 
Comparing equations (2:3) and (2:4), we obtain the equations 
: Tpke Lae nn ee (2°5), 
etc., and cos 05 = —2 4 VAR bo Se ee (2-6), 
§22" §33" 


etc. Equations (2-5) and (2-6) connect the g,; with the geometry of the parallele- 
piped. 
We are now in a position to define two types of components of the vector A. 
These will be described as the “‘contravariant” components A* and the “co- 
variant” components A,. The A“ are defined at once by the three equations 
ee eee ee) ie Sate ns (2-7) 
so that equation (2°5) gives PieP Ri sea, CV es (28). 
Again PN,, the projection of PQ on PR,, is equal to the sum of the projections 
of PR,, PR, PR, on PR, so that 
PN, = PR, + PR,.cos 04. + PR.cos s 


1 


1 g. 1 
= £1" dx, + — Se i » £097 AXy + — ou 1 tea ane 


£11" §22 £1" 833" 

I 
= 4 S1adea™ ere 3 (2:9). 

; &u 

We shall define the covariant component A, of A by the equation 
a Pa ae) ee eT ea oo (2°10). 
It follows from equation (2-9) that 

eee ING eo oo top ess (2:11), 


etc., which give the geometrical definition of A, in a form similar to equation (2:8) 
for A*. Alternative definitions of A* and A, may be given with reference to the 
parallelepiped consisting of planes perpendicular to PR,, PR,, PR;. It appears 
that A* and A, are equally well qualified to be considered as components of the 
vector A. 

When the coordinates are orthogonal the point N, coincides with R,, and PN, 
is the “component” of A as usually defined. Also, when i +7, g;; is zero by equation 
(2:6). It follows from equations (2-8) and (2°11) that 

; gu? Al = PN, = =, Ae a (2°12). 
11 
In Cartesian coordinates g,, = 1 and the distinction between A‘ and A, disappears. 


* See Appendix A. 
18-2 
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§3. SOME FUNDAMENTAL ALGEBRAICAL RESULTS 


At this stage it is convenient to define the most general type of tensor and to 
obtain its laws of transformation, showing that A’, A; as defined above are particular 
cases. We take dx! as the typical contravariant tensor, which evidently satisfies the 


law of transformation ax 
dxt = dx® eae (3:1). 


Ox* 


An alternative definition of a covariant vector U;, is given by the equation 
b=Uzdxee = nana (3-2). 
It may be shown* that the values of U; in any other system of coordinates are 


consistently and uniquely defined by the condition that ¢ is invariant for any 
arbitrary values of dx‘, and that the U; so defined satisfy the law of transformation 


U; > Ox, Uqt lek Lt (3 3). 
Both the laws of transformation (3-1) and (3-3) may be shown to be transitive and 
also invertible, so that e.g. ox. 

U, == ax U; pk trae (3 4) 


The covariant vector U; has still to be identified with the covariant component 4,, 
defined in equation (2°10). 

In addition to the contravariant and covariant vectors having three components 
it is possible to define a more general type of tensor having 9, 27, etc. components 
by an extension of equation (3-2). For example write 

d =— Aw U; U,’ U,’dx* dx’® 
so that A,,“* has 35 components. It may be shown that components A,,%* in any 
other system of coordinates are consistently and uniquely defined by the condition 
that ¢ is invariant for arbitrary values of the U; and dx*, and that the 4,,"* satisfy 
the law of transformation 

A, OX; OX; OX, Oxy OX, A," 


~ Oxy OX, OXy OFy Oy 29 TT (3°6). 


The Aqy‘* are described as the components of a tensor of the fifth rank having 
three elements of contravariant and two of covariant character. 

As a particular example, a comparison of equation (2-3) with equation (3-5) 
above, shows that the g,; are components of a covariant tensor of the second 
rank; it is of fundamental importance, being a function of the coordinate system 
only. 


It may be shown, as a particular case of (3-5), that the “ outer product” of any 
number of vectors or tensors is itself a tensor, e.g. 
A‘. By = C8 
is a tensor of the kind indicated. 


* Levi-Civita, The Absolute Differential Calculus, chap. 4, § 11. 
t Asimple method of remembering this law of transformation is afforded by equation (6:1) below. 
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Contraction. It may be shown* that the result of putting a lower suffix equal 
to an upper suffix in any tensor and performing the indicated summationf, is itself 
a tensor. Thus A,,** is a tensor obtained from the tensor A,;“7 by the process of 
“contraction.” In particular the result of contracting the ‘‘mixed” tensor A,’ is 

“the scalar quantity 4,". 

A further extension of equation (3-5) is that if (e.g.) B,”’ is an arbitrary tensor 
then the condition that any product (contracted or not) of B,” into Ayn,’ is a tensor 
of the kind indicated, defines A,,,,/ uniquely and consistently in any coordinate 
system such that it is a tensor of the kind indicated. 

A very specialised type of mixed tensor is defined as follows: 


di =o when 7+)’ 
Sou ues) Ah ee keroee (38). 
The fact that 5,’ satisfies the correct law of transformation may be verified by direct 


substitution. By virtue of the above rule of contraction it may be verified that the 
contracted product of 8,‘ into e.g. a contravariant vector A®* is given by 
OA A) is Be Pas: (2°90); 
so that 5,* has the property of substituting one suffix letter for another in any tensor. 
Other important properties which may be easily verified are embodied in the 
equation ant 
: nox 
Of = 2? 241 = ae, ee (@sieyy 
Application to geometrical vectors. We may now return to the subject of the 
representation of geometrical vectors by means of the tensor notation. It follows 
from a remark in the section on contraction, that the covariant component as 
defined in equation (2°10) is equivalent to the covariant vector as defined in equa- 
tion (3:2) and satisfies the same law of transformation, equation (3-3). The three 
i i tions 
simultaneous linear equati ey As 


may be solved for the A® in the form 
At=giA; teers (3712); 


where, according to the ordinary rule for solving simultaneous equations, g” is 
the cofactor of g;,; in the determinant having g;, as elements divided by the value g of 
the determinant, and it follows from (3:12) that g is a tensor. Thus the funda- 
mental tensors g,; and g‘/ can be used to convert the contravariant components of 
4 vector into the covariant components and vice versa. 

We are now in a position to show how either the A‘ or A; are convenient for 
representing the vector A in any system of coordinates. The first observation 1s 
that the condition A; = 0 involves the vanishing of A. Again, equations of the type 


A, =8,, 


“hold in any system of coordinates, and involve the identity of the vectors A and B. 


* Levi-Civita, loc. cit. chap. 4, § 9- + Appendix A. } gi is defined in the next section. 
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In view of the linearity of the equations of transformation, it follows that the 


equation Kk, 2A wee Fee (3:13); 
holds in any system of coordinates, and it may be shown geometrically that the 
vector K is equal to the sum of the vectors A, B, C according to the parallelogram 
law. 
§4. GEOMETRICAL EXAMPLES OF TENSORS OF THE SECOND ORDER 
The “outer” product of two vectors. The results of § 3 show that the “outer ” pro- 
duct of the components of two vectors Aand B, e.g. A‘B’ is a tensor of the second rank 
having nine components. Since there are two types of component of each vector, we 
can construct four tensors of the second rank, A‘ B’, A‘ B;, A; B’, A; B; (one contra- 
variant, two mixed, and one covariant), which may be considered as being ‘“‘com- 
ponents” of a tensor (generalised vector) of the second rank. 
Linear vector function of a vector. We consider a vector F which is a function 
of the arbitrary vector A. F (A) is described as a linear function if it satisfies the 
following conditions: (i) if A is fixed in direction, then F is fixed in direction, and 
of magnitude proportional to the magnitude of A; (ii) if A(1), A(2) are any two 
values of A, and A (3) is their vector sum, then F {A (3)} is the vector sum of 
F {A (1)} and F {A (2)}. Under these conditions it may be shown that the com- 
ponents of F are equal to the reduced product of the components of A into a tensor 
T* of the second rank, e.g. P=74,. 7° ~ eee (4:1). 


Proof. Let A (1), A (2), A (3) be three values of A having different directions. 
Substitute these and the corresponding values of F in equation (4-1) and we obtain 
nine simultaneous equations which in virtue of the results of § 3 determine the com- 


ponents of 7 as a tensor in any system of coordinates. Now any value A (4) of A 
may be expressed in the form 


A(4)=m%A(1)+a,A(2)+a,A(Q3) ss (4-2), 
where a, @, @ are scalar constants and it follows from the conditions satisfied by 
F that 

ss F (4) = a, F(1) +a, F(2)+a,F(3) >. sna (4:3), 


and therefore F (4) and A (4) satisfy equation (4:1). Thus equation (4:1) represents 
the necessary and sufficient condition that F is a linear vector function of A. The 


tensor T is of the most general type, unlike the example of the last section (product 
of two vectors). 


§5. CONTRACTION OF TENSORS OF THE SECOND ORDER 


From any tensor 7* of the second rank we can obtain by processes of contraction 


a scalar and a vector which may be described as the “scalar contraction” and the 
“vector contraction” of T respectively. 


Scalar contraction. 'The contracted mixed component 7%* of the tensor T is a 
scalar: the scalar contraction of 7’. 


* By analogy with the outer product we consider any tensor 7, of the second rank, to have the 
four components: T”, T;, T/, Ty. Any three of these may be expressed in terms of the fourth by 


relations identical with those obtainable byapplication of equations (3°11) or (3-12) to the outer product. 
It follows that 7% =T7,", and this contraction will therefore be written TY. 
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Vector contraction. Epsilon tensors. It appears on examination that the anti- 
symmetrical tensor 7“ — T°“ has only three distinct components. It may be reduced 
to a vector by means of the fundamental auxiliary tensor of the third rank F,,.* (in 
space of three dimensions) defined by the conditions 


Bie = ©} 
unless a, 6 and ¢ are all different, when 
1a =o gt ehalstelete (Gage 


the sign being that of the corresponding term of the fundamental three-row deter- 
-minant whose value is g. The corresponding contravariant tensor is defined by 
 Giik pie gibgke Hy tte (52), 
and may be shown to satisfy similar conditions to those of equation (1) but with g 
replaced by g?. An important result which may be easily verified is that 
. eee o, Oh Oe ee atts: (53): 
Now consider the vector L defined by the alternative equations 
Teal eale el) 
SMe Aa Ge oe) ag Gree: (5:4), 
: Li =4E@ (Ta — Ty \eeee Lee oe % (5-5). 
which may be shown to be all equivalent. Multiplying both sides of (5-4) by 
Ete and using (5-3) we have 
Bury, = (8/54 — 818%) T 
SU Led Waly eae oe he Ce (5:6), 
by equation (3-9). By means of the equations (5-5) and (5:6), the vector L (the 
“vector contraction” of 7) may be expressed in terms of the tensor (Tit — T*) 
and wice versa. 

Some additional identities involving epsilon tensors which will be required 
below are given in Appendix B. 

Scalar and vector products. Examples of these two types of contracted tensor 
are given by substituting for the tensor 7 the “outer product” A“ B? of the two 
vectors A and B. These are the scalar product 4 and the vector product V of the 
two vectors defined by the equations 

ie SORTED 0 I te Sek (5°7)s 
and [a Bi ters (58), 
similar to equation (5-4). It may be shown that these are identical with the scalar 
and vector products of two vectors as ordinarily defined. 

If we take the scalar product of V into a third vector Cit is obvious geometrically 
that the result is equal to the volume of the parallelepiped having A, B, C as 
adjacent edges. Alternative expressions for this volume are therefore 

CN ODAC Ey Ai i a ah aire (5°9); 
A, B,C, 
A, B,C3|* 
A; B, C; 


* Levi-Civita, loc. cit. chap. 6, § 8. He uses € in place of E. 


I 
g 
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§6. DIFFERENTIATION OF SCALAR AND VECTOR 
FUNCTIONS OF POSITION 

Derivative of a scalar. The change 6¢ in a scalar function ¢ due to a small 
displacement from the point P to the point Q considered as a vector PO=Bis 
obviously given by 4 

id Se se dx pe (6-1). 
It follows from equation (3-2), that since 6¢ is invariant, 0¢/éx, is the covariant 
component of a vector. This vector is the gradient of the scalar ¢ as usually 
defined. 

Derivative of a vector. Example of a linear vector function. For an arbitrary 
vector function of position A we have to consider the effect on A of a small dis- 
placement from the point P to the point Q, where PQ is considered as an arbitrary 
vector B whose contravariant components are dx’. Let 5A be a vector equal to the 
difference between A, and Ap or more exactly to the difference between Ay and 
the vector through O equal and parallel to Ap. Obviously 5A is a function of the 
vector B and it may be shown geometrically that it is a “linear vector function.” 
Hence it is possible to define the four components A‘:*; A‘ ,; A,»*; A; of the 
second rank “‘tensor derivative’’* of A by equations such as 

sta At Be Alden ee (ay 
In general when the coordinates are curvilinear A‘, is not equal to @4*/@x* nor 
is 0A‘/ox* a tensor. The following expressions for A‘, and A, , are obtained in 
Appendix C, by direct transformation from Cartesian coordinates: 


t = oA ; t / 
A C= ex, Ae De , ee (6 3)> 
0A; 
Hiern =. lJAs. - | =e ee (6-4), 
i — ate (O80 OL%a 0g; 
where I'yn.° = 4d Ox, ae, ox, ee (6-5), 


is known as a Christoffel symbol of the second kind. h 

Derivative of a tensor, Similar expressions may be obtained for the derivative 
of a tensor of any rank, e.g. the tensor 4, has the covariant derivative 

0A b - ~ 
Aw, aes ax, = Poe’ Aai Pact Ag 

It may be shown that the covariant derivatives of the fundamental tensors 83, 4, 
Eq, are all zero so that these tensors behave as constants for the purpose of co- 
variant differentiation. It follows that A, , — Lin A® q. 
Divergence of a vector. The “scalar contraction” of the tensor derivative of A 
is A* , which is identical with the divergence as ordinarily defined. It may be 


* 'The “comma”? notation for the tensor derivative is used throughout the paper. Of the four 
components of the tensor derivative of a vector (see p. 270, footnote), A‘, and A; , are defined by 
equations (6:3) and (6:4), and satisfy the relation 4; .=g,,.4° a While Aba d Ape 
the relations A’ * = 9” At, and Ap*=g® A, y, Pury axis 
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erified that A* , may be expressed in the following form which does not contain 
Christoffel symbols: 4 
div. A= A*,=g4 eo Eri Olek tae el ae eee (6-7). 


Curl of a vector. The “vector contraction” of the tensor derivative of A (with 
ign reversed) is identical with the curl of A as ordinarily defined and may be 
ritten “curl A.” From equation (5:5) we have 


(curl Ay eA, 
See (Ang — Ags) a «a8 ©) cesses (6:8). 


hen the second form is written out in full the Christoffel symbols cancel and we 
ave therefore 


=, 0A 
cure) eo eee 
(curl A) = Bi (6-9), 
form independent of Christoffel symbols. Since the operation of obtaining the 
url changes a vector into a vector it may be repeated any number of times. For 
is purpose it is convenient to obtain formulae connecting corresponding com- 


onents of A and curl A; these are: 


(curl A); = 21, L°"° a AO se (O7rS), 
(curled y= bie se (Get) a eee: (6-12). 


Second derivative. By substituting the first derivative of a vector (equation 6-4) 
n the formula for the first derivative of the corresponding tensor (equation 6-6) we 
btain the second derivative of the vector A of which e.g. the covariant component 
ay be written A, ,,,. It is important to notice that in Euclidean space with which 
lone we are concerned, the order of covariant differentiation is indifferent; this 
act involves a series of complicated identities among the g, and their derivatives. 

Contracted second derivative V?. The contracted second derivative (contra- 
ariant component, A’:* ,) is the component of a vector which is identical with 
2A. This vector may be reduced to a form independent of Christoffel symbols 
by using the identity usually written 


V?A = erad.div.A—curleurlA 2 ws (G12); 


since the divergence, curl, and gradient have all been simplified in this way. A 
proof of this identity by tensor notation is given in Appendix D and an alterna- 
tive method of simplifying V2A in Appendix E. 

Reduced second derivative of a scalar (sey “f. This is identical with the divergence 
of the gradient of ¢ and is also usually written V2¢, but it is important to notice 
that the formulae for V2¢ and V2A only agree in the special case of Cartesian 
coordinates. 


* Levi-Civita, loc. cit. chap. 6, § 7. 
+ See footnote, p. 272. 


274 CON, Hi: look 


§7. VECTORS AND TENSORS REQUIRED IN HYDRODYNAMICS 


Velocity. If a particle of fluid at a point P is displaced to the point O ina sm 
time interval 5¢, then PQ is an infinitesimal vector and the velocity u may b 
defined as the vector 


I 
u= PQ tise Si (7-1); 


the time 8¢ being considered as an independent variable unaffected by change 
coordinates. 

Acceleration. We shall adopt the Eulerian convention by which the velocity u 
at a fixed point P in space is considered as a function of the time and of the co- 
ordinates of P, and shall deal only with the case of fixed axes. In determining t 
acceleration of a given particle of fluid write up (t) and up(¢ + 52) for the velociti 
at the point P at the beginning and end of the small interval 5¢. At the end of 
interval the particle of fluid will have moved to the point O where PO = updt 
the first order and its velocity will then be ug (¢ + 5t). Let f be the vector accele 
tion of the particle, then 

£5¢ = Ug (¢ + dt) — Up (2) 
= Ug (¢)— up(#) + Up(t+ d4)—up() «we (7-2), 
to the first order; in tensor notation (contravariant components) this becomes: __ 
: 2 = TGs 
fst = Sut + Pi i tie. (7°3)- 
By the results of section 6, du‘ the contravariant component of the differen 
between the simultaneous velocities at O and P is given by 


Out = at de® | te a (7-4), 
where dx* the contravariant component of the vector PQ is given by 
x®= u*Ot 80 i RW ee ae (775). 9 
; i i 4a Cut 
Hence we have fi=ut,ut+ = | | See (7:6), 
ae en 
and similarly fi = Uj, gut + ho . @= 26a (7-7). 


A similar formula applied to any vector function of position corresponds to the 
usual formula for differentiation following the motion of the fluid*. i 
When u, ,v* is written out in full it is found that two of the three terms in 

Christoffel symbols cancel, leading to the result: 


Ou; og 
, aie b 
My gi® = 09 ee eee (7°8). 
es Ox; ‘ 


Rate of strain. It is usual in works on hydrodynamics to separate the derivative 
of the velocity w, ; into two components A,; and B,; given by 


Ay =} (ug a) a ee (7-9), : 
B,; = 4 (u; ; mat oe) ee Se ee (7:10), ° 
So that Us; = A;; + Bi; weeee (7-11) » 


* Lamb, Hydrodynamics, § 5. 


a 
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appears from equation (6:8) that the tensor B,; can be expressed in terms of a 

ctor € (the curl of the velocity) given by the equation 
Ete (CULL) os Ltt, go tal Mid Gapste. (elec 

his vector is the ‘‘vorticity” as ordinarily defined; the tensor A;; is described as 
e ‘“‘rate of strain.” The proof given below shows that the stress components are 
nctions of the rate of strain and independent of the rotation so that it becomes 
necessary to make the separation at this stage. The tensor B or the vector € can 
expressed in a form independent of Christoffel symbols by means of equations 
10) or (6-11). 3 

For the tensor A it will be found that when A” is written out in full some of 
e derivatives of the g, in the Christoffel symbols cancel, giving the result 


‘ 


to Pal: 
_ Out _ ous eae oe 
a ee L gia L 4% gid gic ses em os Tie 
OXg T § OXq Ss & OX g (7 3) 
§8. STRESS 


~ Consider an arbitrary small element of area dS considered as a vector defined 
being perpendicular to the plane of the element of area, and equal to it in magni- 
de, and let dF be the resultant force of reaction across this area in the fluid. dF is 
function of the vector dS and it will be shown that it is a linear vector function as 
efined in § 4. Consider the equilibrium of the element of fluid bounded by any 

all tetrahedron. The forces dF on the four sides have zero resultant to the first 
rder since the inertia and external forces are of a higher order of small quantities ; 
e vector elements of area of the four sides have zero resultant by the geometry 
f the tetrahedron. Again, when dS is fixed in direction dF is fixed in direction and 
roportional to dS. Hence dF is a linear vector function of dS, and by § 4 a tensor 
of the second rank may be defined by an equation such as 

Gi PUGS, = cos (8-1). 

he stress tensor P then determines dF for all values of dS. 
+ The scalar contraction of P is P* which may be written 
dF, ¥: dF, dF, dF'. dF? | di® (8:2). 


seeeee 


Fi— 4S, dSe dS, Si) dS*” dS* 
t is usual to write — 1P“ = p the ‘‘mean pressure” in the fluid. 
Equation of moments. Symmetry of P”. An additional relation may be obtained 
rom the equation of moments on the fluid contained in an elementary parallelepiped 
hose adjacent edges are the vectors d& .dn.dt. If dF (€) is the stress across the 
ace bounded by dn, dé, the sum of the moments of the forces on the opposite 
‘aces about the centre of the parallelepiped is the vector M (€) equal to the vector 
roduct of dF (é) by d&, so that 

; Mila Lied (6) abe ee (8°3). 
gain the area dS(£) of the face of the parallelepiped bounded by dy and d€ is given by 

Sn ns (8-4). 
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Adding the moments of the forces on the other pairs of opposite faces, we have fa 


the whole moment M, = M, (é) + M; (n) + M, (0) 
e EjjqP*dS;,d&* + sca 9 Sa ee (8-5), 
= {Eijg Ene + Esp nce + Esse Eran} Pikd&* dn? dl° 
=1E 5 (P* — P™) Eddy al" oe ee (8-6), 


using the result of Appendix B, equation (B-2). Since all other forces are of a high . 
order of small quantities, it follows that M vanishes and hence from equation (8-6) 
pt pa eee (8-7). 
Equation of motion in terms of stress. It is now possible to write down 
equations of motion of the fluid in terms of the stress tensor, by application of 
d’Alembert’s principle to the equilibrium of the fluid contained within the elemen- 
tary parallelepiped specified in the last paragraph. The difference ddF (£) betwee: 
the stresses on the opposite faces bounded by dn, dg is given by 
SdF (€)' = dF (©) dé" 
— P# dS (6),dé 
== Pi! oBg désdqtdl® = saan (8-8). 
Again, the mass acceleration of the element of fluid is 
pf’ Egredé* dn? dL? 
using the expression for the element of volume given in equation (59). Equati ng 
this to the resultant stress (assuming no external forces), we have the equations 
pf'Eqyedé* dy? dle — 3dF (€)' + 34F (n)' + 84F (Q 
a eae os & + En Bice ie PS Bay d&* dn? dt 


= Pi gE ay detdyr dee aaa (8-9), 
using the result of equation (B-1). Hence 
Si = ==3 pf 
=pie au et . 5 Oe 8-10), 
P | ct +u +a u j ( I )> ; 


which is the required equation of motion. 


§9. PROOF OF THE EQUATIONS OF VISCOUS FLOW 

In order to obtain the final equations of motion in terms of the velocities’ 
remains to formulate a relation between the stress components and the velocities, 
Consider the vector dF representing the stress across an arbitrary element of ares 
dS in the fluid. The vector dF must satisfy the following conditions: (i) It must be 
a linear vector function of dS. (ii) It must be a linear function of the first tenso: 
derivative of the velocity vector u. This corresponds to the ordinary assumption 
on which the theory of viscous flow is based. (iii) The function must contain no 
additional variables other than scalar functions of the density, temperature, ete 
since the properties of the fluid are non-directional. 
There is no need at this stage to make use of the additional condition that the 
results must be independent of the vorticity; it will appear that this follows auto- 


matically from the condition proved in equation (8-7) that the stress tensor P is 
symmetrical. | 
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The vector dF must therefore be given by an equation of the form 
dW-=adA + BdB-ydC-in. 981) Whee (9:1), 
here dA, dB, dC, etc. are vector functions of dS and the tensor derivative of u 
hich satisfy the above conditions, and a, f, y, etc. are scalars. 
It appears on examination that there are only four independent vectors satis- 
ing the above conditions. These are 
a em a (O2)e 
dA =u" dS, 
dB = the change of uin a displacement PQ equal to the vector dS 
(i.e. perpendicular to dS and equal to'it). Hence from 
equation (6-2) 


ets the Be 0, 54 4, SBC aria (or )3 
dC = the gradient of the scalar product of uand d$ treating dS as 
a constant. 
NO es (dS,u), (dS, constant), 

URS) i mers es (9:4). 

mare we may Write as — ad St + BdS*u', + ydS,us* 
Se Ba ye dog — 8 eas (9°5); 
here eg GU ge aa Foams (9°6), 


d a&,@,8,y are scalar functions of the temperature, density etc. only. It is 
teresting to notice that the last term will contribute nothing to the equations of 
notion for an incompressible fluid. Comparing with equation (8:1), it follows, 
ince d§ is arbitrary, that 
Pi = ag + Bubi + yuhé 
SF ale els Oscai ) d rrrt (9-7), 

being set equal to f since P“ has been proved to be symmetrical. This proves that 
i involves the pure rate of strain only and is independent of the vorticity. 

We might equally well start with equation (9°7), basing it on the assumption 
hat P# is a linear function of the pure rate of strain. The above development, 
eing based on equations between vectors, appears, however, to be susceptible of 
_more definite physical interpretation. Thus for an arbitrary element of area, the 
rst term adS represents a normal pressure on the element; the second term fdB 
epresents a force parallel and proportional to the rate of change of velocity in a 
‘rection normal to dS; the third term ydC is required to maintain the rotational 
uilibrium of a volume element of fluid. 
~ We can now write down the equations of motion by substituting the values of 
4 from equation (9-7) in equation (8-10). ‘The resulting equation is 

pf’ = Pi4, 


Oe ‘s. de 
=g Fama ees ‘a 


= git 2 fat Bubs} + Bi ee oe, 
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It is here assumed that f is independent of position, and in reducing the last te 
use is made of the condition that the order of differentiation is indifferent. 
It is usual to express @ in terms of the mean pressure p (= — 4P%) defined 

equation (8-2) above. Equation (9-7) may be written 

Pi = of +B tpt ue). ae ees (9:9). 
Contracting by putting 7 = 7 gives 

— 3p = Pi = 30+ 2Bu a 
or a=-—p—$pu, | © "se eet (9°10). 
Substitution in (9°8) (f* being written in full and f identified with j, the coeffici 
of viscosity) gives* 


a 


ou’ ret. ; ; P 
p - + ut gu aa soe 3S {— p+ Zuur a} + pu? g (9°11), 


as the final form of the equations of motion. 
The corresponding expression for P” is 
Pi = — gi(pt+ Zput)tp(ureita) —...... (9-12). 
It is important to notice that in virtue of equations (g-10) and (9-6), p may 
expressed in the form ; 
* —p=&+ (q+ Fr) U0 


where a, a, are functions of density and temperature onlyt. 

The equations of motion for covariant components and the expression f 
mixed and covariant components of the stress tensor can be written down at on 
from equations (9:11) and (9-12). 


§10. REDUCTION TO GENERAL ORTHOGONAL COORDINATES 

It has been shown that the condition that the coordinate system is orthogo 

is that g,; is zero whenz +7. It is convenient to write 
ds* = a* dx* + b® dy? + c*ds* 

so that iu=@,ete, gabe, gierw/@ete 4 (10-2). 
We shall also make use of the ordinary components of velocity, which will 
denoted by U, V, W. It was shown in equation (2-12) that they satisfy the relati 
yi = Vguu = U/a,etc. . == ” eat (10-3), 

and ul, = Vg,,U = aU, etc. ~ al (10°4). 
Again if dS (1), dS (2), dS (3) are actual elements of area of the three coordina 
surfaces, the three ordinary components of stress across dS (1) are written in 
ordinary notation py, dS (1), Py dS (1), pis dS (1). In orthogonal coordinates 


covariant components of dS (1) are {adS (1), 0, o} and the contravariant componen 
of the stress are: 


(t/a) PudS (1), (1/6) Prd (1), (1/C)pis@S (1)} 
by equation (2°12). Hence equation (8-1) gives: 
(1/a) py dS (1) = aP'dS (1) |, 
(1/d) PydS (1) aP"dS (1)| Sai (105%) 
* 'The equation was given in this form by Milne, in Proc. Camb. Phil. Soc. 20 : 
t See Lamb, Hydrodynamics, 4th edition, § 358, equation (1), where pis ny ae a 3 
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c., so that the nine components of stress as usually defined are 

GEE er ab Pte ac Pl! 

Cg aE 02 P53 8 OCP ee ae tn He (10:6). 

Ger CCl * C238 
_ All the expressions occurring in the equations of motion and in the expressions 
r the stresses have been reduced to forms which do not involve Christoffel 
mmbols* ; it is only necessary to substitute the above values to obtain the equations 
motion in orthogonal curvilinear coordinates in full. 
_ For the equations of motion in the ordinary components it is sufficient to 
ultiply the first contravariant equation by a or divide the covariant equation by a; 
r the stresses, the above expression in terms of the contravariant components or 
e corresponding expressions in terms of the mixed or covariant components 
ay be used. It is sufficient to work out a single example in full and we shall 
joose the ‘“‘curl” which is required in reducing u’* , by using the formula 


: (870) “(curl curl aye 9.2... (10°7). 


Ut = gia 
aS OX 


Ve have from equation (6-11) 
 G 
(curl a)? = £4" 256 (x1 U’). 
ence for the first component in orthogonal coordinates 
tthe 123 ret 3 132 ay 2 
(curl uj! = E ax, 60" +E 2, S20 


ince the summation for the g, disappears) or in the notation of this section, 
First component of curl (U, V, W) = a (curl u)" 
a fr WW pone 2 it 


abe dy c oz ~b 
The first of the three equations of viscous flow in general orthogonal coordinates 
ay now be obtained by multiplying the first contravariant equation by a; the 
esulting equation is: 
mero AY WS LV. 0 | Wo U20a V?20b W?dc 
Oke ae ax b By) se) 4 ~— @ dx abdx ac 2a} 
) BOB ge ll a, ab) 
aox 3a 0x abc \ 0x oy Oz 
(dc (obV daU\ 0 Bb ( 7) (16-8) 
be . ab ( Ox = apace 0s). Ox Jy 8 
The typical equations for the stresses may similarly be derived from equations 
7-13) and 9:12) in the form 
au 1 (obcU dcaV , dabW\ | 2p fet i a fee 
e P fo ae Bye a) on ' b dy € Os 


Deg =} | = 4 fe : 55 (=) ceeish (10°10). 


yi? + wht in equation (7°13); u*%,4 in equations 


* 1; ,u*in equation (78); v*,« in equation (6°7); 
6:12), (6°7) and (6:11), or Appendix E. 
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The equations for any special case of orthogonal coordinates such as spheric 
polar coordinates may be obtained at once from the expression for ds? 


ds? = dr® + r2d6? + r? sin? Odd? 
as compared with ds? = a®dx® + b? dy? + c?dz?*, 
and it is only necessary to substitute 
x=r, y=9, z=¢, a=1, b=s7, c=rsin8, U=U,, V=U, W=U, 


to get the first equation, and the proper permutations of a, 4, ¢, U, V, W to ge 
the second and third equations. 


§11. ELIMINATION OF THE PRESSURES 


By the taking of the curl of both sides of the equations of motion (equatiot 
g'11) the pressure p is eliminated by virtue of the obvious result that the curl of 
the gradient of any scalar is zero. The resulting contravariant equation is 


© (curl us)? + {curl (wi ,u’)}* = v (curl u’*,)* —...... (11I-t), 


+. we eS we Gas 4A 


ie : 
wit, = gi? = (u*,)—(eurleurlw)i a. (11-2), 


the last term reduces to 
— vy (curl curl curl x)‘. 


The details of the reduction of the second term are given in Appendix 
Writing € for curl uw, as in equation (7-12), we can show that 


; og* Out. ; 1 
{curl (w?,,u?)}* = u2 = S ra of (ut ght nar Satan (11-3). 
The equation of motion therefore reduces to 


2: 28s OF aaa | 7 
ape + i. Se et (u*.) € = — v (curl curl g)'  ...(11-4)f, 


with — = curl u. 


For an incompressible fluid it is possible to define a “ vector potential” ht by h 
relations 


u=curlh 
Re ge eee (11°5), 
so that §=curl.cun is”) Weg ee ee (11°6), 


' and the last term on the left hand side of equation (11-4) vanishes. 


* : : 

The notation (curl 4‘)‘ stands for a contravariant component of the curl of the vector whose 
contravariant components are A‘, and is identical with (curl 4)‘, The former notation is necessary 
when 4! is an expression such as ué »u°. 


t+ Lamb, Hydrodynamics, 4th edition, § 328, equation (8). } Ibid. § 148. 
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§12. MOTION SYMMETRICAL ABOUT AN AXIS 


The last result leads directly to Stokes’ form for the equation of viscous flow 
hen the motion is symmetrical about an axis. We take for the third coordinate x,, 
he angle of azimuth of the point P with reference to the axis of symmetry, and for 
1, X_ any two-dimensional system of coordinates in planes through the axis. It 
ollows that 

ds? = 91, dx", + 2812d% dX, + Yodxq” + wddx”, 
nd that all the g, are independent of x;, @ being the normal distance of P from 
he axis. ; 

We now assume for the components of the vector potential h that h, = h, = 0; 
hen it may be verified that vu, = o and that ,, wu, are functions of hs which auto- 
atically satisfy the equation of continuity for an incompressible fluid in virtue 
fthe identity 

div curl A=o 
nd therefore the above assumption with regard to h gives the most general motion 
f an incompressible fluid of the kind considered. It may also be shown for 


(curl curl h) = € 
hat & Fy Es == ©) 
nd similarly that (curl curl €), = (curl curl €), = 0, 


ll variables being independent of x,. On examination it appears that the equations 
f motion in the form (11-4) reduce to a single equation in the variable hy. 

In order to identify h, with Stokes’ stream function #, we proceed to obtain 
n expression for the flow across any surface of revolution bounded by a circle 
bout the axis. By Stokes’ theorem and the definition of the vector potential h for 
n incompressible fluid, this is equal to the line integral of h taken round the curve, 
.e. to fh, dx” taken round the curve. In the present case this reduces to 


fs 


2n 
hs aks — 2th, 
/O 


hich agrees with the definition of Stokes’ stream function */*. 

It is possible to reduce the equations of motion in the form (11-4) to a single 
ensor equation in the variables x,, %»:/, (= #), @ and &, being treated as scalars. 
It is sufficient for the present purpose to consider the special case in which the 
o-dimensional coordinates are orthogonal. 

Orthogonal coordinates symmetrical round an axis. Using the notation of § 10, 
utting c = a, and using equation (6-10), we can verify that 


0 

u, = (curl h), = < 5 Stas: (1251). 
bo 

U, = (curl h), = — a ae Tele Saat 5 (12:2), 


Ug = O, 


* Lamb, Hydrodynamics, § 94. 


I 
PHYS. SOC, XLII, 3 9 
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riting yb in place of h,. Hence 
ale ee et 12°3), 
mb oy 
tei eie 2 Custos (12-4). 
wa Ox 
Further, we have &=&=9, 
€, = (curl curl A); 
a (0 b ob 0 a op) kbs, 
ab (oe en +o eae 
= =D) (¥), ’ | 
say. It is of interest to notice how this expression differs from the expression 
V2¢ in two dimensions, where ¢ is a scalar; this takes the form 


866 2 2 ee (12-6). 


rf 
ab \6x adx Cy bey 


V26 = 
The equations of motion (11-4) may now be written down in the form 
Of , w (op @ (&)-2 e (&)t u2—pe hae i 
Gt ab ley dx \w?/ ex ey \w?/} vD (Es) (12°7), 
subject to equation (12-5) above. The vorticity w as usually defined satisfies 
relation w = &/o. 


APPENDIXES 


Appendix A. Summary of notation and conventions. 


A tensor is expressed by the notation 
A;,'* 
in which there may be any number of suffixes. The upper or lower position of the 
sufhx is important as denoting contravariance or covariance. If the suffixes are 2 
different it is understood that a suffix may take any of the values 1, 2, 3 (in 
dimensions); the above expression therefore stands for 3° different components. 
The only case in which the suffix may be repeated is when there are two identical 
suffixes, one upper and one lower. By the swmmation convention an expression such as 
Ag” 
stands for Xg-1,2,3 Ly-1,2,3 dav'®, the signs of summation being omitted for 
venience. It follows that the suffixes a and b} in the above term are d 


Agi® = A,,** 
but A,’ is not a tensor. Apart from repeated suffixes, all terms in a tensor equatiot 
must carry the same suffixes in a similar position, for example the equation 

Ati = BYC,J + My tN, 
is legitimate. These rules extend to the partial differential coefficients in a trans- 
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rmation, the coefficient.of the variable in the numerator of a partial differential 
efficient being treated as an upper suffix, for example 
OX; OXy OX, 


raids a 
OX, Ox; OX, 


be 


jk = 


nd the summation convention still holds. 

A vector is represented when necessary by the symbol A and its tensor com- 
onents by A‘, A;. Expressions containing commas among the suffixes, such as 
in A’a, Av’, Azj,x, etc., stand for tensor derivatives as defined in §6 and 


ppendix C. - 


- Appendix B. Additional properties of epsilon tensors. 


(1) If A,’ is any tensor of the kind indicated, then 
Ag Eine he Ay Ejea ar A4 E sap = APE ae seeree (B-r). 
Proof. If a, b, care all different, both sides are zero. If, say, b = c, right hand 
ide is zero and left hand side is the sum of terms of the type 
(A; Ejnq + Ap Ejan) = 0. 
If a, 6, c are all different and equal to «, 8, y, which are fixed, i.e. they represent 
e numbers 1, 2, 3, in some order, then left hand side is equal to 
A,** Eapy + Ag Epya +. Ay’” Eyap, NOt summed, 
= A,;iEy,, summed. 
(ii) If A* and B” are arbitrary vectors 
{EijaE ne + Eisn Enea + EsseE nav} Ask Bove — Ei5, A® x Eqye BO? 
= BE, (A®* — A™) XK EgcBre wee (B-2). 
Proof. Write left hand side =C =C,+C,+C;. Since both sides vanish if 
j= j, we may puti=a, j= B where «a, f, y are fixed and all different. Hence 
C,=0 unless a=y. For C, +0, Exe $0: if b=y =a, C3=0, C,=—C,, 
Be ofc = y =a, C,=0,C,——G,,C=o;ifk=y=4, 
3 C= EyyA* Eg. BY 
att 
restricted by the condition a= y. Similarly for C, +0, C,+0, C has the above 
value restricted by the conditions b = y, c= y, respectively. Hence the total 


value of C is given by C = Ey, A™. LE ave sens 


unrestricted 
F ‘ = 4 Ex (A ma, A*s) oft. Bue, 
Appendix C. Covariant derivatives. 


“It is required to find expressions for the tensors Ai 4, Aj,q defined by the equa- 


/~ VLE pr al llr (C-1), 
ee eC a ear (C2), 
19-2 
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where A‘, A, are components of the vector A as defined in §2. Let A*= Ae 
components of the vector in the Cartesian system X,. Then 
54, = oot de 
OX 
0X, 0A, 0X 
and therefore 0A; = 9 "OX, 4, 48 © > ag aga (C-3), 
0A; 0 (0X, z) 
aed om “Ox; {oe 2 
0X, 0X, 0A, , O2X, Ox, , 
= on, ‘ ax, ax, ae Ox, ex, ax, A esesee (C 4). 
Comparing 3 and 4, we have 


3A, = lie _ 4. Oxr o2X, } dit ge ae (C5), 


OX. "OX, OxOx, 

0A; r = 
or Ai. = ee CJA. . «a eee (C-6), 
0x, OX, 
OX, Ox, 0x, 


where eo 


To express I',,“ in terms of the derivatives of the g, we have, from the definiti 

in equation (2°2), a hs a ‘a 
Ogi; OX, 0?X® OX, 02X® 
Ox, OX; Ox,Oxg Ox; Ox,0x, 


Write down the other two equations obtainable by interchanging 7k, add tw 
together and subtract the third, and we get 
Ogin . Ofer O8:;\ OX, O2X® 
ee => (28 Rt — Be) = b . 
Cer ens = Ox;  Ox,/ Ox, Omj0%, ~~ ~~ (Co 
This expression is called a Christoffel symbol of the first kind. Multiply equati 
(C7) by 


_ axe aX, 
or = Oxt* Ox, 
and we get va Pc —_ OX, or? nae 7: 
Ox, "Ox,dx, 8" 
, alt follows that Dy" = gt* Tessa 
—— dork (He -1 OLei ends “f) 
a \Ony ~ Bas eae 


Appendix D. Proof of the curl curl identity. 
This takes the form 


(curl curl A), = be (A* ,) — Ay%q 
or curl curl A = grad. div. A— V2A. 


The equations of motion of a viscous fluid in tensor notation 285 
rom equation (6:8) 
(curl A); = gigk®°A, 4, 
(curl curl A); =-piq op Eo E*'A, . 4 
= 2° Ep eA t,s,D 
= g°" (6f57 — 878!) As a0, 
Dy equation (5-3) = g° An — 2" Asn 


0 
i Ox; (A?) a A; Oras 


the order of differentiation being assumed to be indifferent. 


Appendix E. Second method of simplifying the last term in the equation of 
viscous flow. 


The following method does not involve the use of the epsilon tensors, and is 
therefore more direct than that previously used. The resulting term in general 
orthogonal coordinates is more symmetrical than in equation (10°8) but is also con- 
siderably more complicated. The method is most simply applied to the covariant 
components, the term to be evaluated being g?°u; 5, ¢- 


: OUg . 
We have in full lad = > — Pee 
7) 
and Age Aas = Vie Are — Dye" Aar- 
; Ox 
Hence, putting Aq, = Ua,> 
07 Ug , Ou, OUg , Our 
Ua,b,¢ = axrOxe ig axe Wr ant 7 8b Ave 
OL as 
ap sar Pe A eed Dag’ ge Ue mene: (E-1). 
Hence 
02 ug Og LOU, : 
Poly Go maa — g° 1." ory ee lige Ay -+ terms in u, sexe) 


; 0 : 
The first two terms may be combined by analogy with g°° (55) which may be 


written in the alternative forms - 
0 0 
ee ae — ge Dye" ise 


Ox x 
d oe eee Mf 
i on : Vg Ox? Ox? 


I — OU, eat ig Clee : 
Hence Oe Vs — (g Vg =) — 2g°° Dae an + terms 1n W,. 
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In Euclidean space the terms in w, (terms of zero order in u,) may be simplifi 
by the following artifice. It appears that the terms of zero order in equation (E-1) 
are unaltered by interchanging a and b and hence uq,5,- — Us,a,¢ contains no te 
of zero order in u,. Now 

0 be be 
8" Ua,b,0 — Fa (g°° tl», c) = 8°° (ta,v,¢ — Mr, c,a) 
ra Fi (0s ws — Uy ae + Uy a,c — % cab 

and the last two terms cancel since the order of differentiation is indifferent i 


Euclidean space. Hence 


Pls, o,0 = ae {g°u,, -} + terms of 1st and 2nd order 


bye a “ gee Vguy + terms of 1st and and order 
ax" |v g ax? 


=u, s | Te (g?¢ val + terms of 1st and 2nd order. 
Ox IV g Ax? ) 


The first and second order terms have already been evaluated and we have 
be = 12 (gevgts) —2 eT Oly Be Cees. A c 
Bane = Fa aya \E VE Gay) 78h ae aa t MA ep age Bo VED 
or, writing the second term in full, 


as af be Vo 4 — gbtgrs Ou, (2s aa OL se = “fee 0 I a c s/o 
V9 0x" (g : Ox? gs ax? \axe  Ox® xt ri (e** Vg) 


In general orthogonal coordinates in the notation already used this reduces t 
the following form 


x first component of (g°*u; 5, ¢) 


~ atbe (dw a" Ox * dyb" Oy ds c* Os 


~ 5 ee the oe 1 Ca GaU 
a |a@dx° Ox * b2 dy" Oy oe tet 
2 {i558 daU I 0a obV | x 0a dc 
a |a ox dx ° ab oy’ dx ' acds” Ox 
~ 2 eae ge thee ObV . 1 Oc Oc 

a |a@ ox’ Ox © bx" Oy aie ae 


I {re be OaU dcacaU 6 ab st 
a’ 0. 


*aUUinanas (a) +1” & ae (5) * Wimaeae(e)) ae 
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ach line of this equation corresponding to a term of equation (E-3). The sum of 
he first and second lines may also be written 


I {ise ast 0 caoU , 0 aboU 
abc |x a 0x ' dy b dy | Oz TE 


G ani7ca t.0 0 0 0 
U z ( ) (5 =. I I 0a 
4 ae a* 0x T Bay a® oy + dz \at dz 
U (da a (* da 0 (ca\ . 0a O (ab 
a*bc (dx 0x as oy Gte Oz ( ale 
It may be verified by direct substitution that these forms are equivalent to that 
ncluded in equation (10°8). 


Appendix F. Method of reducing the term A’ = {curl (wi ,u°)}" = ie? (uy ,u°) 


We have 
o = 4 {curl grad. (u?u,)} 


a gl" (UU) ac 
= fica (1? ttp, ANE 
= fyica {u® (tao oF dt ee is 
since Up a Uao = Bee by equations (5:6), 6:8), (7°12), 
4 0 = At + (815; — B{8)) WE"), 
using equation (5-3); or ; 
; Ai = (E'u’),, — (uE), 6 
og? ou 

N a Cease ¢ — ye — £6 ~_— 

ow OLE I Sel Ae OX, : OX,’ 


since the terms containing Christoffel symbols cancel. Also é , vanishes identically. 
‘Hence we have finally 


Ai = {curl (ui, ,u®)}# —= yt is éa Ou® 


Hq OKs 


ae Uo Ges 


DISCUSSION 


Prof. W. Witson. Mr Lock has succeeded in showing that the utility of the 
tensor calculus is not restricted to relativistic investigations and that it can be applied 
with advantage to* hydrodynamical and similar three-dimensional problems in 
Euclidean space. This is due to the fact that the equations in these problems are, 
for the most part, statements of the equality of different expressions for the same 
component of some vector or tensor, and these equations retain their validity when 
the coordinate system is changed. It is easy to see why the tensor calculus developed 
- in association with the theory of relativity. The most important principle of that 
theory is that the equations describing physical phenomena have the same form 
whatever may be the system of reference. I should like to congratulate Mr Lock 


on his very interesting paper. 
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Mr T. SmrrH commended tensor methods, and all algebraic methods in which 
a single symbol denotes a number of terms, to the attention of physicists, and said 
that they not only simplify the formulation of problems but also lend themselves 
to very rapid numerical calculation. The distinction between covariants and con 
travariants is important; physicists interested in a particular case, as distinct from 
a more general case, are sometimes misled by the similarity between the coefficients 
in covariant and contravariant expressions and consequently reach a wrong resu 
He thought the paper would be improved if the author would explain some points — 
in the notation, such as the use of a comma to denote differentiation. 


Capt. C. W. Hume. Is it essential to use commas to denote differentiations? 
Formulae printed in this way are rather confusing, especially when an equation 
followed by a punctuational comma. 


AUTHOR’S REPLY. A good example of the necessity for distinguishing be 
covariant and contravariant components is the case of the symbol V? applied to a 
vector and to a scalar, where, as remarked in § 6, the actual formulae, while identical 
in Cartesians, are different in curvilinear coordinates. A third different expression is 
that denoted by D (x) in equation (12-5) for the case of motion symmetrical abo 
an axis. Modifications have been made to the paper in proof to explain more 
clearly the use of the comma notation. In answer to Capt. Hume, the comma 
as denoting covariant differentiation, appears the most convenient of the notations” 
one has seen, and, if it is suggested that the comma might escape notice, it may 
be remarked that Eddington, in his original tracton relativity, used no special 
notation at all, writing 4,, indifferently for a tensor of the second order or for the 
covariant derivative of the vector A,. The ordinary differential notation cannot 
be used, since the covariant derivative is not identical with the ordinary partia 
derivative. 
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DEMONSTRATIONS 


A pH Apparatus. Demonstration given on January 24, 1930, by W. E. Doran, 
.A., The Cambridge Instrument Co., Ltd. 


A solution is either acid, alkaline or neutral, and the degree of effective acidity 
or alkalinity is entirely dependent upon the concentration of undissociated hydrogen- 
ons in the solution. This hydrogen-ion concentration, which may be represented 

pH units, can be accurately determined by immersing a hydrogen electrode in 
he solution and measuring the electric potential set up. The outfit demonstrated 
s designed for this purpose, and gives direct readings in pH units on the scale. 
t consists of a hydrogen electrode, calomel electrode, portable potentiometer, 
eflecting galvanometer and scale, together with a supply of hydrogen. The poten- 
iometer is employed both for standardising the galvanometer deflection and also 
or balancing out the calomel e.m.f., the deflection of the galvanometer then being 
due to the e.m.f., of the hydrogen electrode only. 


Soa 
Hydrogen 


C=Calomel electrode; H=Hydrogen; P=Potentiometer ; R=High resistance. 


At the beginning of the demonstration, the hydrogen electrode was immersed 
sn a solution which was neither acid nor alkaline, and the scale reading was 7, 
which is the pH value of a neutral solution. By the addition of a drop of nitric acid 
the pH reading was reduced, indicating an increase in acidity, and with a little 
sodium hydroxide the pH value was increased until at 7 it was again neutral. A 
further addition of sodium hydroxide caused the liquid to assume an alkaline state. 
In this way it was demonstrated that the acidity of the solution can easily be 
varied and controlled to almost any desired value by means of the apparatus. 

The electrode method of measuring hydrogen-ion concentration is now used 
in nearly every industrial process. As an example of its value, mention was made of 
the importance of measuring the acidity of beer during manufacture. If the beer 
is too acid the effects may be detrimental to the consumer and at least unpleasant 
to the palate; on the other hand, if the acidity is too low it is favourable to the 
growth of bacteria. The actual pH value varies in different types of beer. Among 
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other applications may be mentioned the tanning industry, sugar manufacture 
bread making, preserved food manufacture, textile dyeing—in fact, there are few 
industrial processes in which hydrogen-ion concentration does not play an im 


portant part. 


Model illustrating the Mosaic Theory of the Compound Eye, /ent by he 
Department of Zoology, University of London, King’s College. Demonstra 10 
given on February 14, 1930, by Capt. C. W. HUME. 


The eye of an insect consists of a large number of tubes (ommatidia) formed 
with small lenses at their outer ends and associated with a retina at their innet 
ends, the number of ommatidia varying from 10,000 upwards. It has been suggested 
that each lens forms an independent image, the effects of all the images being 
combined as in the case of binocular vision. According to Alternburg*, however, 
the lenses serve merely as condensers. Light from any point of the object can pass 
through that one of the ommatidia which is pointing directly at it, but is absorbed 
by the walls of the other ommatidia. An erect image is consequently formed or 
the retina by the collocation of the small light-elements admitted by different 
ommatidia. 

Altenburg’s model consists of a large number of contiguous straws, blackened 
to prevent internal reflection, and arranged between a brightly coloured and illu 
minated object and a ground-glass screen. Each straw admits light from a sma 
element of the field of view, and a faint but clear image is formed on the screen. 

In the case of the insect eye the axes of the ommatidia are convergent towards 
the retina, instead of being parallel like the straws. Thus the angular magnitude 
of the object is always the same for an insect, and the definition decreases in direct 
proportion to the distance of the object. The amount of detail given by an eye 
containing 10,000 ommatidia is equivalent to that given in television by a Nipkoy 
disc with 100 apertures, whereas the number of apertures usually employed is 
about 30. 

The illumination of the insect retina is weak, as only a very narrow pencil is 
admitted from each point on the object. The retina is probably more sensitive 
than the retina of a vertebrate, but the weakness of illumination may accour 
for the inactivity of insects on cloudy days. 


* British Journal of Experimental Biology, 4, 38 (1926). 
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REVIEWS OF BOOKS 


ambridge Five-figure Tables, by F. G. Hatt, M.A. and E. K. Riga, M.A. 
Pp. viii + 76. (London: Cambridge University Press.) 35. 6d. 


To the making of books of tables there is no end, but there is still room for a compact 
ook of five-figure tables, published at a reasonable price. The book before us, containing 
ogarithms, natural sines, circular measure, values of e” and e~, useful constants, atomic 

eights, and a list of ‘“‘useful differentials and integrals,’ is full but compact, and can be 
ommended. The type is not all that could be desired, but the page is saved from over- 
crowding by the compilers’ practice of breaking it up into short columns and, after the 
rst five-figure entry, printing,in any one column, only the last three figures. No difference 
columns are necessary. The reviewer confesses to a hankering after the type used in 
Barlow’s tables. 

The table of constants is useful and might with advantage be extended, and the list 
of differential coefficients and integrals will bring joy to the heart of many a toiling chemist. 


ASE: 


Magnetism, by E. C. STONER, Ph.D. Pp. vii + 114. (London: Methuen & Co., 
tc.) 28.64. 


Dr Stoner has produced a small but very useful work on the magnetic properties of 
materials. He deals briefly but capably with the magnetic properties of atoms, diamagnet- 
ism, paramagnetism, ferromagnetism, and the magnetic properties of the elements. 
His book is confessedly one for those to whom the subject is not new, and perhaps plunges 
in medias res a little suddenly for those who like to be dosed ‘“‘with learning put lightly, 
like powder in jam,” but, as the author hopes, it is admirably fitted for those whose desire 
it is to reach a stage at which they may carry out research on the subject. It can be recom- 
mended strongly. A.F 


Magnetic, Meteorological and Seismographic Observations made at Bombay and 
Alibag, 1924, under the direction of S. K. Banerji, D.Sc. Pp. iii + 133. (Cal- 
cutta: Government of India Publication Branch.) 115. 


The delay in the appearance of this volume is of little moment in view of the great 
advance that it marks, an advance on which India is warmly to be congratulated. The 
equipment and programme of the Bombay Observatory was largely modelled on that of 
Kew, but those whose work required hourly values of any of the elements would look in 
vain for them in the first forty-two annual volumes. The breaking down of tradition caused 
by the war led to the publication of two-hourly values of all the important elements in 
1921, though the cost of printing hourly values was then regarded as prohibitive: but 
thanks to support from the Royal Society the hourly scheme has now been adopted for 
the data of pressure, wind, temperature, relative humidity, cloud, the magnetic elements 
and microseisms. No very abnormal winds are recorded, the greatest velocity being only 


27 miles an hour. Cul we 
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Greenwich Magnetic and Meteorological Observations. Pp. D 62 + E46. (London: | 
H.M. Stationery Office, 1927.) 75. 6d. 


The electrification of railways forced the magnetic observatory to be moved from 
Greenwich to near Abinger on the northern slope of Leith Hill, and the present equipment 
includes a Schuster-Smith coil magnetometer for absolute measurements of horizontal 
force as well as an earth inductor for absolute inclination. The meteorological instruments 
still include a black-bulb radiation thermometer, in spite of the difficulties elsewhere founc 
with changes in the sensitiveness of this instrument. : G.T.W 


The Observatories’ Year Book, 1927. Pp. 430. (London: H.M. Stationery Office.) 
Logs 

This volume contains not only data of continuously recorded terrestrial magnetism 
from Lerwick and Eskdalemuir, of atmospheric electricity from Lerwick, Eskdalemuit 
and Kew, and of meteorology from Aberdeen, Eskdalemuir, Cahirciveen (Valentia) am 
Kew; but also of auroras from Lerwick, and of underground water-level, atmosphe: ic 
pollution and seismology from Kew; and there is an aerological section with upper-ail 
data derived from registering-balloons. 

With regard to the microseisms, which are presumably due to waves at sea and so hz 
about seven times the amplitude in winter that they have in summer, it would be in- 
teresting to examine whether very long waves of small amplitude could be used as a 
warning of the existence of distant storms. 

India is now playing her part in the publication of hourly values of her data: we hope 
that other portions of the Empire may see their way to make their geophysical materi: 
more readily available for general use. G.T.W 


